Bonding features in Appel's salt from the orbital-free quantum crystallographic perspective

Bonding properties in the crystal of 4,5-dichloro-l,2,3-dithiazolium chloride (Appel's salt) were studied using a combination of single-crystal high-resolution X-ray diffraction data and the orbital-free quantum crystallography approach. A QTAIM-based topological model shows the proximity of S—C and S—N bonds to the sesquialteral type and establishes the low S—S bond order in the l,2,3-dithiazolium heterocycle. It is found that the electrostatic potential carries the traces of a common positive area on the junction of interatomic zero-flux surfaces of S1 and S2 atomic basins; meanwhile the exchange energy density per particle shows perfectly here two separate minima through which the two bond paths run. Thus, the pair intermolecular interactions Cl−...S1 and Cl−...S2 formed by the common chloride anion placed near the center of the S—S bond are categorized as chalcogen bonds.

Developing orbital-free quantum crystallography: the local potentials and associated partial charge densities

This work extends the orbital-free density functional theory to the field of quantum crystallography. The total electronic energy is decomposed into electrostatic, exchange, Weizsacker and Pauli components on the basis of physically grounded arguments. Then, the one-electron Euler equation is re-written through corresponding potentials, which have clear physical and chemical meaning. Partial electron densities related with these potentials by the Poisson equation are also defined. All these functions were analyzed from viewpoint of their physical content and limits of applicability. Then, they were expressed in terms of experimental electron density and its derivatives using the orbital-free density functional theory approximations, and applied to the study of chemical bonding in a heteromolecular crystal of ammonium hydrooxalate oxalic acid dihydrate. It is demonstrated that this approach allows the electron density to be decomposed into physically meaningful components associated with electrostatics, exchange, and spin-independent wave properties of electrons or with their combinations in a crystal. Therefore, the bonding information about a crystal that was previously unavailable for X-ray diffraction analysis can be now obtained.

Detailed analysis of deformation potentials with application in orbital-free density functional theory

A detailed analysis of the recently published deformation potentials for application in orbital-free density functional theory is given. Since orbital-free density functional theory is a purely density-based description of quantum mechanics, it may in the future provide itself useful in quantum crystallography as it establishes a direct link between experiment and theory via a single meaningful quantity: the electron density. In order to establish this goal, sufficiently accurate approximations for the kinetic energy have to be found. The present work is a further step in this direction. The so-called deformation potentials allow the interaction between the atoms to be taken into account through the help of their electron density only. It is shown that the present ansatz provides a systematic pathway beyond the recently introduced atomic fragment approach.

Further Validation of Quantum Crystallography Approaches

Quantum crystallography is a fast-developing multidisciplinary area of crystallography. In this work, we analyse the influence of different charge density models (i.e., the multipole model (MM), Hirshfeld atom refinement (HAR), and the transferable aspherical atom model (TAAM)), modelling of the thermal motion of hydrogen atoms (anisotropic, isotropic, and with the aid of SHADE or NoMoRe), and the type of radiation used (Mo Kα and Cu Kα) on the final results. To achieve this aim, we performed a series of refinements against X-ray diffraction data for three model compounds and compared their final structures, geometries, shapes of ADPs, and charge density distributions. Our results were also supported by theoretical calculations that enabled comparisons of the lattice energies of these structures. It appears that geometrical parameters are better described (closer to the neutron values) when HAR is used; however, bonds to H atoms more closely match neutron values after MM or TAAM refinement. Our analysis shows the superiority of the NoMoRe method in the description of H-atom ADPs. Moreover, the shapes of the ADPs of H atoms, as well as their electron density distributions, were better described with low-resolution Cu Kα data in comparison to low-resolution Mo Kα data.

Towards Mapping Chemical Bonds in Nanostructured Materials

We introduce a number of techniques in quantitative convergent-beam electron diffraction under development by our group and discuss the basis for measuring interatomic electrostatic potentials (and therefore also electron densities), localised at sub-nanometre scales, with sufficient accuracy and precision to map chemical bonds in and around nanostructures in nanostructured materials. This has never before been possible as experimental measurements of bonding in quantum crystallography have hitherto always been restricted to homogeneous single-phased crystals.

lamaGOET: an interface for quantum crystallography

In quantum crystallography, theoretical calculations and crystallographic refinements are closely intertwined. This means that the employed software must be able to perform both quantum-mechanical calculations and crystallographic least-squares refinements. So far, the program Tonto is the only one able to do that. The lamaGOET interface described herein deals with this issue since it interfaces dedicated quantum-chemical software (the widely used Gaussian package and the specialized ELMOdb program) with the refinement capabilities of Tonto. Three different flavours of quantum-crystallographic refinements of the dipetide glycyl-L-threonine dihydrate are presented to showcase the capabilities of lamaGOET: Hirshfeld atom refinement (HAR), HAR-ELMO, namely HAR coupled with extremely localized molecular orbitals, and X-ray constrained wavefunction fitting.

The advanced treatment of hydrogen bonding in quantum crystallography

Although hydrogen bonding is one of the most important motifs in chemistry and biology, H-atom parameters are especially problematic to refine against X-ray diffraction data. New developments in quantum crystallography offer a remedy. This article reports how hydrogen bonds are treated in three different quantum-crystallographic methods: Hirshfeld atom refinement (HAR), HAR coupled to extremely localized molecular orbitals and X-ray wavefunction refinement. Three different compound classes that form strong intra- or intermolecular hydrogen bonds are used as test cases: hydrogen maleates, the tripeptide L-alanyl-glycyl-L-alanine co-crystallized with water, and xylitol. The differences in the quantum-mechanical electron densities underlying all the used methods are analysed, as well as how these differences impact on the refinement results.