Relativistic X-ray Elastic Scattering Factors for Neutral AtomsZ= 1–54 from Multiconfiguration Dirac–Fock Wavefunctions in the 0–12Å−1sinθ /λ Range, and Six-Gaussian Analytical Expressions in the 0–6Å−1Range

Real-Space Observation of Ligand Hole State in Cubic Perovskite SrFeO3

An anomalously high valence state sometimes shows up in transition-metal oxide compounds. In such systems, holes tend to occupy mainly the ligand p orbitals, giving rise to interesting physical properties such as superconductivity in cuprates and rich magnetic phases in ferrates. However, no one has ever observed the distribution of ligand holes in real space. Here, a successful observation of the spatial distribution of valence electrons in cubic perovskite SrFeO3 by high-energy X-ray diffraction experiments and precise electron density analysis using a core differential Fourier synthesis method is reported. A real-space picture of ligand holes formed by the orbital hybridization of Fe 3d and O 2p is revealed. The anomalous valence state in Fe is attributed to the considerable contribution of the ligand hole, which is related to the metallic nature and the absence of Jahn-Teller distortions in this system.

Applicability of transferable multipole pseudo-atoms for restoring inner-crystal electronic force density fields. Chemical bonding and binding features in the crystal and dimer of 1,3-bis(2-hydroxyethyl)-6-methyluracil

We considered it timely to test the applicability of transferable multipole pseudo-atoms for restoring inner-crystal electronic force density fields. The procedure was carried out on the crystal of 1,3-bis(2-hydroxyethyl)-6-methyluracil, and some derived properties of the scalar potential and vector force fields were compared with those obtained from the experimental multipole model and from the aspherical pseudo-atom model with parameters fitted to the calculated structure factors. The procedure was shown to accurately replicate the general vector-field behavior, the peculiarities of the quantum potentials and the characteristics of the force-field pseudoatoms, such as charge, shape and volume, as well as to reproduce the relative arrangement of atomic and pseudoatomic zero-flux surfaces along internuclear regions. It was found that, in addition to the quantum-topological atoms, the force-field pseudoatoms are spatially reproduced within a single structural fragment and similar environment. In addition, the classical and nonclassical hydrogen bonds in the uracil derivative crystal, as well as the H...O, N...O and N...C interactions in the free π-stacked dimer of the uracil derivative molecules, were studied using the potential and force fields within the concepts of interatomic charge transfer and electron lone pair donation-acceptance. Remarkably, the nitrogen atoms in the N...O and N...C interactions behave rather like a Lewis base and an electron contributor. At the same time, the hydrogen atom in the H...O interaction, being a Lewis acid, also participates in the interatomic electron transfer by acting as a contributor. Thus, it has been argued that, when describing polar interatomic interactions within orbital-free considerations, it makes more physical sense to identify electronegative (electron occupier) and electropositive (electron contributor) atoms or subatomic fragments rather than nucleophilic and electrophilic sites.

Revisited relativistic Dirac-Hartree-Fock X-ray scattering factors. I. Neutral atoms with Z = 2-118

