Role of Charge Density and Surface Area of Tailored Ionic Porous Organic Polymers for Adsorption and Antibacterial Actions

The development of high-performance adsorbents for environmental remediation is a current need, and ionic porous organic polymers (iPOPs), due to their high physicochemical stability, high surface area, added electrostatic interaction, and easy reusability, have already established themselves as a better adsorbent. However, research on the structural design of high-performance iPOP-based adsorbents is still nascent. This study explored the building blocks' role in optimizing the polymers' charge density and surface area to develop better polymeric adsorbents. Among the three synthesized polymers, iPOP-ZN1, owing to its high surface area and high charge density in its active sites, proved to be the best adsorbent for adsorbing inorganic and organic pollutants in an aqueous medium. The polymers were efficient enough to capture and store iodine vapor in the solid state. Further, this study tried to address using iodine-loaded polymers in antibacterial action. Iodine-loaded iPOPs show impressive antibacterial behavior against E. coli, B. subtilis, and H. pylori.

On the structure refinement of metal complexes against 3D electron diffraction data using multipolar scattering factors

This study examines various methods for modelling the electron density and, thus, the electrostatic potential of an organometallic complex for use in crystal structure refinement against 3D electron diffraction (ED) data. It focuses on modelling the scattering factors of iron(III), considering the electron density distribution specific for coordination with organic linkers. We refined the structural model of the metal-organic complex, iron(III) acetylacetonate (FeAcAc), using both the independent atom model (IAM) and the transferable aspherical atom model (TAAM). TAAM refinement initially employed multipolar parameters from the MATTS databank for acetylacetonate, while iron was modelled with a spherical and neutral approach (TAAM ligand). Later, custom-made TAAM scattering factors for Fe-O coordination were derived from DFT calculations [TAAM-ligand-Fe(III)]. Our findings show that, in this compound, the TAAM scattering factor corresponding to Fe3+ has a lower scattering amplitude than the Fe3+ charged scattering factor described by IAM. When using scattering factors corresponding to the oxidation state of iron, IAM inaccurately represents electrostatic potential maps and overestimates the scattering potential of the iron. In addition, TAAM significantly improved the fitting of the model to the data, shown by improved R1 values, goodness-of-fit (GooF) and reduced noise in the Fourier difference map (based on the residual distribution analysis). For 3D ED, R1 values improved from 19.36% (IAM) to 17.44% (TAAM-ligand) and 17.49% (TAAM-ligand-Fe3+), and for single-crystal X-ray diffraction (SCXRD) from 3.82 to 2.03% and 1.98%, respectively. For 3D ED, the most significant R1 reductions occurred in the low-resolution region (8.65-2.00 Å), dropping from 20.19% (IAM) to 14.67% and 14.89% for TAAM-ligand and TAAM-ligand-Fe(III), respectively, with less improvement in high-resolution ranges (2.00-0.85 Å). This indicates that the major enhancements are due to better scattering modelling in low-resolution zones. Furthermore, when using TAAM instead of IAM, there was a noticeable improvement in the shape of the thermal ellipsoids, which more closely resembled those of an SCXRD-refined model. This study demonstrates the applicability of more sophisticated scattering factors to improve the refinement of metal-organic complexes against 3D ED data, suggesting the need for more accurate modelling methods and highlighting the potential of TAAM in examining the charge distribution of large molecular structures using 3D ED.

An orbital strategy for regulating the Jahn-Teller effect

The Jahn-Teller effect (JTE) arising from lattice-electron coupling is a fascinating phenomenon that profoundly affects important physical properties in a number of transition-metal compounds. Controlling JT distortions and their corresponding electronic structures is highly desirable to tailor the functionalities of materials. Here, we propose a local coordinate strategy to regulate the JTE through quantifying occupancy in the [Formula: see text] and [Formula: see text] orbitals of Mn and scrutinizing the symmetries of the ligand oxygen atoms in MnO6 octahedra in LiMn2O4 and Li0.5Mn2O4. The effectiveness of such a strategy has been demonstrated by constructing P2-type NaLi x Mn1 - x O2 oxides with different Li/Mn ordering schemes. In addition, this strategy is also tenable for most 3d transition-metal compounds in spinel and perovskite frameworks, indicating the universality of local coordinate strategy and the tunability of the lattice-orbital coupling in transition-metal oxides. This work demonstrates a useful strategy to regulate JT distortion and provides useful guidelines for future design of functional materials with specific physical properties.

Dynamical refinement with multipolar electron scattering factors

