Fragment Orbitals Extracted from First-Principles Plane-Wave Calculations.

Decomposing extended structures into smaller, molecular, even functional groups or simple fragments has a long tradition in chemistry because it allows for understanding certain electronic peculiarities in truly chemical terms. By doing so, invaluable property information is chemically accessible, for example, needed to rationalize catalytic or magnetic or optical nature. In order to also follow that train of thought for periodic materials, we have developed a tool which in a straightforward manner derives fragment molecular orbitals from plane-wave electronic-structure data of whatever kind of solid-state material. We here report on the mathematical apparatus of the method dubbed linear combination of fragment orbitals (LCFO) used for that purpose, implemented within the LOBSTER code. The method is illustrated from various sorts of molecular entities contained in such crystalline materials, together with an assessment of both accuracy and robustness of the new tool.

Na2C3O2: The First Crystalline Compound with the Elusive -C≡C-COO- Anion.

By reaction of sodium electride or NaC2H with the anhydrous sodium salt of propiolic acid, Na(OOC-C≡C-H), in liquid ammonia crystalline powders of Na2C3O2 were obtained. The structure analysis based on synchrotron powder diffraction data revealed that Na2C3O2 crystallizes in a monoclinic unit cell (I2/a, Z=4) exhibiting the elusive Y-shaped -C≡C-COO- anion, which is unprecedented in a crystalline compound up to now. IR/Raman and solid-state NMR spectroscopic investigations with assignments supported by DFT-based ab initio calculations confirm this finding. Reaction with sodium electride led to a higher crystallinity of the product, but additionally a by-product apparently due to decomposition and polymerization of Na2C3O2 was formed. No such by-product was observed in the reaction with NaC2H, which turned out to be a milder metalation route. However, the product of the latter reaction is less crystalline.

Switchable Chemical-Bond Reorganization for the Stable Charge Trapping in Amorphous Silicon Nitride.

Despite the widespread use of charge-trap flash (CTF) memory, the atomistic mechanism behind the exceptionally stable charge storage at the localized trap sites is still controversial. Herein, by combining first-principles calculations and orbital interaction analysis, a charge-dependent switchable chemical-bond reorganization is elucidated as the underpinning chemistry in the working mechanism of CTF. Especially, positively charged fourfold-coordinated nitrogen (dubbed N+ center), unappreciated until now, is the decisive component of the entire process; once an electron occupies this site, the N+ center disappears by breaking one N─Si bond, simultaneously forming a new Si─Si bond with a nearby Si atom which, in turn, creates fivefold coordinated Si. As a result, the electron is stored in a multi-center orbital belonging to multiple atoms including the newly formed Si─Si bond. It is also observed that hole trapping accompanies the creation of an N+ center by forming a new N─Si bond, which represents the reverse process. To further support and validate this model by means of core-level calculations, it is also shown that an N+ center's 1s core level is 1.0-2.5 eV deeper in energy than those of the threefold coordinated N atoms, in harmony with experimental X-ray photoelectron spectroscopy data.

Quantum Chemical Similarities of Bonding in Polyiodides and Phase-Change Materials.

Covalent chemical bonding beyond the two-center two-electron (2c-2e) bond is well-known for (inter)halogenic compounds, in particular, electron-rich multicenter (or hypervalent) bonding of the three-center four-electron (3c-4e) type to explain both their structure and stability. In the present work, we examine different solid-state polyiodides by combining both local orbital wave function and projected force constant analysis in order to numerically quantify the influence of multicenter (hypervalent) bonding based on periodic density functional theory (DFT) calculations. After linking our findings to established qualitative theories on multicenter bonding, particularly, Alcock's "secondary" bonding, we relate the bonding behavior in polyiodides to industrially relevant phase-change materials of the Ge-Sb-Te class, finding further evidence for the same underlying cause as regards chemical bonding in both material classes.

Interpenetration Phenomena via Anion Template Effects in Fe(II) and Co(II) Coordination Networks with a Bis-(1,2,4-triazole) Ligand.

