Can X-ray constrained Hartree-Fock wavefunctions retrieve electron correlation?

Dynamic Covalent Bonds in the Ebselen Class of Antioxidants Probed by X-ray Quantum Crystallography

Dynamic bonds are essential structural ingredients of dynamic covalent chemistry that involve reversible cleavage and formation of bonds. Herein, we explore the electronic characteristics of Se-N bonds in the organo-selenium antioxidant ebselen and its derivatives for their propensity to function as dynamic covalent bonds by employing high-resolution X-ray quantum crystallography and complementary computational studies. An analysis of the experimentally reconstructed X-ray wavefunctions reveals the salient electronic features of the Se-N bonds with very low electron density localized at the bonding region and a positive Laplacian value at the bond critical point. Bond orders and percentage covalency and ionicity estimated from the X-ray wavefunctions, along with localized orbital locator (LOL) and electron localization function (ELF) analyses show that the Se-N bond is unique in its closed shell-like features, despite being a covalent bond. Time-dependent DFT calculations simulate the cleavage of Se-N bonds in ebselen in the excited state, further substantiating their nature as dynamic bonds.

Introduction of a weighting scheme for the X-ray restrained wavefunction approach: advantages and drawbacks

In a quite recent study [Genoni et al. (2017). IUCrJ, 4, 136-146], it was observed that the X-ray restrained wavefunction (XRW) approach allows a more efficient and larger capture of electron correlation effects on the electron density if high-angle reflections are not considered in the calculations. This is due to the occurrence of two concomitant effects when one uses theoretical X-ray diffraction data corresponding to a single-molecule electron density in a large unit cell: (i) the high-angle reflections are generally much more numerous than the low- and medium-angle ones, and (ii) they are already very well described at unrestrained level. Nevertheless, since high-angle data also contain important information that should not be disregarded, it is not advisable to neglect them completely. For this reason, based on the results of the previous investigation, this work introduces a weighting scheme for XRW calculations to up-weight the contribution of low- and medium-angle reflections, and, at the same time, to reasonably down-weight the importance of the high-angle data. The proposed strategy was tested through XRW computations with both theoretical and experimental structure-factor amplitudes. The tests have shown that the new weighting scheme works optimally if it is applied with theoretically generated X-ray diffraction data, while it is not advantageous when traditional experimental X-ray diffraction data (even of very high resolution) are employed. This also led to the conclusion that the use of a specific external parameter λ_J for each resolution range might not be a suitable strategy to adopt in XRW calculations exploiting experimental X-ray data as restraints.

X-ray constrained wavefunctions based on Hirshfeld atoms. II. Reproducibility of electron densities in crystals of α-oxalic acid dihydrate

