Topology of the electron density of d0 transition metal compounds at subatomic resolution

The Effects of Chemical Bonding at Subatomic Resolution: A Case Study on α-Boron

Similar to classical asphericity shifts, aspherical deformations of the electron density in the atomic core region can result in core asphericity shifts in refinements using a Hansen-Coppens multipolar model (HCM), especially when highly precise experimental datasets with resolutions far beyond sin(θ)/λ ≤ 1.0 Å⁻¹ are employed. These shifts are about two orders of magnitude smaller than their counterparts caused by valence shell deformations, and their underlying deformations are mainly of dipolar character for 1st row atoms. Here, we analyze the resolution dependence of core asphericity shifts in α-boron. Based on theoretical structure factors, an appropriate Extended HCM (EHCM) is developed, which is tested against experimental high-resolution (sin(θ)/λ ≤ 1.6 Å⁻¹) single-crystal diffraction data. Bond length deviations due to core asphericity shifts of α-boron in the order of 4–6·10⁻⁴ Å are small but significant at this resolution and can be effectively compensated by an EHCM, although the correlation of the additional model parameters with positional parameters prevented a free refinement of all core model parameters. For high quality, high resolution data, a proper treatment with an EHCM or other equivalent methods is therefore highly recommended.

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.

Bonding and metastability for Group 12 dications

Electronic structure and bonding properties of the Group 12 dications M₂ ²⁺ (M = Zn, Cd, Hg) are investigated and electron density-derived quantities are used to characterize the metastability of these species. Of particular interest are the complementary descriptions afforded by the Laplacian of the electron density ∇² ρ(r) and the one-electron Bohm quantum potential (Q =  ∇ 2 ρ r / 2 ρ r ) along the bond path. Further, properties derived from the pair density including the localization-delocalization matrices (LDMs) and the interacting quantum atoms (IQA) energies are analyzed within the framework of the quantum theory of atoms in molecules (QTAIM). From the crossing points of the singlet (ground) and triplet (excited) potential energy curves, the barriers for dissociation (BFD) are estimated to be 25.2 kcal/mol (1.09 eV) for Zn₂ ²⁺ , 22.8 kcal/mol (0.99 eV) for Cd₂ ²⁺ , and 26.4 kcal/mol (1.14 eV) for Hg₂ ²⁺ . For comparison and benchmarking purposes, the case of N₂ ²⁺ is considered as a texbook example of metastability. At the equilibrium geometries, LDMs, which are used here as an electronic fingerprinting tool, discriminate and group together Group 12 M₂ ²⁺ from its isoelectronic Group 11 M₂ . While "classical" bonding indices are inconclusive in establishing regions of metastability in the bonding, it is shown that the one-electron Bohm quantum potential is promising in this regard.

Charge densities in actinide compounds: strategies for data reduction and model building

The data quality requirements for charge density studies on actinide compounds are extreme. Important steps in data collection and reduction required to obtain such data are summarized and evaluated. The steps involved in building an augmented Hansen-Coppens multipole model for an actinide pseudo-atom are provided. The number and choice of radial functions, in particular the definition of the core, valence and pseudo-valence terms are discussed. The conclusions in this paper are based on a re-examination and improvement of a previously reported study on [PPh4][UF6]. Topological analysis of the total electron density shows remarkable agreement between experiment and theory; however, there are significant differences in the Laplacian distribution close to the uranium atoms which may be due to the effective core potential employed for the theoretical calculations.

Valence-Shell Electron-Pair Repulsion Theory Revisited: An Explanation for Core Polarization

Valence-shell electron-pair repulsion (VSEPR) theory constitutes one of the pillars of theoretical predictive chemistry. It was proposed even before the advent of the concept of "spin", and it is still a very useful tool in chemistry. In this article we propose an extension of VSEPR theory to understand the core structure and predict core polarization in the main-group elements. We show from first principles (Electron Localization Function analysis) how the inner- and outer-core shells are organized. In particular, electrons in these regions are structured following the shape of the dual polyhedron of the valence shell (3^rd period) or the equivalent polyhedron (4^th and 5^th periods). We interpret these results in terms of "hard" and "soft" core character. All the studied systems follow this trend, providing a framework for predicting electron distribution in the core. We also show that lone pairs behave as "standard ligands" in terms of core polarization. The predictive character of the model was tested by proposing the core polarization in different systems not included in the original set (such as XeF₄ and [Fe(CN)₆ ]^3- ) and checking the hypothesis by means of a posteriori calculations. From the experimental point of view, the extension of VSEPR to the core region has consequences for current crystallography research. In particular, it explains the core polarization revealed by high resolution X-ray experiments.

Molecular QTAIM Topology Is Sensitive to Relativistic Corrections

The topology of the molecular electron density of benzene dithiol gold cluster complex Au₄ -S-C₆ H₄ -S'-Au'₄ changed when relativistic corrections were made and the structure was close to a minimum of the Born-Oppenheimer energy surface. Specifically, new bond paths between hydrogen atoms on the benzene ring and gold atoms appeared, indicating that there is a favorable interaction between these atoms at the relativistic level. This is consistent with the observation that gold becomes a better electron acceptor when relativistic corrections are applied. In addition to relativistic effects, here, we establish the sensitivity of molecular topology to basis sets and convergence thresholds for geometry optimization.

Linear MgCp*2 vs Bent CaCp*2: London Dispersion, Ligand-Induced Charge Localizations, and Pseudo-Pregostic C-H···Ca Interactions