Seven new coordination networks, [Fe(tbbt)3](BF4)2 (1), [Co(tbbt)3](BF4)2 (2), [Fe(tbbt)3](ClO4)2 (3), [Co(tbbt)3](ClO4)2 (4), [Fe(NCS)2(tbbt)2] (5), [Co(NCS)2(tbbt)2] (6), and [Fe(H2O)2(tbbt)2]Br2·2H2O (7), were synthesized with the linker 1,1'-(trans-2-butene-1,4-diyl)bis-1,2,4-triazole (tbbt) and structurally investigated. The structure of complexes 1-4 is composed of three interpenetrating, symmetry-related 3D networks. Each individual 3D network forms a primitive, nearly cubic lattice (pcu) with BF4- or ClO4- anions present in the interstitial spaces. The structure of compounds 5 and 6 is composed of two-dimensional sql layers, which are parallel to each other in the AB stacking type. These layers are interpenetrated by one-dimensional chains, both having the same formula unit, [M(NCS)2(tbbt)2] (M = Fe, Co). The structure of compound 7 consists of parallel, two-dimensional sql layers in the ABCD stacking type. The interpenetration in 1-6 is not controlled by π-π-interactions between the triazole rings or C=C bonds, as could have been expected, but by (triazole)C-H⋯F4B, C-H⋯O4Cl, and C-H⋯SCN anion hydrogen bonds, which suggests a template effect of the respective non-coordinated or coordinated anion for the interpenetration. In 7, the (triazole)C-H⋯Br anion interactions are supplemented by O-H⋯O and O-H⋯Br hydrogen bonds involving the aqua ligand and crystal water molecules. It is evident that the coordinated and non-coordinated anions play an essential role in the formation of the networks and guide the interpenetration. All iron(II) coordination networks are colorless, off-white to yellow-orange, and have the metal ions in the high-spin state down to 77 K. Compound 5 stays in the high spin state even at temperatures down to 10 K.

Synthesis and Characterization of Iron Bispyridine Bisdicyanamide, Fe[C5H5N]2[N(CN)2]2.

Fe[C5H5N]2[N(CN)2]2 (1) was synthesized from a reaction of stoichiometric amounts of NaN(CN)2 and FeCl2·4H2O in a methanol/pyridine solution. Single-crystal and powder diffraction show that 1 crystallizes in the monoclinic space group I2/m (no. 12), different from Mn[C5H5N]2[N(CN)2]2 (P21/c, no. 14) due to tilted pyridine rings, with a = 7.453(7) Å, b = 13.167(13) Å, c = 8.522(6) Å, ß = 114.98(6)° and Z = 2. ATR-IR, AAS, and CHN measurements confirm the presence of dicyanamide and pyridine. Thermogravimetric analysis shows that π-stacking interactions of the pyridine rings play an important role in structural stabilization. Based on DFT-optimized structures, a chemical bonding analysis was performed using a local-orbital framework by projection from a plane-wave basis. The resulting bond orders and atomic charges are in good agreement with the expectations based on the structure analysis. SQUID magnetic susceptibility measurements show a high-spin state FeII compound with predominantly antiferromagnetic exchange interactions at lower temperatures.

Embedded Localized Molecular-Orbital Representations for Periodic Wave Functions.

To more straightforwardly provide local chemical-bonding reasoning in crystalline matter, we introduce a new approach to generate a real-space analogue of periodic electronic structures using "exact" top-down frozen-density embedding calculations. Based on the obtained real-space electronic structure, we then construct localized molecular orbitals and evidence that our technique compares favorably against the commonly used Wannier method, both in terms of numerical efficiency and details of chemical bonding. The new method has been implemented into the LOBSTER software package and designed as a black-box approach, digesting any periodic electronic structure from the currently supported codes, i.e., VASP, Quantum ESPRESSO, and ABINIT.

POWTEX visits POWGEN.

