X-ray Constrained Wavefunction Research Articles

The effects of experimentally obtained electron correlation and polarization on electron densities and exchange-correlation potentials.

In X-ray constrained wavefunction (XCW) fitting, external information, such as electron correlation and polarization, is included into a single-determinantal isolated-molecule wavefunction. In a first step, we show that the extraction of these two physical effects by XCW fitting is complete and accurate by comparing to theoretical reference calculations. In a second step, we show that fitting to data from single-crystal x-ray diffraction measurements provides the same results qualitatively and how the physical effects can be separated, although always inherently convolved in the experiment. We further demonstrate that exchange-correlation potentials are systematically affected by XCW fitting in a physically meaningful way, which could be exploited for method development in quantum chemistry, subject to some remaining challenges that we also outline.

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.

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.

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 HgH2 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.

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.

X-ray Constrained Spin-Coupled Wavefunction: a New Tool to Extract Chemical Information from X-ray Diffraction Data.

The X-ray constrained wavefunction (XCW) approach is a reliable and widely used method of quantum crystallography that allows the determination of wavefunctions compatible with X-ray diffraction data. So far, all the existing XCW techniques have been developed in the framework of molecular orbital theory and, consequently, provide only pictures of the "experimental" electronic structures that are far from the traditional chemical perception. Here a new strategy is proposed that, by combining the XCW philosophy with the spin-coupled method of valence bond theory, enables direct extraction of traditional chemical information (e.g., weights of resonance structures) from X-ray diffraction measurements. Preliminary results have shown that the new technique is really able to efficiently capture the effects of the crystal environment on the electronic structure, and can be considered as a new useful tool to perform chemically sound analyses of the X-ray diffraction data.

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 wavefunction ansatz. 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.

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

The X-ray constrained wavefunction (XC-WF) method proposed by Jayatilaka [Jayatilaka & Grimwood (2001) ▸, Acta Cryst. A57, 76-86] has attracted much attention because it represents a possible third way of theoretically studying the electronic structure of atoms and molecules, combining features of the more popular wavefunction- and DFT-based approaches. In its original formulation, the XC-WF technique extracts statistically plausible wavefunctions from experimental X-ray diffraction data of molecular crystals. A weight is used to constrain the pure Hartree-Fock solution to the observed X-ray structure factors. Despite the wavefunction being a single Slater determinant, it is generally assumed that its flexibility could guarantee the capture, better than any other experimental model, of electron correlation effects, absent in the Hartree-Fock Hamiltonian but present in the structure factors measured experimentally. However, although the approach has been known for long time, careful testing of this fundamental hypothesis is still missing. Since a formal demonstration is impossible, the validation can only be done heuristically and, to accomplish this task, X-ray constrained Hartree-Fock calculations have been performed using structure factor amplitudes computed at a very high correlation level (coupled cluster) for selected molecules in isolation, in order to avoid the perturbations due to intermolecular interactions. The results show that a single-determinant XC-WF is able to capture the electron correlation effects only partially. The largest amount of electron correlation is extracted when: (i) a large external weight is used (much larger than what has normally been used in XC-WF calculations using experimental data); and (ii) the high-order reflections, which carry less information on the electron correlation, are down-weighted (or even excluded), otherwise they would bias the fitting towards the unconstrained Hartree-Fock wavefunction.

Importance of Relativistic Effects and Electron Correlation in Structure Factors and Electron Density of Diphenyl Mercury and Triphenyl Bismuth.

This study investigates the possibility of detecting relativistic effects and electron correlation in single-crystal X-ray diffraction experiments using the examples of diphenyl mercury (HgPh2) and triphenyl bismuth (BiPh3). In detail, the importance of electron correlation (ECORR), relativistic effects (REL) [distinguishing between total, scalar and spin-orbit (SO) coupling relativistic effects] and picture change error (PCE) on the theoretical electron density, its topology and its Laplacian using infinite order two component (IOTC) wave functions is discussed. This is to develop an understanding of the order of magnitude and shape of these different effects as they manifest in the electron density. Subsequently, the same effects are considered for the theoretical structure factors. It becomes clear that SO and PCE are negligible, but ECORR and scalar REL are important in low- and medium-order reflections on absolute and relative scales-not in the high-order region. As a further step, Hirshfeld atom refinement (HAR) and subsequent X-ray constrained wavefunction (XCW) fitting have been performed for the compound HgPh2 with various relativistic and nonrelativistic wave functions against the experimental structure factors. IOTC calculations of theoretical structure factors and relativistic HAR as well as relativistic XCW fitting are presented for the first time, accounting for both scalar and spin-orbit relativistic effects.