In this first of a series of publications, the X-ray scattering factors for neutral atoms are revisited. Using the recently developed DBSR_HF program [Zatsarinny & Froese Fischer (2016). Comput. Phys. Comm. 202, 287-303] the fully relativistic Dirac-Hartree-Fock ground-state wavefunctions for all atoms with Z = 2-118 (He-Og) have been calculated using the extended average level scheme and including both the Breit interaction correction to the electronic motion due to magnetic and retardation effects, and the Fermi distribution function for the description of the nuclear charge density. The comparison of our wavefunctions with those obtained in several previous studies in terms of the total and orbital (spinor) electronic energies, and a number of local and integrated total and orbital properties, confirmed the quality of the generated wavefunctions. The employed dense radial grid combined with the DBSR_HF's B-spline representation of the relativistic one-electron orbitals allowed for a precise integration of the X-ray scattering factors using a newly developed Fortran program SF. Following the established procedure [Maslen et al. (2006). International Tables for Crystallography, Vol. C, Section 6.1.1, pp. 554-589], the resulting X-ray scattering factors have been interpolated in the 0 ≤ sin θ/λ ≤ 2 Å^-1 and 2 ≤ sin θ/λ ≤ 6 Å^-1 ranges using the recommended analytical functions with both the four- (which is a current convention) and five-term expansions. An exhaustive comparison of the newly generated X-ray scattering factors with the International Union of Crystallography recommended values and those from a number of previous studies showed an overall good agreement and allowed identification of a number of typos and inconsistencies in the recommended quantities. A detailed analysis of the results suggests that the newly derived values may represent an excellent compromise among all the previous studies. The determined conventional interpolating functions for the two sin θ/λ intervals show, on average, the same accuracy as the recommended parametrizations. However, an extension of each expansion by only a single term provides a significant improvement in the accuracy of the interpolated values for an overwhelming majority of the atoms. As such, an updated set of the fully relativistic X-ray scattering factors and the interpolating functions for neutral atoms with Z = 2-118 can be easily incorporated into the existing X-ray diffraction software with only minor modifications. The outcomes of the undertaken research should be of interest to members of the crystallographic community who push the boundaries of the accuracy and precision of X-ray diffraction studies.

Anomalous temperature dependence of the experimental x-ray structure factor of supercooled water

The structural changes of water upon deep supercooling were studied through wide-angle x-ray scattering at SwissFEL. The experimental setup had a momentum transfer range of 4.5 Å^-1, which covered the principal doublet of the x-ray structure factor of water. The oxygen-oxygen structure factor was obtained for temperatures down to 228.5 ± 0.6 K. Similar to previous studies, the second diffraction peak increased strongly in amplitude as the structural change accelerated toward a local tetrahedral structure upon deep supercooling. We also observed an anomalous trend for the second peak position of the oxygen-oxygen structure factor (q₂). We found that q₂ exhibits an unprecedented positive partial derivative with respect to temperature for temperatures below 236 K. Based on Fourier inversion of our experimental data combined with reference data, we propose that the anomalous q₂ shift originates from that a repeat spacing in the tetrahedral network, associated with all peaks in the oxygen-oxygen pair-correlation function, gives rise to a less dense local ordering that resembles that of low-density amorphous ice. The findings are consistent with that liquid water consists of a pentamer-based hydrogen-bonded network with low density upon deep supercooling.

Coexistence and Coupling of Multiple Charge Orderings and Spin States in Hexagonal Ferrite

The coupling between charge and spin orderings in strongly correlated systems plays a crucial role in fundamental physics and device applications. As a candidate of multiferroic materials, LuFe₂O₄ with a nominal Fe^2.5+ valence state has the potential for strong charge-spin interactions; however, these interactions have not been fully understood until now. Here, combining complementary characterization methods with theoretical calculations, two types of charge orderings with distinct magnetic properties are revealed. The ground states of LuFe₂O₄ are decided by the parallel/antiparallel coupling of both charge and spin orderings in the adjacent FeO double layers. Whereas the ferroelectric charge ordering remains ferrimagnetic below 230 K, the antiferroelectric ordering undergoes antiferromagnetic-ferrimagnetic-paramagnetic transitions from 2 K to room temperature. This study demonstrates the unique aspects of strong spin-charge coupling within LuFe₂O₄. Our results shed light on the coexistence and competing nature of orderings in quantum materials.

Orbital-Free Quantum Crystallographic View on Noncovalent Bonding: Insights into Hydrogen Bonds, π⋅⋅⋅π and Reverse Electron Lone Pair⋅⋅⋅π Interactions