The high-intensity time-of-flight (TOF) neutron diffractometer POWTEX for powder and texture analysis is currently being built prior to operation in the eastern guide hall of the research reactor FRM II at Garching close to Munich, Germany. Because of the world-wide 3He crisis in 2009, the authors promptly initiated the development of 3He-free detector alternatives that are tailor-made for the requirements of large-area diffractometers. Herein is reported the 2017 enterprise to operate one mounting unit of the final POWTEX detector on the neutron powder diffractometer POWGEN at the Spallation Neutron Source located at Oak Ridge National Laboratory, USA. As a result, presented here are the first angular- and wavelength-dependent data from the POWTEX detector, unfortunately damaged by a 50g shock but still operating, as well as the efforts made both to characterize the transport damage and to successfully recalibrate the voxel positions in order to yield nonetheless reliable measurements. Also described is the current data reduction process using the PowderReduceP2D algorithm implemented in Mantid [Arnold et al. (2014). Nucl. Instrum. Methods Phys. Res. A, 764, 156-166]. The final part of the data treatment chain, namely a novel multi-dimensional refinement using a modified version of the GSAS-II software suite [Toby & Von Dreele (2013). J. Appl. Cryst.46, 544-549], is compared with a standard data treatment of the same event data conventionally reduced as TOF diffraction patterns and refined with the unmodified version of GSAS-II. This involves both determining the instrumental resolution parameters using POWGEN's powdered diamond standard sample and the refinement of a friendly-user sample, BaZn(NCN)2. Although each structural parameter on its own looks similar upon comparing the conventional (1D) and multi-dimensional (2D) treatments, also in terms of precision, a closer view shows small but possibly significant differences. For example, the somewhat suspicious proximity of the a and b lattice parameters of BaZn(NCN)2 crystallizing in Pbca as resulting from the 1D refinement (0.008 Å) is five times less pronounced in the 2D refinement (0.038 Å). Similar features are found when comparing bond lengths and bond angles, e.g. the two N-C-N units are less differently bent in the 1D results (173 and 175°) than in the 2D results (167 and 173°). The results are of importance not only for POWTEX but also for other neutron TOF diffractometers with large-area detectors, like POWGEN at the SNS or the future DREAM beamline at the European Spallation Source.

Salt-flux growth of HoCuMg4 single crystals.

Polycrystalline samples of the magnesium-rich intermetallic compounds RECuMg4 (RE = Dy, Ho, Er, Tm) were synthesized by reaction of the elements in sealed tantalum ampoules heated in a high-frequency induction furnace. Phase purity of the RECuMg4 phases was ascertained by powder X-ray diffraction patterns. Well-shaped single crystals of HoCuMg4 could be grown in a NaCl/KCl salt flux and the crystal structure was refined from single crystal X-ray diffraction data: TbCuMg4 structure-type, space group Cmmm, a = 1361.4(2), b = 2039.3(4), c = 384.62(6) pm. The crystal structure of the RECuMg4 phases can be understood as a complex intergrowth variant of CsCl and AlB2 related slabs. The remarkable crystal chemical motif concerns the orthorhombically distorted bcc-like magnesium cubes with Mg-Mg distances ranging from 306 to 334 pm. At high temperatures DyCuMg4 and ErCuMg4 are Curie-Weiss paramagnets with paramagnetic Curie-Weiss temperatures of -15 K and -2 K for RE = Dy and Er, respectively. The effective magnetic moments, 10.66µB for RE = Dy and 9.65µB for RE = Er prove stable trivalent ground states for the rare earth cations. Magnetic susceptibility and heat capacity measurements reveal long-range antiferromagnetic ordering at low temperatures (<21 K). Whereas DyCuMg4 exhibits two subsequent antiferromagnetic transitions at TN = 21 and 7.9 K which successively remove half of the entropy of a doublet crystal field ground state of Dy, ErCuMg4 shows a single, possibly broadened, antiferromagnetic transition at 8.6 K. The successive antiferromagnetic transitions are discussed with respect to magnetic frustration in the tetrameric units present in the crystal structure.

The Myth of "Metavalency" in Phase-Change Materials.