Dynamical refinement is a well established method for refining crystal structures against 3D electron diffraction (ED) data and its benefits have been discussed in the literature [Palatinus, Petříček & Corrêa, (2015). Acta Cryst. A71, 235-244; Palatinus, Corrêa et al. (2015). Acta Cryst. B71, 740-751]. However, until now, dynamical refinements have only been conducted using the independent atom model (IAM). Recent research has shown that a more accurate description can be achieved by applying the transferable aspherical atom model (TAAM), but this has been limited only to kinematical refinements [Gruza et al. (2020). Acta Cryst. A76, 92-109; Jha et al. (2021). J. Appl. Cryst. 54, 1234-1243]. In this study, we combine dynamical refinement with TAAM for the crystal structure of 1-methyluracil, using data from precession ED. Our results show that this approach improves the residual Fourier electrostatic potential and refinement figures of merit. Furthermore, it leads to systematic changes in the atomic displacement parameters of all atoms and the positions of hydrogen atoms. We found that the refinement results are sensitive to the parameters used in the TAAM modelling process. Though our results show that TAAM offers superior performance compared with IAM in all cases, they also show that TAAM parameters obtained by periodic DFT calculations on the refined structure are superior to the TAAM parameters from the UBDB/MATTS database. It appears that multipolar parameters transferred from the database may not be sufficiently accurate to provide a satisfactory description of all details of the electrostatic potential probed by the 3D ED experiment.

Comparative study of conventional and synchrotron X-ray electron densities on molecular crystals

Five different electron density datasets obtained from conventional and synchrotron single crystal X-ray diffraction experiments are compared. The general aim of the study is to investigate the quality of data for electron density analysis from current state-of-the-art conventional sources, and to see how the data perform in comparison with high-quality synchrotron data. A molecular crystal of melamine was selected as the test compound due to its ability to form excellent single crystals, the light atom content, and an advantageous suitability factor of 3.6 for electron density modeling. These features make melamine an optimal system for conventional X-ray diffractometers since the inherent advantages of synchrotron sources such as short wavelength and high intensity are less critical in this case. Data were obtained at 100 K from new in-house diffractometers Rigaku Synergy-S (Mo and Ag source, HyPix100 detector) and Stoe Stadivari (Mo source, EIGER2 1M CdTe detector), and an older Oxford Diffraction Supernova (Mo source, Atlas CCD detector). The synchrotron data were obtained at 25 K from BL02B1 beamline at SPring-8 in Japan (λ = 0.2480 Å, Pilatus3 X 1M CdTe detector). The five datasets were compared on general quality parameters such as resolution, ⟨I/σ⟩, redundancy and R factors, as well as the more model specific fractal dimension plot and residual density maps. Comparison of the extracted electron densities reveals that all datasets can provide reliable multipole models, which overall convey similar chemical information. However, the new laboratory X-ray diffractometers with advanced pixel detector technology clearly measure data with significantly less noise and much higher reliability giving densities of higher quality, compared to the older instrument. The synchrotron data have higher resolution and lower measurement temperature, and they allow for finer details to be modeled (e.g. hydrogen κ parameters).

Revisiting the Nature of Chemical Bonding in Chalcogenides to Explain and Design their Properties

Quantum chemical bonding descriptors have recently been utilized to design materials with tailored properties. Their usage to facilitate a quantitative description of bonding in chalcogenides as well as the transition between different bonding mechanisms is reviewed. More importantly, these descriptors can also be employed as property predictors for several important material characteristics, including optical and transport properties. Hence, these quantum chemical bonding descriptors can be utilized to tailor material properties of chalcogenides relevant for thermoelectrics, photovoltaics, and phase-change memories. Relating material properties to bonding mechanisms also shows that there is a class of materials, which are characterized by unconventional properties such as a pronounced anharmonicity, a large chemical bond polarizability, and strong optical absorption. This unusual property portfolio is attributed to a novel bonding mechanism, fundamentally different from ionic, metallic, and covalent bonding, which is called "metavalent." In the concluding section, a number of promising research directions are sketched, which explore the nature of the property changes upon changing bonding mechanism and extend the concept of quantum chemical property predictors to more complex compounds.

Electron density and thermal motion of diamond at elevated temperatures