A detailed analysis of a complete set of the local potentials that appear in the Euler equation for electron density is carried out for noncovalent interactions in the crystal of a uracil derivative using experimental X-ray charge density. The interplay between the quantum theory of atoms in molecules and crystals and the local potentials and corresponding inner-crystal electronic forces of electrostatic and kinetic origin is explored. Partitioning of crystal space into atomic basins and atomic-like potential basins led us to the definite description of interatomic interaction and charge transfer. Novel physically grounded bonding descriptors derived within the orbital-free quantum crystallography provided the detailed examination of π-stacking and intricate C=O⋅⋅⋅π interactions and nonclassical hydrogen bonds present in the crystal. The donor-acceptor character of these interactions is revealed by analysis of Pauli and von Weizsäcker potentials together with well-known functions, e. g., deformation electron density and electron localization function. In this way, our analysis throws light on aspects of these closed-shell interactions hitherto hidden from the description.

Direct Observation of Molecular Orbitals Using Synchrotron X-ray Diffraction

The physical properties of molecular crystals are governed by the frontier orbitals of molecules. A molecular orbital, which is formed by superposing the atomic orbitals of constituent elements, has complicated degrees of freedom in the crystal because of the influence of electron correlation and crystal field. Therefore, in general, it is difficult to experimentally observe the whole picture of a frontier orbital. Here, we introduce a new method called “core differential Fourier synthesis” (CDFS) using synchrotron X-ray diffraction to observe the valence electron density in materials. By observing the valence electrons occupied in molecular orbitals, the orbital state can be directly determined in a real space. In this study, we applied the CDFS method to molecular materials such as diamond, C60 fullerene, (MV)I2, and (TMTTF)2X. Our results not only demonstrate the typical orbital states in some materials, but also provide a new method for studying intramolecular degrees of freedom.

Experimental validation of bifurcated hydrogen bond of 2,5-lutidinium bromanilate and its charge density distribution

Abstract

2,5-Lutidinium bromanilate is a molecular complex that consists of bromanilic acid and 2,5-lutidine in which hydrogen-bonding interactions occur between them, producing a charge-assisted bifurcated N–H…O hydrogen bond. Bond characteristics are determined from the experimental charge density distribution of the molecular complex using the Hansen–Coppens model. The electron density, topological properties, electrostatic potential and atomic charges of the molecule have been investigated to better understand the atomic, molecular and electronic properties in a detailed manner. The electronic nature of the significantly important charge-assisted bifurcated hydrogen bond has been analyzed with the help of topological properties.

Graphic abstract

An ultrafast symmetry switch in a Weyl semimetal

Topological quantum materials exhibit fascinating properties^1-3, with important applications for dissipationless electronics and fault-tolerant quantum computers^4,5. Manipulating the topological invariants in these materials would allow the development of topological switching applications analogous to switching of transistors⁶. Lattice strain provides the most natural means of tuning these topological invariants because it directly modifies the electron-ion interactions and potentially alters the underlying crystalline symmetry on which the topological properties depend^7-9. However, conventional means of applying strain through heteroepitaxial lattice mismatch¹⁰ and dislocations¹¹ are not extendable to controllable time-varying protocols, which are required in transistors. Integration into a functional device requires the ability to go beyond the robust, topologically protected properties of materials and to manipulate the topology at high speeds. Here we use crystallographic measurements by relativistic electron diffraction to demonstrate that terahertz light pulses can be used to induce terahertz-frequency interlayer shear strain with large strain amplitude in the Weyl semimetal WTe₂, leading to a topologically distinct metastable phase. Separate nonlinear optical measurements indicate that this transition is associated with a symmetry change to a centrosymmetric, topologically trivial phase. We further show that such shear strain provides an ultrafast, energy-efficient way of inducing robust, well separated Weyl points or of annihilating all Weyl points of opposite chirality. This work demonstrates possibilities for ultrafast manipulation of the topological properties of solids and for the development of a topological switch operating at terahertz frequencies.

Aspherical scattering factors for SHELXL - model, implementation and application

A new aspherical scattering factor formalism has been implemented in the crystallographic least-squares refinement program SHELXL. The formalism relies on Gaussian functions and can optionally complement the independent atom model to take into account the deformation of electron-density distribution due to chemical bonding and lone pairs. Asphericity contributions were derived from the electron density obtained from quantum-chemical density functional theory computations of suitable model compounds that contain particular chemical environments, as defined by the invariom formalism. Thanks to a new algorithm, invariom assignment for refinement in SHELXL is automated. A suitable parameterization for each chemical environment within the new model was achieved by metaheuristics. Figures of merit, precision and accuracy of crystallographic least-squares refinements improve significantly upon using the new model.