Phase-change memory materials (PCMs) have unusual properties and important applications, and recent efforts to find improved materials have focused on their bonding mechanisms. "Metavalent bonding" or "metavalency," intermediate between "metallic" and "covalent" bonding and comprising single-electron bonds, has been proposed as a fundamentally new mechanism that is relevant both here and for halide perovskite materials. However, it is shown that PCMs, which violate the octet rule, have two types of covalent bond: two-center, two-electron (2c-2e) bonds, and electron-rich, multicenter bonds (3c-4e bonds, hyperbonds) involving lone-pair electrons. The latter have bond orders less than one and are examples of the century-old concept of "partial" bonds.

Structural Transformation in LnHS (Ln = La, Nd, Gd, and Er) with Coordination Change between an S-Centered Octahedron and a Trigonal Prism.

Lanthanide hydride chalcogenides LnHSe and LnHTe (Ln = lanthanides) crystallize in two polymorphs, 2H and 1H structures (ZrBeSi-type and filled-WC-type structures, respectively), but the chemical origin of the structural selection is unknown. Here, we have expanded the LnHCh (Ch = O, Se, and Te) family to include LnHS (Ln = La, Nd, Gd, and Er) using high-pressure synthesis. LnHS adopts the 2H structure for large Ln (La, Nd, and Gd) and the 1H structure for small Er. We compared the two polymorphs using anion-centered polyhedra and found that in the compounds with large ionicity, the 2H structure with ChLn6 octahedra is stabilized over the 1H structure with ChLn6 trigonal prisms due to relatively small electrostatic repulsion, supported by analysis of Madelung energy, crystal orbital Hamilton population (COHP), and density of energy (DOE).

Synthesis, crystal structures and semiconducting properties of new hexacyanidometallates.

We describe the synthesis, crystal structure and semiconducting properties of a number of hexacyanidometallates with the formula A2[MFe(CN)6]·xH2O (A = Na, K; M = Mg, Ca, Sr and Ba). All crystal structures were studied via single-crystal or powder X-ray diffraction. The unexpectedly low-symmetric structures in these ferrocyanides are described and contrasted with analogous transition-metal compounds which have been reported to be strictly or nearly cubic. The amount of crystal water in the structure for powder samples was determined by the thermogravimetric analysis (TGA), supported by IR and Raman spectroscopy. Electronic-structure calculations of K2[MgFe(CN)6] and K2[CaFe(CN)6] are compared with experimental UV-Vis measurements. The large band gaps by advanced theory indicate that the smaller experimental band gaps are due to surface effects of impurity states. Mott-Schottky curves of K2[MgFe(CN)6], K2[CaFe(CN)6] and K2[BaFe(CN)6]·3H2O exhibit positive slopes, which characterizes these compounds as n-type semiconductors.

Computational Design and Theoretical Properties of WC3N6, an H-Free Melaminate and Potential Multifunctional Material.

By means of first-principles theory, existence, synthetic conditions, and structural as well as physicochemical properties have been predicted for the first hydrogen-free melaminate salt of the composition WC3N6. We find at least two energetically favorable polymorphs adopting space groups P1 and P3, both of which are layer-like porous materials. In addition to sizable Madelung fields stabilizing saltlike WC3N6, the complex C3N66- anions are connected via perfectly optimized W-N bonds, forming WN5 in the P1 and WN6 coordination polyhedra in the P3 polymorphs. The band gaps of the P1 and P3 phases are HSE-predicted as 2.25 and 1.21 eV, respectively, significantly smaller than those of g-C3N4 and WO3. Moreover, both phases have suitable band-edge potentials that may provide sufficient driving force for photocatalytic water splitting; at least for the P1 phase, there is also a reasonable chance for reduced electron-hole recombination. In addition, the polymorphs's large optical absorption coefficients should greatly enhance the photocatalytic performance. WC3N6 defines a new class of compounds and has unique structural characteristics, mirrored from its electrical and optical properties, and it should provide another chemical path for preparing efficient photocatalysts and optoelectronic devices.