The electron density and thermal motion of diamond are determined at nine temperatures between 100 K and 1000 K via synchrotron powder X-ray diffraction (PXRD) data collected on a high-accuracy detector system. Decoupling of the thermal motion from the thermally smeared electron density is performed via an iterative Wilson-Hansen-Coppens-Rietveld procedure using theoretical static structure factors from density functional theory (DFT) calculations. The thermal motion is found to be harmonic and isotropic in the explored temperature range, and excellent agreement is observed between experimental atomic displacement parameters (ADPs) and those obtained via theoretical harmonic phonon calculations (HPC), even at 1000 K. The Debye temperature of diamond is determined experimentally to be Θ_D = 1883 (35) K. A topological analysis of the electron density explores the temperature dependency of the electron density at the bond critical point. The properties are found to be constant throughout the temperature range. The robustness of the electron density confirms the validity of the crystallographic convolution approximation for diamond in the explored temperature range.

Real-space measurement of orbital electron populations for Li1-xCoO2

The operation of lithium-ion batteries involves electron removal from and filling into the redox orbitals of cathode materials, experimentally probing the orbital electron population thus is highly desirable to resolve the redox processes and charge compensation mechanism. Here, we combine quantitative convergent-beam electron diffraction with high-energy synchrotron powder X-ray diffraction to quantify the orbital populations of Co and O in the archetypal cathode material LiCoO₂. The results indicate that removing Li ions from LiCoO₂ decreases Co t_2g orbital population, and the intensified covalency of Co-O bond upon delithiation enables charge transfer from O 2p orbital to Co e_g orbital, leading to increased Co e_g orbital population and oxygen oxidation. Theoretical calculations verify these experimental findings, which not only provide an intuitive picture of the redox reaction process in real space, but also offer a guidance for designing high-capacity electrodes by mediating the covalency of the TM-O interactions.

X‐ray Electron Density Study of the Chemical Bonding Origin of Glass Formation in Metal–Organic Frameworks**

Glass‐forming metal–organic frameworks (MOFs) have novel applications, but the origin of their peculiar melting behavior is unclear. Here, we report synchrotron X‐ray diffraction electron densities of two zeolitic imidazolate frameworks (ZIFs), the glass‐forming Zn‐ZIF‐zni and the isostructural thermally decomposing Co‐ZIF‐zni. Electron density analysis shows that the Zn−N bonds are more ionic than the Co−N bonds, which have distinct covalent features. Variable‐temperature Raman spectra reveal the onset of significant imidazolate bond weakening in Co‐ZIF‐zni above 673 K. Melting can be controlled by tuning the metal–ligand and imidazole bonding strength as shown from thermal analysis of nine solid‐solution Co_xZn_1−x‐ZIF‐zni (x=0.3 to 0.003) MOFs, and a mere 4 % Co‐doping into Zn‐ZIF‐zni results in thermal decomposition instead of melting. The present findings demonstrate the key role of the metal–ligand bonds and imidazolate bonds in controlling the delicate balance between melting and decomposition processes in this class of ZIF compounds.

Large scale dataset of real space electronic charge density of cubic inorganic materials from density functional theory (DFT) calculations

Driven by the big data science, material informatics has attracted enormous research interests recently along with many recognized achievements. To acquire knowledge of materials by previous experience, both feature descriptors and databases are essential for training machine learning (ML) models with high accuracy. In this regard, the electronic charge density ρ(r), which in principle determines the properties of materials at their ground state, can be considered as one of the most appropriate descriptors. However, the systematic electronic charge density ρ(r) database of inorganic materials is still in its infancy due to the difficulties in collecting raw data in experiment and the expensive first-principles based computational cost in theory. Herein, a real space electronic charge density ρ(r) database of 17,418 cubic inorganic materials is constructed by performing high-throughput density functional theory calculations. The displayed ρ(r) patterns show good agreements with those reported in previous studies, which validates our computations. Further statistical analysis reveals that it possesses abundant and diverse data, which could accelerate ρ(r) related machine learning studies. Moreover, the electronic charge density database will also assists chemical bonding identifications and promotes new crystal discovery in experiments.

Estimating electron density at the bond critical point through atomic Hirshfeld surfaces

Interatomic surfaces often carry information related to the electron distribution in a molecule or crystal, not only being a visual aid but also enabling quantitative analyses. Under certain conditions, atomic Hirshfeld surfaces present a high resemblance to the interatomic surfaces obtained through the quantum theory of atoms in molecules (QTAIM), with the advantage of being easily calculated, even for crystal structures determined at low resolutions (i.e. when a charge-density refinement is not performed). Here an empirical relation between the curvedness property of the Hirshfeld surfaces and the electron density at the bond critical point for certain types of covalent and coordination interactions involving carbon atoms has been obtained. The exponential function was tested to estimate the electron density in different crystalline systems, and the highest deviation from reference values obtained through QTAIM was just 16%. Additionally, fine details of this fit may be salient to the difference in electronegativity of the atoms involved in the bond.

