High-Throughput Identification of Electrides from All Known Inorganic
Total Page:16
File Type:pdf, Size:1020Kb
Article Cite This: Chem. Mater. XXXX, XXX, XXX−XXX pubs.acs.org/cm High-Throughput Identification of Electrides from All Known Inorganic Materials Lee A. Burton, Francesco Ricci, Wei Chen, Gian-Marco Rignanese, and Geoffroy Hautier* Institute of Condensed Matter and Nanoscience, UniversitéCatholique de Louvain, 1348 Louvain-la-Neuve, Belgium *S Supporting Information ABSTRACT: In this paper, we present the results of a large- scale, high-throughput computational search for electrides among all known inorganic materials. Analyzing a database of density functional theory results on more than 60 000 compounds, we identify 65 new electride candidates and recover 4 already known. We report on all these candidates and discuss the structural and chemical factors leading to electride formation. Among these candidates, our work identifies the first partially filled 3d transition-metal-containing fi electrides Ba3CrN3 and Sr3CrN3, an unexpected nding that contravenes conventional chemistry. ■ INTRODUCTION large energy of a unbound electron.15 For example, the 2D 2+ 3− Electrides are rare ionic compounds in which an electron does electride Ca2N has two Ca cations and one N anion not occupy an atomic orbital but rather acts as an anion. Such (according to commonly held concepts of oxidation state). an electron is expected to behave differently from those The electropositive Ca cations wish to donate four electrons, occupying the valence state of standard materials, making but the anion can only accept three, leading to an excess electrides desirable as electron emitters,1 nonlinear optical electron occupying free space in the crystal structure, thereby switches,2 superconductors,3 battery anodes,4 and catalysts for satisfying all three criteria. The need for these three factors to combine adequately makes electrides a rare occurrence among applications ranging from compound synthesis to CO2 splitting.5,6 solid compounds with only a handful electrides known so far. Chemically independent electrons were first reported from Quantum-mechanical calculations have been essential in the salts of metal−ammonia solutions.7 Solid-state organic quest for new electrides. They are commonly used in electrides were reported in 1990 but were unstable at room conjunction with experimental characterization (e.g., catalytic temperature even in inert atmosphere.8 In 2012, a break- activity or transport measurements) to demonstrate the through came with the first inorganic electride stable at room electride nature of a material. By observing the localization temperature and in air: Ca Al O in the mayenite structure.9 of electrons resulting from computations in the density 12 14 32 functional theory (DFT) framework, one can determine if Downloaded via UNIV CATHOLIQUE DE LOUVAIN on October 22, 2018 at 14:10:08 (UTC). Now, various stable electrides are known. They are often See https://pubs.acs.org/sharingguidelines for options on how to legitimately share published articles. fi the electrons occupy space away from a nuclei, the definition of classi ed according to the localization of the anionic electron 16,17 in the lattice: 0-dimensional (0D)10 when it resides in a cavity, an electride. Ab initio techniques have been used 11 extensively in the past to predict electride-like behavior in 1-dimensional (1D) when it is enclosed in a channel, and 2- 18−21 dimensional (2D)12 when it is confined to a layer. These high-pressure elements. However, no large-scale, broad compounds are now being tested in multiple industrial computationally driven search for electrides has yet been 13 performed among known, stable, inorganic compounds at applications from NH3 catalysis, to light emitting diode (LED) devices.14 As such, electride chemistry has become an ambient pressure. active field of solid-state chemistry and materials science with Here, we search for potential electrides among a large large efforts directed toward their discovery, synthesis, and database containing DFT electronic structure data on materials originating mainly from those experimentally reported in the characterization. 22 Several simple rules have been used to identify possible Inorganic Crystal Structure Database (ICSD). From this electride candidates. First, the structure must contain free database of more than 60 000 compounds, we identify 65 new space where the anionic electron can reside,8 be this a cavity, a potential electrides and recover 4 already known, without channel, or a layer. Second, a compound must have an excess electron to contribute to the lattice, as determined by the sum Received: June 15, 2018 of oxidation states of the components.15 Finally, the compound Revised: October 9, 2018 must contain a strongly donating cation to offset the naturally Published: October 10, 2018 © XXXX American Chemical Society A DOI: 10.1021/acs.chemmater.8b02526 Chem. Mater. XXXX, XXX, XXX−XXX Chemistry of Materials Article relying on any of the design principles associated with their NdGa, CaSi, SrSn, BaGe, SrGe, CaGe, BaSi, SrSi, and BaSn, usual discovery. We discuss the 3 most popular groups of which are orthorhombic with the space group Cmcm (No. 63). candidates, which incidentally form 0D, 1D, and 2D electrides. The structures consist of bonded atomic bilayers sandwiched Furthermore, we highlight what we believe to be the only between 2D free electron layers (see Figure 1) and are referred known instances of electrides containing a redox active fi element: Ba3CrN3 and Sr3CrN3 (signi cant as one would expect a redox active transition metal to accept the anionic electron). Finally, we re-evaluate the previously proposed guiding principles for electride formation in light of our large- scale screening. ■ METHODS Our computational screening began with the Materials Project,22 a database of 69 640 materials computed using the Vienna ab initio Simulation Package (VASP) using the GGA-PBE implementation of DFT.23,24 Of these, we exclude materials of just one component, i.e., elements, and those with a calculated energy above the convex hull of more than 24 meV per atom, which has been shown to be an accurate measure of stability in compounds.25 We focused on materials with Figure 1. Electron density around the Fermi level for the CrB class the Fermi level (EF) crossing several bands of the compounds that have calculated band structures, i.e., metals. The first screening step candidates clearly indicative of a 2D electride material (NdGa in this was using the projections on atomic orbitals for the eigenstates case). around EF in a range of 0.01 eV. If these projections summed to less than 50% of the total electrons present, we expect a high chance for to as CrB-type structures in the original reports. Of these 11, 8 the electrons to occupy nonatomic sites and, thus, be an electride. In a are II−IV compounds, complying with the principle of second step, we performed a Bader charge analysis to identify electron including an electropositive element (group II) in the distributions with a relative charge of at least 0.1 located 2.2 Å or 26−28 compound. However, the rule of an excess electron from the more away from the nearest atom, within 0.2 eV of the EF for a 2000 k points per Å−3 resolution calculation,29 similar to alternate oxidation states does not apply. To the knowledge of the work in the field.30,31 The parameters were chosen to return as few authors, these compounds have never been considered candidates as possible while still recovering the four known electrides electrides but show 2D behavior similar to that of the well- 32 33 16 34 Ca2N, Y2C, Y5Si3, and LaH2, identifying 69 nonelemental studied Ca2N. “ ” candidates in total. We note that only as-is electrides will be found Of the 69 systems, 9 adopt the so-called Mn5Si3 structure fi through our approach and that we do not expect to nd mayenite, for type, including the known electride Y5Si3, which is reported to instance, which requires a chemical treatment (oxygen deintercala- 16 35 be chemically stable even in the presence of moisture. These tion) to produce its electride form. are Ba As ,CaSb ,NdGe ,SrAs ,YSi ,SrSb ,SrBi , The quasiparticle self-consistent GW (QSGW) calculations were 5 3 5 3 5 3 5 3 5 3 5 3 5 3 36 Ba5Bi3, and Ba5Sb3, all of which are hexagonal with the space performed using ABINIT. The dielectric function and quasiparticle − self-energies are evaluated on a 4 × 4 × 6 k-point mesh and with 1000 group P63/mcm (No. 193). Of the 9 compounds, 7 are II V bands (see Supporting Information for the convergence study). The compounds. In the unit cells, 3 of 5 cations sit at the center of energy cutoff for the dielectric function is 160 eV. The semicore sp face-sharing octahedra, while the remaining 2 coordinate electrons are included among valence for the Ba, Sr, and Cr atoms.37 toward a 1D channel populated by the anionic electron, see In addition to all the valence states, we include 80 unoccupied states Figure 2. for which the wave functions and eigenvalues are iterated to self- consistency. The virtue of QSGW (compared to one-shot G0W0)is that when self-consistency is reached, all remnants of the DFT starting point are eliminated. ■ RESULTS Our screening begins with the Materials Project database, which includes DFT electronic structure data for a large set of known materials. By analyzing the charge density around the Fermi level for all stable metallic materials within this database, we can detect the compounds presenting the off-nuclei electron localization characteristic of an electride (see Methods). Our search identifies 69 compounds with 4 32 33 previously reported electrides among them (Ca2N, Y2C, Figure 2. Structure and electron density channel (perpendicular to Y Si ,16 and LaH 34). We stress that, contrary to previous 5 3 2 the page) for the Mn5Si3 compounds (Ba5As3 in this case).