X-ray Scattering and O-O Pair-Distribution Functions of Amorphous Ices

The structure factor and oxygen–oxygen pair-distribution functions of amorphous ices at liquid nitrogen temperature (T = 77 K) have been derived from wide-angle X-ray scattering (WAXS) up to interatomic distances of r = 23 Å, where local structure differences between the amorphous ices can be seen for the entire range. The distances to the first coordination shell for low-, high-, and very-high-density amorphous ice (LDA, HDA, VHDA) were determined to be 2.75, 2.78, and 2.80 Å, respectively, with high accuracy due to measurements up to a large momentum transfer of 23 Å^–1. Similarities in pair-distribution functions between LDA and supercooled water at 254.1 K, HDA and liquid water at 365.9 K, and VHDA and high-pressure liquid water were found up to around 8 Å, but beyond that at longer distances, the similarities were lost. In addition, the structure of the high-density amorphous ices was compared to high-pressure crystalline ices IV, IX , and XII, and conclusions were drawn about the local ordering.

Ionic scattering factors of atoms that compose biological molecules

Ionic scattering factors of atoms that compose biological molecules have been computed by the multi-configuration Dirac-Fock method. These ions are chemically unstable and their scattering factors had not been reported except for O-. Yet these factors are required for the estimation of partial charges in protein molecules and nucleic acids. The electron scattering factors of these ions are particularly important as the electron scattering curves vary considerably between neutral and charged atoms in the spatial-resolution range explored in structural biology. The calculated X-ray and electron scattering factors have then been parameterized for the major scattering curve models used in X-ray and electron protein crystallography and single-particle cryo-EM. The X-ray and electron scattering factors and the fitting parameters are presented for future reference.

Exploring the rare S—H...S hydrogen bond using charge density analysis in isomers of mercaptobenzoic acid

Experimental and theoretical charge density analyses on isomers of mercaptobenzoic acid have been carried out to quantify the hydrogen bonding of the hitherto less explored thiols, to assess the strength of the interactions using the topological features of the electron density. The electron density study offers interesting insights into the nature of the S—H...S interaction. The interaction energy is comparable with that of a weak hydrogen bond. The strength and directionality of the S—H...S hydrogen bond is demonstrated to be mainly due to the conformation locking potential of the intramolecular S...O chalcogen bond in 2-mercaptobenzoic acid and is stronger than in 3-mercaptobenzoic acid, which lacks the intramolecular S...O bond. Thepara-substituted mercaptobenzoic acid depicts a type I S...S interaction.

Testing the tools for revealing and characterizing the iodine-iodine halogen bond in crystals

To understand what tools are really suitable to identify and classify the iodine-iodine non-covalent interactions in solid organic polyiodides, we have examined the anisotropy of the electron density within the iodine atomic basin along and across the iodine-iodine halogen bond using the Laplacian of electron density, one-electron potential and electron localization function produced by Kohn-Sham calculations with periodic boundary conditions. The Laplacian of electron density exhibits the smallest anisotropy and yields a vague picture of the outermost electronic shells. The one-electron potential does not show such a deficiency and reveals that the valence electron shell for the halogen-bond acceptor iodine is always wider than that for the halogen-bond donor iodine along its σ-hole direction. We have concluded that the one-electron potential is the most suitable for classification of the iodine-iodine bonds and interactions in complicated cases, while the electron localization function allows to distinguish the diiodine molecule bonded with the monoiodide anion from the typical triiodide anion.

Extracting water and ion distributions from solution x-ray scattering experiments