The explicit role of electron exchange in the hydrogen bonded molecular complexes

We applied a set of advanced bonding descriptors to establish the hidden electron density features and binding energy characteristics of intermolecular DH∙∙∙A hydrogen bonds (OH∙∙∙O, NH∙∙∙O and SH∙∙∙O) in 150 isolated and solvated molecular complexes. The exchange-correlation and Pauli potentials as well as corresponding local one-electron forces allowed us to explicitly ascertain how electron exchange defines the bonding picture in the proximity of the H-bond critical point. The electron density features of DH∙∙∙A interaction are governed by alterations in the electron localization in the H-bond region displaying itself in the exchange hole. At that, they do not depend on the variations in the exchange hole mobility. The electrostatic interaction mainly defines the energy of H-bonds of different types, whereas the strengthening/weakening of H-bonds in complexes with varying substituents depends on the barrier height of the exchange potential near the bond critical point. Energy variations between H-bonds in isolated and solvated systems are also caused the electron exchange peculiarities as follows from the corresponding potential and the interacting quantum atom analyses complemented by electron delocalization index calculations. Our approach is based on the bonding descriptors associated with the characteristics of the observable electron density and can be recommended for in-depth studies of non-covalent bonding.

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.

Topological analysis of chemical bonding in the layered FePSe3 upon pressure-induced phase transitions

Two pressure-induced phase transitions have been theoretically studied in the layered iron phosphorus triselenide (FePSe₃ ). Topological analysis of chemical bonding in FePSe₃ has been performed based on the results of first-principles calculations within the periodic linear combination of atomic orbitals (LCAO) method with hybrid Hartree-Fock-DFT B3LYP functional. The first transition at about 6 GPa is accompanied by the symmetry change from R 3 ¯ to C2/m, whereas the semiconductor-to-metal transition (SMT) occurs at about 13 GPa leading to the symmetry change from C2/m to P 3 ¯ 1 m . We found that the collapse of the band gap at about 13 GPa occurs due to changes in the electronic structure of FePSe₃ induced by relative displacements of phosphorus or selenium atoms along the c-axis direction under pressure. The results of the topological analysis of the electron density and its Laplacian demonstrate that the pressure changes not only the interatomic distances but also the bond nature between the intralayer and interlayer phosphorus atoms. The interlayer P-P interactions are absent in two non-metallic FePSe₃ phases while after SMT the intralayer P-P interactions weaken and the interlayer P-P interactions appear.

Quantification of the Magnetic Anisotropy of a Single-Molecule Magnet from the Experimental Electron Density

Reported here is an entirely new application of experimental electron density (EED) in the study of magnetic anisotropy of single-molecule magnets (SMMs). Among those SMMs based on one single transition metal, tetrahedral Co^II-complexes are prominent, and their large zero-field splitting arises exclusively from coupling between the d x 2 - y 2 and d_xy orbitals. Using very low temperature single-crystal synchrotron X-ray diffraction data, an accurate electron density (ED) was obtained for a prototypical SMM, and the experimental d-orbital populations were used to quantify the d_xy-d x 2 - y 2 coupling, which simultaneously provides the composition of the ground-state Kramers doublet wave function. Based on this experimentally determined wave function, an energy barrier for magnetic relaxation in the range 193-268 cm^-1 was calculated, and is in full accordance with the previously published value of 230 cm^-1 obtained from near-infrared spectroscopy. These results provide the first clear and direct link between ED and molecular magnetic properties.

Deep Learning the Hohenberg-Kohn Maps of Density Functional Theory

A striking consequence of the Hohenberg-Kohn theorem of density functional theory is the existence of a bijection between the local density and the ground-state many-body wave function. Here we study the problem of constructing approximations to the Hohenberg-Kohn map using a statistical learning approach. Using supervised deep learning with synthetic data, we show that this map can be accurately constructed for a chain of one-dimensional interacting spinless fermions in different phases of this model including the charge ordered Mott insulator and metallic phases and the critical point separating them. However, we also find that the learning is less effective across quantum phase transitions, suggesting an intrinsic difficulty in efficiently learning nonsmooth functional relations. We further study the problem of directly reconstructing complex observables from simple local density measurements, proposing a scheme amenable to statistical learning from experimental data.

Chemical Bonding in Colossal Thermopower FeSb2

