Computationally Predicted High-Throughput Free-Energy Phase Diagrams for the Discovery of Solid-State Hydrogen Storage Reactions.

The design of multinary solid-state material systems that undergo reversible phase changes via changes in temperature and pressure provides a potential means of safely storing hydrogen. However, fully mapping the stabilities of known or newly targeted compounds relative to competing phases at reaction conditions has previously required many stringent experiments or computationally demanding calculations of each compound's change in Gibbs energy with respect to temperature, G(T). In this work, we have extended the approach of constructing chemical potential phase diagrams based on ΔGf(T) to enable the analysis of phase stability at non-zero temperatures. We first performed density functional theory calculations to compute the formation enthalpies of binary, ternary, and quaternary compounds within several compositional spaces of current interest for solid-state hydrogen storage. Temperature effects on solid compound stability were then accounted for using our recently introduced machine learned descriptor for the temperature-dependent contribution Gδ(T) to the Gibbs energy G(T). From these Gibbs energies, we evaluated each compound's stability relative to competing compounds over a wide range of conditions and show using chemical potential and composition phase diagrams that the predicted stable phases and H2 release reactions are consistent with experimental observations. This demonstrates that our approach rapidly computes the thermochemistry of hydrogen release reactions for compounds at sufficiently high accuracy relative to experiment to provide a powerful framework for analyzing hydrogen storage materials. This framework based on G(T) enables the accelerated discovery of active materials for a variety of technologies that rely on solid-state reactions involving these materials.

Stabilizing Hydrogen Adsorption through Theory-Guided Chalcogen Substitution in Chevrel-Phase Mo6X8 (X=S, Se, Te) Electrocatalysts.

In this work, we implement a facile microwave-assisted synthesis method to yield three binary Chevrel-Phase chalcogenides (Mo6X8; X = S, Se, Te) and investigate the effect of increasing chalcogen electronegativity on hydrogen evolution catalytic activity. Density functional theory predictions indicate that increasing chalcogen electronegativity in these materials will yield a favorable electronic structure for proton reduction. This is confirmed experimentally via X-ray absorption spectroscopy as well as traditional electrochemical analysis. We have identified that increasing the electronegativity of X in Mo6X8 increases the hydrogen adsorption strength owing to a favorable shift in the p-band position as well as an increase in the Lewis basicity of the chalcogen, thereby improving hydrogen evolution reaction energetics. We find that Mo6S8 exhibits the highest hydrogen evolution activity of the Mo6X8 series of catalysts, requiring an overpotential of 321 mV to achieve a current density of 10 mA cm-2ECSA, a Tafel slope of 74 mV per decade, and an exchange current density of 6.01 × 10-4 mA cm-2ECSA. Agreement between theory and experiment in this work indicates that the compositionally tunable Chevrel-Phase chalcogenide family is a promising framework for which electronic structure can be predictably modified to improve catalytic small-molecule reduction reactivity.

Inorganic Halide Double Perovskites with Optoelectronic Properties Modulated by Sublattice Mixing.