In the family of metallocenes, MgCp*₂ (Cp* = pentamethylcyclopentadienyl) exhibits a regular linear sandwich structure, whereas CaCp*₂ is bent in both the gas phase and solid state. Bending is typically observed for metal ions which possess a lone pair. Here, we investigate which electronic differences cause the bending in complexes lacking lone pairs at the metal atoms. The bent gas-phase geometry of CaCp*₂ suggests that the bending must have an intramolecular origin. Geometry optimizations with and without dispersion effects/d-type polarization functions on MCp₂ and MCp*₂ gas-phase complexes (M = Ca, Mg) establish that attractive methyl···methyl London dispersion interactions play a decisive role in the bending in CaCp*₂. A sufficient polarizability of the metal to produce a shallow bending potential energy curve is a prerequisite but is not the reason for the bending. Concomitant ligand-induced charge concentrations and localizations at the metal atoms are studied in further detail, for which real-space bonding and orbital-based descriptors are used. Low-temperature crystal structures of MgCp*₂ and CaCp*₂ were determined which facilitated the identification and characterization of intermolecular pseudo-pregostic interactions, C-H···Ca, in the CaCp*₂ crystal structure.

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.

Accurate charge densities from powder X-ray diffraction - a new version of the Aarhus vacuum imaging-plate diffractometer

In recent years powder X-ray diffraction has proven to be a valuable alternative to single-crystal X-ray diffraction for determining electron-density distributions in high-symmetry inorganic materials, including subtle deformation in the core electron density. This was made possible by performing diffraction measurements in vacuum using high-energy X-rays at a synchrotron-radiation facility. Here we present a new version of our custom-built in-vacuum powder diffractometer with the sample-to-detector distance increased by a factor of four. In practice this is found to give a reduction in instrumental peak broadening by approximately a factor of three and a large improvement in signal-to-background ratio compared to the previous instrument. Structure factors of silicon at room temperature are extracted using a combined multipole-Rietveld procedure and compared with ab initio calculations and the results from the previous diffractometer. Despite some remaining issues regarding peak asymmetry, the new diffractometer yields structure factors of comparable accuracy to the previous diffractometer at low angles and improved accuracy at high angles. The high quality of the structure factors is further assessed by modelling of core electron deformation with results in good agreement with previous investigations.

Experimental charge-density studies: data reduction and model quality: the more the better?

In this review, recent developments concerning data and model quality in experimental charge-density investigations from a personal view-point are described. Data quality is not only achieved by the high resolution, high I/σ(I) values, low merging R values and high multiplicity. The quality of the innermost reflections especially is crucial for mapping the density distribution of the outermost valence electrons and can be monitored by (I/σ)^asymptotic. New detector technologies seem to be promising improvements. Empirical corrections to correct for low-energy contamination of mirror-focused X-ray data and for resolution- and temperature-dependent errors caused by factors such as thermal diffuse scattering are described. Shashlik-like residual density patterns can indicate the need for an anharmonic description of the thermal motion of individual atoms. The physical reliability of the derived model must be thoroughly analysed. The derived probability density functions for the mean-squared atomic vibrational displacements especially should have only small negative values. The treatment of H atoms has been improved by methods to estimate anisotropic thermal motion. For very high resolution data, the polarization of the core density cannot be neglected. Several tools to detect systematic errors are described. A validation tool is presented that easily detects when the refinement of additional parameters yields a real improvement in the model or simply overfits the given data. In all investigated structures, it is proved that the multipole parameters of atoms with a comparable chemical environment should be constrained to be identical. The use of restraints could be a promising alternative.

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.

Communication: Determination of relativistic effects from X-ray structure factors

In this communication, a procedure is presented which allows for the determination of the scalar-relativistic contraction of individual electronic shells of transition metal atoms from X-ray structure factor data. The procedure is verified and benchmarked employing theoretical and experimental F(hkl) data, revealing an overall good agreement between the experimentally determined results and the theoretical reference values. From the experimental data, the relativistic contraction of the n = 2 shell of a cerium atom is, for example, determined as 0.097 pm, compared to a theoretical reference value of 0.116 pm. It is further demonstrated that the reproducibility of the results is excellent when comparing different experimental data sets. Finally, the dependency of the according results on the data resolution of the structure factor data is investigated.

Experimental and theoretical charge-density analysis of 1,4-bis(5-hexyl-2-thienyl)butane-1,4-dione: applications of a virtual-atom model

The experimental and theoretical charge densities of 1,4-bis(5-hexyl-2-thienyl)butane-1,4-dione, a precursor in the synthesis of thiophene-based semiconductors and organic solar cells, are presented. A dummy bond charges spherical atom model is applied besides the multipolar atom model. The results show that the dummy bond charges model is accurate enough to calculate electrostatic-derived properties which are comparable with those obtained by the multipolar atom model. The refinement statistics and the residual electron density values are found to be intermediate between the independent atom and the multipolar formalisms.

Synchrotron powder diffraction of silicon: high-quality structure factors and electron density

Crystalline silicon is an ideal compound to test the current state of experimental structure factors and corresponding electron densities. High-quality structure factors have been measured on crystalline silicon with synchrotron powder X-ray diffraction. They are in excellent agreement with benchmark Pendellösung data having comparable accuracy and precision, but acquired in far less time and to a much higher resolution (sin θ/λ < 1.7 Å(-1)). The extended data range permits an experimental modelling of not only the valence electron density but also the core deformation in silicon, establishing an increase of the core density upon bond formation in crystalline silicon. Furthermore, a physically sound procedure for evaluating the standard deviation of powder-derived structure factors has been applied. Sampling statistics inherently account for contributions from photon counts as well as the limited number of diffracting particles, where especially the latter are particularly difficult to handle.

A charge density study of π-delocalization and intermolecular interactions

The location of bond critical points (red dots) and its associated bond path (black line) provide the evidence on the existence of the weak intermolecular interactions of the π–π interactions between triazole rings of atrz molecules in crystal with the close ring distance of 3.17 Å.