FeSb₂ exhibits a colossal Seebeck coefficient ( S ) and a record-breaking high thermoelectric power factor. It also has an atypical shift from diamagnetism to paramagnetism with increasing temperature, and the fine details of its electron correlation effects have been widely discussed. The extraordinary physical properties must be rooted in the nature of the chemical bonding, and indeed, the chemical bonding in this archetypical marcasite structure has been heavily debated on a theoretical basis since the 1960s. The two prevalent models for describing the bonding interactions in FeSb₂ are based on either ligand-field stabilization of Fe or a network structure of Sb hosting Fe ions. However, neither model can account for the observed properties of FeSb₂ . Herein, an experimental electron density study is reported, which is based on analysis of synchrotron X-ray diffraction data measured at 15 K on a minute single crystal to limit systematic errors. The analysis is supplemented with density functional theory calculations in the experimental geometry. The experimental data are at variance with both the additional single-electron Sb-Sb bond implied by the covalent model, and the large formal charge and expected d-orbital splitting advocated by the ionic model. The structure is best described as an extended covalent network in agreement with expectations based on electronegativity differences.

Accurate high-resolution single-crystal diffraction data from a Pilatus3 X CdTe detector

Hybrid photon-counting detectors are widely established at third-generation synchrotron facilities and the specifications of the Pilatus3 X CdTe were quickly recognized as highly promising in charge-density investigations. This is mainly attributable to the detection efficiency in the high-energy X-ray regime, in combination with a dynamic range and noise level that should overcome the perpetual problem of detecting strong and weak data simultaneously. These benefits, however, come at the expense of a persistent problem for high diffracted beam flux, which is particularly problematic in single-crystal diffraction of materials with strong scattering power and sharp diffraction peaks. Here, an in-depth examination of data collected on an inorganic material, FeSb2, and an organic semiconductor, rubrene, revealed systematic differences in strong intensities for different incoming beam fluxes, and the implemented detector intensity corrections were found to be inadequate. Only significant beam attenuation for the collection of strong reflections was able to circumvent this systematic error. All data were collected on a bending-magnet beamline at a third-generation synchrotron radiation facility, so undulator and wiggler beamlines and fourth-generation synchrotrons will be even more prone to this error. On the other hand, the low background now allows for an accurate measurement of very weak intensities, and it is shown that it is possible to extract structure factors of exceptional quality using standard crystallographic software for data processing (SAINT-Plus, SADABS and SORTAV), although special attention has to be paid to the estimation of the background. This study resulted in electron-density models of substantially higher accuracy and precision compared with a previous investigation, thus for the first time fulfilling the promise of photon-counting detectors for very accurate structure factor measurements.

Expression and interactions of stereochemically active lone pairs and their relation to structural distortions and thermal conductivity

In chemistry, stereochemically active lone pairs are typically described as an important non-bonding effect, and recent interest has centred on understanding the derived effect of lone pair expression on physical properties such as thermal conductivity. To manipulate such properties, it is essential to understand the conditions that lead to lone pair expression and provide a quantitative chemical description of their identity to allow comparison between systems. Here, density functional theory calculations are used first to establish the presence of stereochemically active lone pairs on antimony in the archetypical chalcogenide MnSb2O4. The lone pairs are formed through a similar mechanism to those in binary post-transition metal compounds in an oxidation state of two less than their main group number [e.g. Pb(II) and Sb(III)], where the degree of orbital interaction (covalency) determines the expression of the lone pair. In MnSb2O4 the Sb lone pairs interact through a void space in the crystal structure, and their their mutual repulsion is minimized by introducing a deflection angle. This angle increases significantly with decreasing Sb-Sb distance introduced by simulating high pressure, thus showing the highly destabilizing nature of the lone pair interactions. Analysis of the chemical bonding in MnSb2O4 shows that it is dominated by polar covalent interactions with significant contributions both from charge accumulation in the bonding regions and from charge transfer. A database search of related ternary chalcogenide structures shows that, for structures with a lone pair (SbX 3 units), the degree of lone pair expression is largely determined by whether the antimony-chalcogen units are connected or not, suggesting a cooperative effect. Isolated SbX 3 units have larger X-Sb-X bond angles and therefore weaker lone pair expression than connected units. Since increased lone pair expression is equivalent to an increased orbital interaction (covalent bonding), which typically leads to increased heat conduction, this can explain the previously established correlation between larger bond angles and lower thermal conductivity. Thus, it appears that for these chalcogenides, lone pair expression and thermal conductivity may be related through the degree of covalency of the system.