Small-angle X-ray scattering measurements can provide valuable information about the solvent environment around biomolecules, but it can be difficult to extract solvent-specific information from observed intensity profiles. Intensities are proportional to the square of scattering amplitudes, which are complex quantities. Amplitudes in the forward direction are real, and the contribution from a solute of known structure (and from the waters it excludes) can be estimated from theory; hence, the amplitude arising from the solvent environment can be computed by difference. We have found that this "square root subtraction scheme" can be extended to non-zero q values, out to 0.1 Å(-1) for the systems considered here, since the phases arising from the solute and from the water environment are nearly identical in this angle range. This allows us to extract aspects of the water and ion distributions (beyond their total numbers), by combining experimental data for the complete system with calculations for the solutes. We use this approach to test molecular dynamics and integral-equation (3D-RISM (three-dimensional reference interaction site model)) models for solvent structure around myoglobin, lysozyme, and a 25 base-pair duplex DNA. Comparisons can be made both in Fourier space and in terms of the distribution of interatomic distances in real space. Generally, computed solvent distributions arising from the MD simulations fit experimental data better than those from 3D-RISM, even though the total small-angle X-ray scattering patterns are very similar; this illustrates the potential power of this sort of analysis to guide the development of computational models.

Spatial distribution of electrons near the Fermi level in the metallic LaB6 through accurate X-ray charge density study

Charge densities of iso-structural metal hexaborides, a transparent metal LaB₆ and a semiconductor BaB₆, have been determined using the d > 0.22 Å ultra-high resolution synchrotron radiation X-ray diffraction data by a multipole refinement and a maximum entropy method (MEM). The quality of the experimental charge densities was evaluated by comparison with theoretical charge densities. The strong inter-octahedral and relatively weak intra-octahedral boron-boron bonds were observed in the charge densities. A difference of valence charge densities between LaB₆ and BaB₆ was calculated to reveal a small difference between isostructural metal and semiconductor. The weak electron lobes distributed around the inter B₆ octahedral bond were observed in the difference density. We found the electron lobes are the conductive π-electrons in LaB₆ from the comparison with the theoretical valence charge density. We successfully observed a spatial distribution of electrons near the Fermi level from the X-ray charge density study of the series of iso-structural solids.

Molten uranium dioxide structure and dynamics

Uranium dioxide (UO2) is the major nuclear fuel component of fission power reactors. A key concern during severe accidents is the melting and leakage of radioactive UO2 as it corrodes through its zirconium cladding and steel containment. Yet, the very high temperatures (>3140 kelvin) and chemical reactivity of molten UO2 have prevented structural studies. In this work, we combine laser heating, sample levitation, and synchrotron x-rays to obtain pair distribution function measurements of hot solid and molten UO2. The hot solid shows a substantial increase in oxygen disorder around the lambda transition (2670 K) but negligible U-O coordination change. On melting, the average U-O coordination drops from 8 to 6.7 ± 0.5. Molecular dynamics models refined to this structure predict higher U-U mobility than 8-coordinated melts.

Experimental and database-transferred electron-density analysis and evaluation of electrostatic forces in coumarin-102 dye

The electron-density distribution of a new crystal form of coumarin-102, a laser dye, has been investigated using the Hansen-Coppens multipolar atom model. The charge density was refined versus high-resolution X-ray diffraction data collected at 100 K and was also constructed by transferring the charge density from the Experimental Library of Multipolar Atom Model (ELMAM2). The topology of the refined charge density has been analysed within the Bader `Atoms In Molecules' theory framework. Deformation electron-density peak heights and topological features indicate that the chromen-2-one ring system has a delocalized π-electron cloud in resonance with the N (amino) atom. The molecular electrostatic potential was estimated from both experimental and transferred multipolar models; it reveals an asymmetric character of the charge distribution across the molecule. This polarization effect is due to a substantial charge delocalization within the molecule. The molecular dipole moments derived from the experimental and transferred multipolar models are also compared with the liquid and gas-phase dipole moments. The substantial molecular dipole moment enhancements observed in the crystal environment originate from the crystal field and from intermolecular charge transfer induced and controlled by C-H···O and C-H···N intermolecular hydrogen bonds. The atomic forces were integrated over the atomic basins and compared for the two electron-density models.