All-inorganic halide double perovskites have emerged as a promising class of materials that are potentially more stable and less toxic than lead-containing hybrid organic-inorganic perovskite optoelectronic materials. In this work, 311 cesium chloride double perovskites (Cs2BB'Cl6) were selected from a set of 903 compounds as likely being stable on the basis of a statistically learned tolerance factor (τ) for perovskite stability. First-principles calculations on these 311 double perovskites were then performed to assess their stability and identify candidates with band gaps appropriate for optoelectronic applications. We predict that 261 of the 311 Cs2BB'Cl6 compounds are likely synthesizable on the basis of a thermodynamic analysis of their decomposition to competing compounds (decomposition enthalpy <0.05 eV/atom). Of these 261 likely synthesizable compounds, 47 contain no toxic elements and have direct or nearly direct (within 100 meV) band gaps between 1 and 3 eV, as computed with hybrid density functional theory (HSE06). Within this set, we identify the triple-alkali perovskites Cs2[Alk]+[TM]3+Cl6, where Alk is a group 1 alkali cation and TM is a transition-metal cation, as a class of Cs2BB'Cl6 double perovskites with remarkable optical properties, including large and tunable exciton binding energies as computed by the GW-Bethe-Salpeter equation (GW-BSE) method. We attribute the unusual electronic structure of these compounds to the mixing of the Alk-Cl and TM-Cl sublattices, leading to materials with small band gaps, large exciton binding energies, and absorption spectra that are strongly influenced by the identity of the transition metal. The role of the double-perovskite structure in enabling these unique properties is probed through an analysis of the electronic structures and chemical bonding of these compounds in comparison with other transition-metal and alkali transition-metal halides.

High-Throughput Experimental Study of Wurtzite Mn1-x Zn x O Alloys for Water Splitting Applications.

We used high-throughput experimental screening methods to unveil the physical and chemical properties of Mn1-x Zn x O wurtzite alloys and identify their appropriate composition for effective water splitting application. The Mn1-x Zn x O thin films were synthesized using combinatorial pulsed laser deposition, permitting for characterization of a wide range of compositions with x varying from 0 to 1. The solubility limit of ZnO in MnO was determined using the disappearing phase method from X-ray diffraction and X-ray fluorescence data and found to increase with decreasing substrate temperature due to kinetic limitations of the thin-film growth at relatively low temperature. Optical measurements indicate the strong reduction of the optical band gap down to 2.1 eV at x = 0.5 associated with the rock salt-to-wurtzite structural transition in Mn1-x Zn x O alloys. Transmission electron microscopy results show evidence of a homogeneous wurtzite alloy system for a broad range of Mn1-x Zn x O compositions above x = 0.4. The wurtzite Mn1-x ZnxO samples with the 0.4 < x < 0.6 range were studied as anodes for photoelectrochemical water splitting, with a maximum current density of 340 µA cm-2 for 673 nm-thick films. These Mn1-x Zn x O films were stable in pH = 10, showing no evidence of photocorrosion or degradation after 24 h under water oxidation conditions. Doping Mn1-x Zn x O materials with Ga dramatically increases the electrical conductivity of Mn1-x Zn x O up to ∼1.9 S/cm for x = 0.48, but these doped samples are not active in water splitting. Mott-Schottky and UPS/XPS measurements show that the presence of dopant atoms reduces the space charge region and increases the number of mid-gap surface states. Overall, this study demonstrates that Mn1-x Zn x O alloys hold promise for photoelectrochemical water splitting, which could be enhanced with further tailoring of their electronic properties.

A map of the inorganic ternary metal nitrides.

Exploratory synthesis in new chemical spaces is the essence of solid-state chemistry. However, uncharted chemical spaces can be difficult to navigate, especially when materials synthesis is challenging. Nitrides represent one such space, where stringent synthesis constraints have limited the exploration of this important class of functional materials. Here, we employ a suite of computational materials discovery and informatics tools to construct a large stability map of the inorganic ternary metal nitrides. Our map clusters the ternary nitrides into chemical families with distinct stability and metastability, and highlights hundreds of promising new ternary nitride spaces for experimental investigation-from which we experimentally realized seven new Zn- and Mg-based ternary nitrides. By extracting the mixed metallicity, ionicity and covalency of solid-state bonding from the density functional theory (DFT)-computed electron density, we reveal the complex interplay between chemistry, composition and electronic structure in governing large-scale stability trends in ternary nitride materials.

Rational Design of Efficient Amine Reductant Initiators for Amine-Peroxide Redox Polymerization.

Amine-peroxide redox polymerization (APRP) has been highly prevalent in industrial and medical applications since the 1950s, yet the initiation mechanism of this radical polymerization process is poorly understood so that innovations in the field are largely empirically driven and incremental. Through a combination of computational prediction and experimental analysis, we elucidate the mechanism of this important redox reaction between amines and benzoyl peroxide for the ambient production of initiating radicals. Our calculations show that APRP proceeds through SN2 attack by the amine on the peroxide but that homolysis of the resulting intermediate is the rate-determining step. We demonstrate a correlation between the computationally predicted initiating rate and the experimentally measured polymerization rate with an R2 = 0.80. The new mechanistic understanding was then applied to computationally predict amine reductant initiators with faster initiating kinetics. This led to our discovery of N-(4-methoxyphenyl)pyrrolidine (MPP) as amine reductant, which we confirmed significantly outperforms current state-of-the-art tertiary aromatic amines by ∼20-fold, making it the most efficient amine-peroxide redox initiator to date. The application of amines with superior kinetics such as MPP in APRP could greatly accelerate existing industrial processes, facilitate new industrial manufacturing methods, and improve biocompatibility in biomedical applications conducted with reduced initiator concentrations yet higher overall efficiency.

High-Throughput Equilibrium Analysis of Active Materials for Solar Thermochemical Ammonia Synthesis.

Solar thermochemical ammonia (NH3) synthesis (STAS) is a potential route to produce NH3 from air, water, and concentrated sunlight. This process involves the chemical looping of an active redox pair that cycles between a metal nitride and its complementary metal oxide to yield NH3. To identify promising candidates for STAS cycles, we performed a high-throughput thermodynamic screening of 1,148 metal nitride/metal oxide pairs. This data-driven screening was based on Gibbs energies of crystalline metal oxides and nitrides at elevated temperatures, G(T), calculated using a recently introduced statistically learned descriptor and 0 K DFT formation energies tabulated in the Materials Project database. Using these predicted G(T) values, we assessed the viability of each of the STAS reactions-hydrolysis of the metal nitride, reduction of the metal oxide, and nitrogen fixation to reform the metal nitride-and analyzed a revised cycle that directly converts between metal oxides and nitrides, which alters the thermodynamics of the STAS cycle. For all 1148 redox pairs analyzed and each of the STAS-relevant reactions, we implemented a Gibbs energy minimization scheme to predict the equilibrium composition and yields of the STAS cycle, which reveals new active materials based on B, V, Fe, and Ce that warrant further investigation for their potential to mediate the STAS cycle. This work details a high-throughput approach to assessing the relevant temperature-dependent thermodynamics of thermochemical redox processes that leverages the wealth of publicly available temperature-independent thermodynamic data calculated using DFT. This approach is readily adaptable to discovering optimal materials for targeted thermochemical applications and enabling the predictive synthesis of new compounds using thermally controlled solid-state reactions.

Physical descriptor for the Gibbs energy of inorganic crystalline solids and temperature-dependent materials chemistry.

The Gibbs energy, G, determines the equilibrium conditions of chemical reactions and materials stability. Despite this fundamental and ubiquitous role, G has been tabulated for only a small fraction of known inorganic compounds, impeding a comprehensive perspective on the effects of temperature and composition on materials stability and synthesizability. Here, we use the SISSO (sure independence screening and sparsifying operator) approach to identify a simple and accurate descriptor to predict G for stoichiometric inorganic compounds with ~50 meV atom-1 (~1 kcal mol-1) resolution, and with minimal computational cost, for temperatures ranging from 300-1800 K. We then apply this descriptor to ~30,000 known materials curated from the Inorganic Crystal Structure Database (ICSD). Using the resulting predicted thermochemical data, we generate thousands of temperature-dependent phase diagrams to provide insights into the effects of temperature and composition on materials synthesizability and stability and to establish the temperature-dependent scale of metastability for inorganic compounds.

Novel phase diagram behavior and materials design in heterostructural semiconductor alloys.

Structure and composition control the behavior of materials. Isostructural alloying is historically an extremely successful approach for tuning materials properties, but it is often limited by binodal and spinodal decomposition, which correspond to the thermodynamic solubility limit and the stability against composition fluctuations, respectively. We show that heterostructural alloys can exhibit a markedly increased range of metastable alloy compositions between the binodal and spinodal lines, thereby opening up a vast phase space for novel homogeneous single-phase alloys. We distinguish two types of heterostructural alloys, that is, those between commensurate and incommensurate phases. Because of the structural transition around the critical composition, the properties change in a highly nonlinear or even discontinuous fashion, providing a mechanism for materials design that does not exist in conventional isostructural alloys. The novel phase diagram behavior follows from standard alloy models using mixing enthalpies from first-principles calculations. Thin-film deposition demonstrates the viability of the synthesis of these metastable single-phase domains and validates the computationally predicted phase separation mechanism above the upper temperature bound of the nonequilibrium single-phase region.

Dihydropteridine/Pteridine as a 2H+/2e- Redox Mediator for the Reduction of CO2 to Methanol: A Computational Study.

Conflicting experimental results for the electrocatalytic reduction of CO2 to CH3OH on a glassy carbon electrode by the 6,7-dimethyl-4-hydroxy-2-mercaptopteridine have been recently reported [ J. Am. Chem. Soc. 2014 , 136 , 14007 - 14010 , J. Am. Chem. Soc. 2016 , 138 , 1017 - 1021 ]. In this connection, we have used computational chemistry to examine the issue of this molecule's ability to act as a hydride donor to reduce CO2. We first determined that the most thermodynamically stable tautomer of this aqueous compound is its oxothione form, termed here PTE. It is argued that this species electrochemically undergoes concerted 2H+/2e- transfers to first form the kinetic product 5,8-dihydropteridine, followed by acid-catalyzed tautomerization to the thermodynamically more stable 7,8-dihydropteridine PTEH2. While the overall conversion of CO2 to CH3OH by three successive hydride and proton transfers from this most stable tautomer is computed to be exergonic by 5.1 kcal/mol, we predict high activation free energies (ΔG‡HT) of 29.0 and 29.7 kcal/mol for the homogeneous reductions of CO2 and its intermediary formic acid product by PTE/PTEH2, respectively. These high barriers imply that PTE/PTEH2 is unable, by this mechanism, to homogeneously reduce CO2 on a time scale of hours at room temperature.

Mechanisms of LiCoO2 Cathode Degradation by Reaction with HF and Protection by Thin Oxide Coatings.

Reactions of HF with uncoated and Al and Zn oxide-coated surfaces of LiCoO2 cathodes were studied using density functional theory. Cathode degradation caused by reaction of HF with the hydroxylated (101̅4) LiCoO2 surface is dominated by formation of H2O and a LiF precipitate via a barrierless reaction that is exothermic by 1.53 eV. We present a detailed mechanism where HF reacts at the alumina coating to create a partially fluorinated alumina surface rather than forming AlF3 and H2O and thus alumina films reduce cathode degradation by scavenging HF and avoiding H2O formation. In contrast, we find that HF etches monolayer zinc oxide coatings, which thus fail to prevent capacity fading. However, thicker zinc oxide films mitigate capacity loss by reacting with HF to form a partially fluorinated zinc oxide surface. Metal oxide coatings that react with HF to form hydroxyl groups over H2O, like the alumina monolayer, will significantly reduce cathode degradation.

Intrinsic Material Properties Dictating Oxygen Vacancy Formation Energetics in Metal Oxides.

Oxygen vacancies (V(O)) in oxides are extensively used to manipulate vital material properties. Although methods to predict defect formation energies have advanced significantly, an understanding of the intrinsic material properties that govern defect energetics lags. We use first-principles calculations to study the connection between intrinsic (bulk) material properties and the energy to form a single, charge neutral oxygen vacancy (E(V)). We investigate 45 binary and ternary oxides and find that a simple model which combines (i) the oxide enthalpy of formation (ΔH(f)), (ii) the midgap energy relative to the O 2p band center (E(O 2p) + (1/2)E(g)), and (iii) atomic electronegativities reproduces calculated E(V) within ∼0.2 eV. This result provides both valuable insights into the key properties influencing E(V) and a direct method to predict E(V). We then predict the E(V) of ∼1800 oxides and validate the predictive nature of our approach against direct defect calculations for a subset of 18 randomly selected materials.

Solvent Control of Surface Plasmon-Mediated Chemical Deposition of Au Nanoparticles from Alkylgold Phosphine Complexes.

Bottom-up approaches to nanofabrication are of great interest because they can enable structural control while minimizing material waste and fabrication time. One new bottom-up nanofabrication method involves excitation of the surface plasmon resonance (SPR) of a Ag surface to drive deposition of sub-15 nm Au nanoparticles from MeAuPPh3. In this work we used density functional theory to investigate the role of the PPh3 ligands of the Au precursor and the effect of adsorbed solvent on the deposition process, and to elucidate the mechanism of Au nanoparticle deposition. In the absence of solvent, the calculated barrier to MeAuPPh3 dissociation on the bare surface is <20 kcal/mol, making it facile at room temperature. Once adsorbed on the surface, neighboring MeAu fragments undergo ethane elimination to produce Au adatoms that cluster into Au nanoparticles. However, if the sample is immersed in benzene, we predict that the monolayer of adsorbed solvent blocks the adsorption of MeAuPPh3 onto the Ag surface because the PPh3 ligand is large compared to the size of the exposed surface between adsorbed benzenes. Instead, the Au-P bond of MeAuPPh3 dissociates in solution (Ea = 38.5 kcal/mol) in the plasmon heated near-surface region followed by the adsorption of the MeAu fragment on Ag in the interstitial space of the benzene monolayer. The adsorbed benzene forces the Au precursor to react through the higher energy path of dissociation in solution rather than dissociatively adsorbing onto the bare surface. This requires a higher temperature if the reaction is to proceed at a reasonable rate and enables the control of deposition by the light induced SPR heating of the surface and nearby solution.

Catalytic Reduction of CO2 by Renewable Organohydrides.

Dihydropyridines are renewable organohydride reducing agents for the catalytic reduction of CO2 to MeOH. Here we discuss various aspects of this important reduction. A centerpiece, which illustrates various general principles, is our theoretical catalytic mechanism for CO2 reduction by successive hydride transfers (HTs) and proton transfers (PTs) from the dihydropyridine PyH2 obtained by 1H(+)/1e(-)/1H(+)/1e(-) reductions of pyridine. The Py/PyH2 redox couple is analogous to NADP(+)/NADPH in that both are driven to effect HTs by rearomatization. High-energy radical intermediates and their associated high barriers/overpotentials are avoided because HT involves a 2e(-) reduction. A HT-PT sequence dictates that the reduced intermediates be protonated prior to further reduction for ultimate MeOH formation; these protonations are aided by biased cathodes that significantly lower the local pH. In contrast, cathodes that efficiently reduce H(+) such as Pt and Pd produce H2 and create a high interfacial pH, both obstructing dihydropyridine production and formate protonation and thus ultimately CO2 reduction by HTPTs. The role of water molecule proton relays is discussed. Finally, we suggest future CO2 reduction strategies by organic (photo)catalysts.

Reduction of CO2 to methanol catalyzed by a biomimetic organo-hydride produced from pyridine.

We use quantum chemical calculations to elucidate a viable mechanism for pyridine-catalyzed reduction of CO2 to methanol involving homogeneous catalytic steps. The first phase of the catalytic cycle involves generation of the key catalytic agent, 1,2-dihydropyridine (PyH2). First, pyridine (Py) undergoes a H(+) transfer (PT) to form pyridinium (PyH(+)), followed by an e(-) transfer (ET) to produce pyridinium radical (PyH(0)). Examples of systems to effect this ET to populate PyH(+)'s LUMO (E(0)(calc) ∼ -1.3 V vs SCE) to form the solution phase PyH(0) via highly reducing electrons include the photoelectrochemical p-GaP system (E(CBM) ∼ -1.5 V vs SCE at pH 5) and the photochemical [Ru(phen)3](2+)/ascorbate system. We predict that PyH(0) undergoes further PT-ET steps to form the key closed-shell, dearomatized (PyH2) species (with the PT capable of being assisted by a negatively biased cathode). Our proposed sequential PT-ET-PT-ET mechanism for transforming Py into PyH2 is analogous to that described in the formation of related dihydropyridines. Because it is driven by its proclivity to regain aromaticity, PyH2 is a potent recyclable organo-hydride donor that mimics important aspects of the role of NADPH in the formation of C-H bonds in the photosynthetic CO2 reduction process. In particular, in the second phase of the catalytic cycle, which involves three separate reduction steps, we predict that the PyH2/Py redox couple is kinetically and thermodynamically competent in catalytically effecting hydride and proton transfers (the latter often mediated by a proton relay chain) to CO2 and its two succeeding intermediates, namely, formic acid and formaldehyde, to ultimately form CH3OH. The hydride and proton transfers for the first of these reduction steps, the homogeneous reduction of CO2, are sequential in nature (in which the formate to formic acid protonation can be assisted by a negatively biased cathode). In contrast, these transfers are coupled in each of the two subsequent homogeneous hydride and proton transfer steps to reduce formic acid and formaldehyde.

Roles of the Lewis acid and base in the chemical reduction of CO2 catalyzed by frustrated Lewis pairs.

We employ quantum chemical calculations to discover how frustrated Lewis pairs (FLP) catalyze the reduction of CO2 by ammonia borane (AB); specifically, we examine how the Lewis acid (LA) and Lewis base (LB) of an FLP activate CO2 for reduction. We find that the LA (trichloroaluminum, AlCl3) alone catalyzes hydride transfer (HT) to CO2 while the LB (trimesitylenephosphine, PMes3) actually hinders HT; inclusion of the LB increases the HT barrier by ∼8 kcal/mol relative to the reaction catalyzed by LAs only. The LB hinders HT by donating its lone pair to the LUMO of CO2, increasing the electron density on the C atom and thus lowering its hydride affinity. Although the LB hinders HT, it nonetheless plays a crucial role by stabilizing the active FLP·CO2 complex relative to the LA dimer, free CO2, and free LB. This greatly increases the concentration of the reactive complex in the form FLP·CO2 and thus increases the rate of reaction. We expect that the principles we describe will aid in understanding other catalytic CO2 reductions.

Bulk and surface tunneling hydrogen defects in alumina.

We perform ab initio calculations of hydrogen-based tunneling defects in alumina to identify deleterious two-level systems (TLS) in superconducting qubits. The defects analyzed include bulk hydrogenated Al vacancies, bulk hydrogen interstitial defects, and a surface OH rotor. The formation energies of the defects are first computed for an Al- and O-rich environment to give the likelihood of defect occurrence during growth. The potential energy surfaces are then computed and the corresponding dipole moments are evaluated to determine the coupling of the defects to an electric field. Finally, the tunneling energy is computed for the hydrogen defect and the analogous deuterium defect, providing an estimate of the TLS energy and the corresponding frequency for photon absorption. We predict that hydrogenated cation vacancy defects will form a significant density of GHz frequency TLSs in alumina.

Mechanism of homogeneous reduction of CO2 by pyridine: proton relay in aqueous solvent and aromatic stabilization.

We employ quantum chemical calculations to investigate the mechanism of homogeneous CO(2) reduction by pyridine (Py) in the Py/p-GaP system. We find that CO(2) reduction by Py commences with PyCOOH(0) formation where: (a) protonated Py (PyH(+)) is reduced to PyH(0), (b) PyH(0) then reduces CO(2) by one electron transfer (ET) via nucleophilic attack by its N lone pair on the C of CO(2), and finally (c) proton transfer (PT) from PyH(0) to CO(2) produces PyCOOH(0). The predicted enthalpic barrier for this proton-coupled ET (PCET) reaction is 45.7 kcal/mol for direct PT from PyH(0) to CO(2). However, when PT is mediated by one to three water molecules acting as a proton relay, the barrier decreases to 29.5, 20.4, and 18.5 kcal/mol, respectively. The water proton relay reduces strain in the transition state (TS) and facilitates more complete ET. For PT mediated by a three water molecule proton relay, adding water molecules to explicitly solvate the core reaction system reduces the barrier to 13.6-16.5 kcal/mol, depending on the number and configuration of the solvating waters. This agrees with the experimentally determined barrier of 16.5 ± 2.4 kcal/mol. We calculate a pK(a) for PyH(0) of 31 indicating that PT preceding ET is highly unfavorable. Moreover, we demonstrate that ET precedes PT in PyCOOH(0) formation, confirming PyH(0)'s pK(a) as irrelevant for predicting PT from PyH(0) to CO(2). Furthermore, we calculate adiabatic electron affinities in aqueous solvent for CO(2), Py, and Py·CO(2) of 47.4, 37.9, and 66.3 kcal/mol respectively, indicating that the anionic complex PyCOO(-) stabilizes the anionic radicals CO(2)(-) and Py(-) to facilitate low barrier ET. As the reduction of CO(2) proceeds through ET and then PT, the pyridine ring becomes aromatic, and thus Py catalyzes CO(2) reduction by stabilizing the PCET TS and the PyCOOH(0) product through aromatic resonance stabilization. Our results suggest that Py catalyzes the homogeneous reductions of formic acid and formaldehyde en route to formation of CH(3)OH through a series of one-electron reductions analogous to the PCET reduction of CO(2) examined here, where the electrode only acts to reduce PyH(+) to PyH(0).