On the atomistic origin of the polymorphism and the dielectric physical properties of beryllium oxide.

Polymorphic beryllium oxide has been theoretically investigated from first principles as regards orbital occupancies, chemical bonding, polarization, as well as dielectric properties. By means of Crystal-Orbital Bond Index (COBI) analysis, the important role of the 2p orbitals on beryllium has been elucidated, in particular in terms of the correlation between polarization and beryllium-atom displacement, including the impact of the latter on the covalency of the BeO bond. In addition, several structural possibilities for a Bex Mg1-x O solid solution have been investigated for a Be content between 6% and 22%; for those, dynamically stable structures have been found, displaying large polarization values, more covalent BeO bonds, and a tendency for tetrahedral Be coordination. The dynamically unstable structures, however, resemble rock-salt BeO in their local structural properties around the Be atom. High dielectric constants and band gaps indicating insulating behavior have been found for those.

Li2[SeC2Se]·2NH3: A Crystalline Ammoniate with a -Se-C≡C-Se- Dianion.

Reaction of Li2C2 with elemental selenium in a molar ratio of 1:2 in liquid ammonia led to the formation of the ammoniate Li2[SeC2Se]·2NH3. Its crystal structure was solved and refined from high-resolution synchrotron powder diffraction data (P21/c, Z = 4). It contains the -Se-C≡C-Se- anion, unprecedented in a crystalline material, whose existence was corroborated by IR/Raman spectra and electronic-structure theory, showing an almost perfect agreement with calculated spectra. Elaborated magnetic-bottle and velocity-map imaging photoelectron spectroscopic investigations show that the -Se-C≡C-Se• radical anion can be transferred to the gas phase, where it was analyzed by NIPE (Negative Ion Photoelectron) and VMI (Velocity-Map Imaging) spectra, which correlate nicely with simulated spectra based on 2Πu → 3Σg- and 2Πu → 1Σg+ transitions including spin-orbit couplings.

Synthetic Engineering in Na2MSn2(NCN)6 (M = Mn, Fe, Co, and Ni) Based on Electronic Structure Theory.

Quaternary transition metal cyanamides Na2MSn2(NCN)6 with M = Mn, Fe, Co, and Ni were prepared via solid-state metathesis reactions between Na2Sn(NCN)3 and binary transition metal fluorides MF2 in a 2:1 molar ratio. All phases crystallize isotypically in [NiAs]-derived structures (P3̅1m) with inter- and intra-layer cation ordering over the octahedral sites. This leads to a highly asymmetric coordination of the NCN2- anion, resulting in a strong degree of cyanamide character, which is confirmed via IR measurements. Intriguingly, the optical properties of Na2MSn2(NCN)6 phases change markedly as the nature of the divalent transition metal is varied, and UV-vis measurements evidence a band gap reduction from Mn (3.43 eV) via Fe (1.90 eV) to Co (1.75 eV), which broadly mirrors the DFT+U calculated energetic interval from the Fermi level to the unoccupied 3d states. Mott-Schottky analysis then goes on to characterize Na2FeSn2(NCN)6 and Na2CoSn2(NCN)6 as n-type semiconductors with flat-band potentials of 0.46 and -0.24 eV, respectively, vs RHE. This study demonstrates the utility of transition metal substitutions, within a flexible cyanamide framework, to electronically tune this growing family of pseudo-oxides.

Automated Bonding Analysis with Crystal Orbital Hamilton Populations.

Invited for this month's cover are researchers from Bundesanstalt für Materialforschung und -prüfung (Federal Institute for Materials Research and Testing) in Germany, Friedrich Schiller University Jena, Université catholique de Louvain, University of Oregon, Science & Technology Facilities Council, RWTH Aachen University, Hoffmann Institute of Advanced Materials, and Dartmouth College. The cover picture shows a workflow for automatic bonding analysis with Python tools (green python). The bonding analysis itself is performed with the program LOBSTER (red lobster). The starting point is a crystal structure, and the results are automatic assessments of the bonding situation based on Crystal Orbital Hamilton Populations (COHP), including automatic plots and text outputs. Coordination environments and charges are also assessed. More information can be found in the Research Article by J. George, G. Hautier, and co-workers.