The Hirshfeld atom-based X-ray constrained wavefunction fitting (HA-XCW) procedure is tested for its reproducibility, and the information content of the fitted wavefunction is critically assessed. Fourteen different α-oxalic acid dihydrate data sets are used for this purpose, and the first joint fitting to 12 of these data sets is reported. There are systematic features in the electron density obtained from all data sets which agree with higher level benchmark calculations, but there are also many other strong systematic features which disagree with the reference calculations, most notably those associated with the electron density near the nuclei. To enhance reproducibility, three new protocols are described and tested to address the halting problem of XCW fitting, namely: an empirical power-function method, which is useful for estimating the accuracy of the structure factor uncertainties; an asymptotic extrapolation method based on ideas from density functional theory; and a `conservative method' whereby the smallest value of the regularization parameter is chosen from a series of data sets, or subsets.

Remarks on X-ray constrained/restrained wavefunction fitting

X-ray constrained/restrained wavefunctions (XCWs/XRWs) result from a combination of theory and experiment and are therefore affected by experimental errors and model uncertainties. The present XCW/XRW procedure does not take this into account, thus limiting the meaning and significance of the obtained wavefunctions.

X-ray constrained wavefunctions based on Hirshfeld atoms. I. Method and review

The X-ray constrained wavefunction (XCW) procedure for obtaining an experimentally reconstructed wavefunction from X-ray diffraction data is reviewed. The two-center probability distribution model used to perform nuclear-position averaging in the original paper [Grimwood & Jayatilaka (2001). Acta Cryst. A57, 87-100] is carefully distinguished from the newer one-center probability distribution model. In the one-center model, Hirshfeld atoms are used, and the Hirshfeld atom based X-ray constrained wavefunction (HA-XCW) procedure is described for the first time, as well as its efficient implementation. In this context, the definition of the related X-ray wavefunction refinement (XWR) method is refined. The key halting problem for the XCW method - the procedure by which one determines when overfitting has occurred - is named and work on it reviewed.

π-Hole bonding in a new co-crystal hydrate of gallic acid and pyrazine: static and dynamic charge density analysis

A new cocrystal hydrate of gallic acid with pyrazine (4GA, Py, 4H₂O; GA₄PyW₄) was obtained and characterized by single crystal X-ray diffraction. In addition to structure determination, experimental charge density analysis was carried out in terms of Multipole Modelling (MP), X-ray wavefunction refinement (XWR) and maximum entropy method (MEM). As a part of XWR, the structural refinement via Hirshfeld atom refinement was carried out and resulted in O-H bond lengths close to values from neutron diffraction. A systematic comparison of molecular conformations and aromatic interactions in this new cocrystal hydrate was performed with other existing polymorphs of gallic acid. In GA₄PyW₄, the two symmetry-independent gallic acid molecules have a syn COOH orientation and form the common (COOH)₂ dimeric synthon. The carboxyl C atom displays the characteristics of π-holes with electropositive regions above and below the molecular plane and engages in acceptor-donor interactions with oxygen atoms of acidic O-H groups and phenol groups of neighbouring gallic acid molecules. The signature of the π-hole was identified from experimental charge density analysis, both in static density maps in MP and XWR as well as dynamic density in MEM, but it cannot be pinned down to a specific atom-atom interaction. This study presents the first comparison between an XWR and a MEM experimental electron-density determination.

Structure Factors and Charge Density Description of Aluminum: A Quantum Crystallographic Study

Static structure factors and charge density for metallic aluminum were investigated by periodic calculations using atom-centered Gaussian-type basis sets with the Perdew-Burke-Ernzerhof (PBE) functional implemented in the CRYSTAL14 package and X-ray constrained wave function (XCW) fitting. The effects of additional diffuse d and f basis functions on structure factors were compared with synchrotron powder X-ray diffraction and quantitative convergent electron beam diffraction data. Changes in structure factors from an independent atom model at 022, 113, and 222 reflections introduced d and f basis functions similar to those of the experimental data. The XCW fitting was applied to different sizes of aluminum clusters. The charge density features for a 50-atom cluster clearly demonstrated electron accumulation at tetrahedral sites and electron depletion at octahedral sites. The resolution dependence of the XCW study suggests that structure factors of the five lowest resolution reflections with 0.1% accuracy were indispensable for determining the detailed bonding description in the case of metallic aluminum.

New venues in electron density analysis

We provide a comprehensive overview of the chemical information within the electron density: how to extract information, but also how to obtain and how to assess the quality of the...

Charge density studies of multicentre two-electron bonding of an anion radical at non-ambient temperature and pressure

The variation of charge density of two-electron multicentre bonding (pancake bonding) between semi-quinone radicals with pressure and temperature was studied on a salt of 5,6-di-chloro-2,3-di-cyano-semi-quinone radical anion (DDQ) with 4-cyano-N-methyl-pyridinium cation (4-CN) using the Transferable Aspheric Atom Model (TAAM) refinement. The pancake-bonded radical dimers are stacked by non-bonding π-interactions. With rising pressure, the covalent character of interactions between radicals increases, and above 2.55 GPa, the electron density indicates multicentric covalent interactions throughout the stack. The experimental charge densities were verified and corroborated by periodic DFT computations. The TAAM approach has been tested and validated for atomic resolution data measured at ambient pressure; this work shows this approach can also be applied to diffraction data obtained at pressures up to several gigapascals.

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.

Measuring Density Functional Parameters from Electron Diffraction Patterns

We integrate density functional theory (DFT) into quantitative convergent-beam electron diffraction (QCBED) to create a synergy between experiment and theory called QCBED-DFT. This synergy resides entirely in the electron density which, in real materials, gives rise to the experimental CBED patterns used by QCBED-DFT to refine DFT model parameters. We use it to measure the Hubbard energy U for two strongly correlated electron systems, NiO and CeB_{6} (U=7.4±0.6  eV for d orbitals in NiO and U=3.0±0.6  eV for f orbitals in CeB_{6}), and the boron position parameter x for CeB_{6} (x=0.1992±0.0003). In verifying our measurements, we demonstrate an accuracy test for any modeled electron density.

Similarities and Differences between Crystal and Enzyme Environmental Effects on the Electron Density of Drug Molecules

The crystal interaction density is generally assumed to be a suitable measure of the polarization of a low-molecular weight ligand inside an enzyme, but this approximation has seldomly been tested and has never been quantified before. In this study, we compare the crystal interaction density and the interaction electrostatic potential for a model compound of loxistatin acid (E64c) with those inside cathepsin B, in solution, and in vacuum. We apply QM/MM calculations and experimental quantum crystallography to show that the crystal interaction density is indeed very similar to the enzyme interaction density. Less than 0.1 e are shifted between these two environments in total. However, this difference has non-negligible consequences for derived properties.

HgH2 meets relativistic quantum crystallography. How to teach relativity to a non-relativistic wavefunction

The capability of X-ray constrained wavefunction (XCW) fitting to introduce relativistic effects into a non-relativistic wavefunction is tested. It is quantified how much of the reference relativistic effects can be absorbed in the non-relativistic XCW calculation when fitted against relativistic structure factors of a model HgH₂ molecule. Scaling of the structure-factor sets to improve the agreement statistics is found to introduce a significant systematic error into the XCW fitting of relativistic effects.

Quantum Crystallography in the Last Decade: Developments and Outlooks

In this review article, we report on the recent progresses in the field of quantum crystallography that has witnessed a massive increase of production coupled with a broadening of the scope in the last decade. It is shown that the early thoughts about extracting quantum mechanical information from crystallographic experiments are becoming reality, although a century after prediction. While in the past the focus was mainly on electron density and related quantities, the attention is now shifting toward determination of wavefunction from experiments, which enables an exhaustive determination of the quantum mechanical functions and properties of a system. Nonetheless, methods based on electron density modelling have evolved and are nowadays able to reconstruct tiny polarizations of core electrons, coupling charge and spin models, or determining the quantum behaviour at extreme conditions. Far from being routine, these experimental and computational results should be regarded with special attention by scientists for the wealth of information on a system that they actually contain.

TAAM: a reliable and user friendly tool for hydrogen-atom location using routine X-ray diffraction data

Hydrogen is present in almost all of the molecules in living things. It is very reactive and forms bonds with most of the elements, terminating their valences and enhancing their chemistry. X-ray diffraction is the most common method for structure determination. It depends on scattering of X-rays from electron density, which means the single electron of hydrogen is difficult to detect. Generally, neutron diffraction data are used to determine the accurate position of hydrogen atoms. However, the requirement for good quality single crystals, costly maintenance and the limited number of neutron diffraction facilities means that these kind of results are rarely available. Here it is shown that the use of Transferable Aspherical Atom Model (TAAM) instead of Independent Atom Model (IAM) in routine structure refinement with X-ray data is another possible solution which largely improves the precision and accuracy of X-H bond lengths and makes them comparable to averaged neutron bond lengths. TAAM, built from a pseudoatom databank, was used to determine the X-H bond lengths on 75 data sets for organic molecule crystals. TAAM parametrizations available in the modified University of Buffalo Databank (UBDB) of pseudoatoms applied through the DiSCaMB software library were used. The averaged bond lengths determined by TAAM refinements with X-ray diffraction data of atomic resolution (d_min ≤ 0.83 Å) showed very good agreement with neutron data, mostly within one single sample standard deviation, much like Hirshfeld atom refinement (HAR). Atomic displacements for both hydrogen and non-hydrogen atoms obtained from the refinements systematically differed from IAM results. Overall TAAM gave better fits to experimental data of standard resolution compared to IAM. The research was accompanied with development of software aimed at providing user-friendly tools to use aspherical atom models in refinement of organic molecules at speeds comparable to routine refinements based on spherical atom model.

On the use of the Obara–Saika recurrence relations for the calculation of structure factors in quantum crystallography

Modern methods of quantum crystallography are techniques firmly rooted in quantum chemistry and, as in many quantum chemical strategies, electron densities are expressed as two-centre expansions that involve basis functions centred on atomic nuclei. Therefore, the computation of the necessary structure factors requires the evaluation of Fourier transform integrals of basis function products. Since these functions are usually Cartesian Gaussians, in this communication it is shown that the Fourier integrals can be efficiently calculated by exploiting an extension of the Obara–Saika recurrence formulas, which are successfully used by quantum chemists in the computation of molecular integrals. Implementation and future perspectives of the technique are also discussed.

Hydrogen atoms in bridging positions from quantum crystallographic refinements: influence of hydrogen atom displacement parameters on geometry and electron density

Hydrogen atom positions can be obtained accurately from X-ray diffraction data of hydrogen maleate salts via Hirshfeld atom refinement.

X-ray constrained spin-coupled technique: theoretical details and further assessment of the method

One of the well-established methods of modern quantum crystallography is undoubtedly the X-ray constrained wavefunction (XCW) approach, a technique that enables the determination of wavefunctions which not only minimize the energy of the system under examination, but also reproduce experimental X-ray diffraction data within the limit of the experimental errors. Initially proposed in the framework of the Hartree–Fock method, the strategy has been gradually extended to other techniques of quantum chemistry, but always remaining limited to a single-determinant ansatz for the wavefunction to extract. This limitation has been recently overcome through the development of the novel X-ray constrained spin-coupled (XCSC) approach [Genoni et al. (2018). Chem. Eur. J. 24, 15507–15511] which merges the XCW philosophy with the traditional spin-coupled strategy of valence bond theory. The main advantage of this new technique is the possibility of extracting traditional chemical descriptors (e.g. resonance structure weights) compatible with the experimental diffraction measurements, without the need to introduce information a priori or perform analyses a posteriori. This paper provides a detailed theoretical derivation of the fundamental equations at the basis of the XCSC method and also introduces a further advancement of its original version, mainly consisting in the use of molecular orbitals resulting from XCW calculations at the Hartree–Fock level to describe the inactive electrons in the XCSC computations. Furthermore, extensive test calculations, which have been performed by exploiting high-resolution X-ray diffraction data for salicylic acid and by adopting different basis sets, are presented and discussed. The computational tests have shown that the new technique does not suffer from particular convergence problems. Moreover, all the XCSC calculations provided resonance structure weights, spin-coupled orbitals and global electron densities slightly different from those resulting from the corresponding unconstrained computations. These discrepancies can be ascribed to the capability of the novel strategy to capture the information intrinsically contained in the experimental data used as external constraints.

Relativistic quantum crystallography of diphenyl- and dicyanomercury. Theoretical structure factors and Hirshfeld atom refinement

Quantum crystallographic refinement of heavy-element-containing compounds is a challenge, because many physical effects have to be accounted for adequately. Here, the impact and magnitude of relativistic effects are compared with those of electron correlation, polarization through the environment, choice of basis set and treatment of thermal motion effects on the structure factors of diphenylmercury(II) [Hg(Ph)₂] and dicyanomercury(II) [Hg(CN)₂]. Furthermore, the individual atomic contributions to the structure factors are explored in detail (using Mulliken population analysis and the exponential decay of atomic displacement parameters) to compare the contributions of lighter atoms, especially hydrogen atoms, against mercury. Subsequently, relativistic Hirshfeld atom refinement (HAR) is validated against theoretical structure factors of Hg(Ph)₂ and Hg(CN)₂, starting from perturbed geometries, to test if the relativistic variant of HAR leads to multiple solutions. Generally, relativistic HAR is successful, leading to a perfect match with the reference geometries, but some limitations are pointed out.

Revisiting a Historical Concept by Using Quantum Crystallography: Are Phosphate, Sulfate and Perchlorate Anions Hypervalent?

There are many examples of atoms in molecules that violate Lewis' octet rule, because they have more than four electron pairs assigned to their valence. These atoms are referred to as hypervalent. However, hypervalency may be regarded as an artifact arising from Lewis' description of molecules, which is based on the assumption that electrons are localized in two-center two-electron bonds and lone pairs. In the present paper, the isoelectronic phosphate (PO₄ ^3- ), sulfate (SO₄ ^2- ) and perchlorate (ClO₄ ^- ) anions were examined with respect to the concept of hypervalency. Lewis formulas containing a hypervalent central atom exist for all three anions. Based on X-ray wavefunction refinements of high-resolution X-ray diffraction data of representative crystal structures (MgNH₄ PO₄ ⋅6 H₂ O, Li₂ SO₄ ⋅H₂ O, and KClO₄ ), complementary bonding analyses were performed. In this way, experimental information from the new field of quantum crystallography validate long-known facts, or refute long-standing misunderstandings. It is shown that the P-O and S-O bonds are highly polarized covalent bonds and, thus, the increase in the valence population following three-center four-electron bonding is not sufficient to yield hypervalent phosphorus or sulfur atoms, respectively. However, for the highly covalent Cl-O bond, most bonding indicators imply a hypervalent chlorine atom.

Libraries of Extremely Localized Molecular Orbitals. 3. Construction and Preliminary Assessment of the New Databanks

The fast and reliable determination of wave functions and electron densities of macromolecules has been one of the goals of theoretical chemistry for a long time, and in this context, several linear scaling techniques have been successfully devised over the years. Different approaches have been adopted to tackle this problem, and one of them exploits the fact that, according to the traditional chemical perception, molecules can be seen as constituted of recurring units (e.g., functional groups) with well-defined chemical features. This has led to the development of methods in which the global wave functions or electron densities of macromolecules are obtained by simply transferring density matrices or fuzzy electron densities associated with molecular fragments. In this context, we propose an alternative strategy that aims at quickly reconstructing wave functions and electron densities of proteins through the transfer of extremely localized molecular orbitals (ELMOs), which are orbitals strictly localized on small molecular units and, for this reason, easily transferable from molecule to molecule. To accomplish this task we have constructed original libraries of ELMOs that cover all the possible elementary fragments of the 20 natural amino acids in all their possible protonation states and forms. Our preliminary test calculations have shown that, compared to more traditional methods of quantum chemistry, the transfers from the novel ELMO databanks allow to obtain wave function and electron densities of large polypeptides and proteins at a significantly reduced computational cost. Furthermore, notwithstanding expected discrepancies, the obtained electron distributions and electrostatic potentials are in very good agreement with those obtained at Hartree-Fock and density functional theory (DFT) levels. Therefore, the results encourage to use the new libraries as alternatives to the popular pseudoatom-databases of crystallography in the refinement of crystallographic structures of macromolecules. In particular, in this context, we have already envisaged the coupling of the ELMO databanks with the promising Hirshfeld atom refinement technique to extend the applicability of the latter to very large systems.

Quantum Crystallography: Current Developments and Future Perspectives

Crystallography and quantum mechanics have always been tightly connected because reliable quantum mechanical models are needed to determine crystal structures. Due to this natural synergy, nowadays accurate distributions of electrons in space can be obtained from diffraction and scattering experiments. In the original definition of quantum crystallography (QCr) given by Massa, Karle and Huang, direct extraction of wavefunctions or density matrices from measured intensities of reflections or, conversely, ad hoc quantum mechanical calculations to enhance the accuracy of the crystallographic refinement are implicated. Nevertheless, many other active and emerging research areas involving quantum mechanics and scattering experiments are not covered by the original definition although they enable to observe and explain quantum phenomena as accurately and successfully as the original strategies. Therefore, we give an overview over current research that is related to a broader notion of QCr, and discuss options how QCr can evolve to become a complete and independent domain of natural sciences. The goal of this paper is to initiate discussions around QCr, but not to find a final definition of the field.

Probing the accuracy and precision of Hirshfeld atom refinement with HARt interfaced with Olex2

Hirshfeld atom refinement (HAR) is a novel X-ray structure refinement technique that employs aspherical atomic scattering factors obtained from stockholder partitioning of a theoretically determined tailor-made static electron density. HAR overcomes many of the known limitations of independent atom modelling (IAM), such as too short element-hydrogen distances, r(X-H), or too large atomic displacement parameters (ADPs). This study probes the accuracy and precision of anisotropic hydrogen and non-hydrogen ADPs and of r(X-H) values obtained from HAR. These quantities are compared and found to agree with those obtained from (i) accurate neutron diffraction data measured at the same temperatures as the X-ray data and (ii) multipole modelling (MM), an established alternative method for interpreting X-ray diffraction data with the help of aspherical atomic scattering factors. Results are presented for three chemically different systems: the aromatic hydro-carbon rubrene (orthorhombic 5,6,11,12-tetra-phenyl-tetracene), a co-crystal of zwitterionic betaine, imidazolium cations and picrate anions (BIPa), and the salt potassium hydrogen oxalate (KHOx). The non-hydrogen HAR-ADPs are as accurate and precise as the MM-ADPs. Both show excellent agreement with the neutron-based values and are superior to IAM-ADPs. The anisotropic hydrogen HAR-ADPs show a somewhat larger deviation from neutron-based values than the hydrogen SHADE-ADPs used in MM. Element-hydrogen bond lengths from HAR are in excellent agreement with those obtained from neutron diffraction experiments, although they are somewhat less precise. The residual density contour maps after HAR show fewer features than those after MM. Calculating the static electron density with the def2-TZVP basis set instead of the simpler def2-SVP one does not improve the refinement results significantly. All HARs were performed within the recently introduced HARt option implemented in the Olex2 program. They are easily launched inside its graphical user interface following a conventional IAM.

Validation of X-ray Wavefunction Refinement

In this work, the quality of the electron density in crystals reconstructed by the multipolar model (MM) and by X-ray wavefunction refinement (XWR) is tested on a set of high-resolution X-ray diffraction data sets of four amino acids and six tripeptides. It results in the first thorough validation of XWR. Agreement statistics, figures of merit, residual- and deformation-density maps, as well as atomic displacement parameters are used to measure the quality of the reconstruction relative to the measured structure factors. Topological analysis of the reconstructed density is carried out to obtain atomic and bond-topological properties, which are subsequently compared to the values derived from benchmarking periodic DFT geometry optimizations. XWR is simultaneously in better agreement than the MM with both benchmarking theory and the measured diffraction pattern. In particular, the obvious problems with the description of polar bonds in the MM are significantly reduced by using XWR. Similarly, modeling of electron density in the vicinity of hydrogen atoms with XWR is visibly improved.

Exploring charge density analysis in crystals at high pressure: data collection, data analysis and advanced modelling

The possibility to determine electron-density distribution in crystals has been an enormous breakthrough, stimulated by a favourable combination of equipment for X-ray and neutron diffraction at low temperature, by the development of simplified, though accurate, electron-density models refined from the experimental data and by the progress in charge density analysis often in combination with theoretical work. Many years after the first successful charge density determination and analysis, scientists face new challenges, for example: (i) determination of the finer details of the electron-density distribution in the atomic cores, (ii) simultaneous refinement of electron charge and spin density or (iii) measuring crystals under perturbation. In this context, the possibility of obtaining experimental charge density at high pressure has recently been demonstrated [Casatiet al.(2016).Nat. Commun.7, 10901]. This paper reports on the necessities and pitfalls of this new challenge, focusing on the speciessyn-1,6:8,13-biscarbonyl[14]annulene. The experimental requirements, the expected data quality and data corrections are discussed in detail, including warnings about possible shortcomings. At the same time, new modelling techniques are proposed, which could enable specific information to be extracted, from the limited and less accurate observations, like the degree of localization of double bonds, which is fundamental to the scientific case under examination.

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?

Quantum crystallography

Approximate wavefunctions can be improved by constraining them to reproduce observations derived from diffraction and scattering experiments. Conversely, charge density models, incorporating electron-density distributions, atomic positions and atomic motion, can be improved by supplementing diffraction experiments with quantum chemically calculated, tailor-made electron densities (form factors). In both cases quantum chemistry and diffraction/scattering experiments are combined into a single, integrated tool. The development of quantum crystallographic research is reviewed. Some results obtained by quantum crystallography illustrate the potential and limitations of this field.

A first-prototype multi-determinant X-ray constrained wavefunction approach: the X-ray constrained extremely localized molecular orbital–valence bond method

All the current variants of Jayatilaka's X-ray constrained wavefunction (XCW) approach work within the framework of the single-determinant wavefunctionansatz. In this paper, a first-prototype multi-determinant XCW technique is proposed. The strategy assumes that the desired XCW is written as a valence-bond-like expansion in terms of pre-determined single Slater determinants constructed with extremely localized molecular orbitals. The method, which can be particularly suitable to investigate systems with a multi-reference character, has been applied to determine the weights of the resonance structures of naphthalene at different temperatures by exploiting experimental high-resolution X-ray diffraction data. The results obtained have shown that the explicit consideration of experimental structure factors in the determination of the resonance structure weights may lead to results significantly different compared with those resulting only from the simple energy minimization.

Developments in the structural science of materials

Recent developments in the structural science of materials and the growing power of computational methods in this field are discussed.