Between Elemental Match and Mismatch: From K12 Ge3.5 Sb6 to Salts of (Ge2 Sb2 )2- , (Ge4 Sb12 )4- , and (Ge4 Sb14 )4.

The solid mixture "K2 GeSb" was shown to comprise single-crystalline K12 Ge3.5 Sb6 (1), a double salt of K5 [GeSb3 ] with carbonate-like [GeSb3 ]5- anions, and the metallic Zintl phase K2 Ge1.5 . Extraction of 1 with ethane-1,2-diamine in the presence of crypt-222 afforded [K(crypt-222)]+ salts of several novel binary Zintl anions: (Ge2 Sb2 )2- (in 2), (Ge4 Sb12 )4- (in 3), and in the presence of [AuMePPh3 ] also (Ge4 Sb14 )4- (in 4). The anion in 2 represents a predicted, yet heretofore missing pseudo-tetrahedral anion. 4 comprises a cluster analogous to (Ge4 Bi14 )4- and (Ga2 Bi16 )4- , and thus one of the most Sb-rich binary p-block anions. The unprecedented cluster topology in 3 can be viewed as a defect-version of the one in 4 upon following a "dead end" of cluster growth. The findings indicate that Ge and Sb atoms are at the border of a well-matching and a mismatch elemental combination. We discuss the syntheses, the geometric structures, and the electronic structures of the new compounds.

Automated Bonding Analysis with Crystal Orbital Hamilton Populations.

Understanding crystalline structures based on their chemical bonding is growing in importance. In this context, chemical bonding can be studied with the Crystal Orbital Hamilton Population (COHP), allowing for quantifying interatomic bond strength. Here we present a new set of tools to automate the calculation of COHP and analyze the results. We use the program packages VASP and LOBSTER, and the Python packages atomate and pymatgen. The analysis produced by our tools includes plots, a textual description, and key data in a machine-readable format. To illustrate those capabilities, we have selected simple test compounds (NaCl, GaN), the oxynitrides BaTaO2 N, CaTaO2 N, and SrTaO2 N, and the thermoelectric material Yb14 Mn1 Sb11 . We show correlations between bond strengths and stabilities in the oxynitrides and the influence of the Mn-Sb bonds on the magnetism in Yb14 Mn1 Sb11 . Our contribution enables high-throughput bonding analysis and will facilitate the use of bonding information for machine learning studies.

Chemical Reactions and Phase Stabilities in the Si-Te System at High Pressures and High Temperatures.

Chemical reactions and phase stabilities in the Si-Te system at high pressures were explored using in situ angle-dispersive synchrotron powder diffraction in a large-volume multianvil press together with density functional theory-based calculations. Cubic and rhombohedrally distorted clathrates, with the general formula Te8@(Si38Te8) and wide compositional range, preceded by a hexagonal phase with the composition Si0.14Te, are formed for different mixtures of Si and Te as starting materials. Si0.14Te, with the structural formula Te2(Te0.74Si0.26)3(Te0.94Si0.06)3, is the very first chalcogenide with the Mn5Si3-type structure. Silicon sesquitelluride α-Si2Te3 decomposes into a mixture of phases that includes the clathrate and hexagonal phases at high pressures and high temperatures. The higher the pressure, the lower the temperature for the two phases to occur. Regardless of the starting compositions, only the clathrate is quenched to atmospheric conditions, while the hexagonal phase amorphizes on decompression. The rhombohedral clathrates Te8@(Si38Te8) form on quenching of the cubic phases to ambient conditions. There is a high degree of interchangeability of Si and Te not only in the clathrates but also in the Mn5Si3-type structure. The theoretical calculations of enthalpies indicate that the reported decomposition of α-Si2Te3 is energetically favorable over its transformation to another polymorph of the A2X3 type at extreme conditions.