Chemical Potential Analysis as an Alternative to the van't Hoff Method: Hypothetical Limits of Solar Thermochemical Hydrogen.

The van't Hoff method is a standard approach for determining reaction enthalpies and entropies, e.g., in the thermochemical reduction of oxides, which is an important process for solar thermochemical fuels and numerous other applications. However, by analyzing the oxygen partial pressure pO2, e.g., as measured by thermogravimetric analysis (TGA), this method convolutes the properties of the probe gas with the solid-state properties of the examined oxides, which define their suitability for specific applications. The "chemical potential method" is here proposed as an alternative. Using the oxygen chemical potential ΔµO instead of pO2 for the analysis, this method does not only decouple gas-phase and solid-state contributions but also affords a simple and transparent approach to extracting the temperature dependence of the reduction enthalpy and entropy, which carries important information about the defect mechanism. For demonstration of the approach, this work considers three model systems; (1) a generic oxide with noninteracting, charge-neutral oxygen vacancy defects, (2) Sr0.86Ce0.14MnO3(1-δ) alloys with interacting vacancies, and (3) a model for charged vacancy formation in CeO2, which reproduces the extensive experimental TGA data available in the literature. The reduction behavior of these model systems obtained from the chemical potential method is correlated with simulated results for the thermochemical water splitting cycle, highlighting the exceptional behavior of CeO2, which originates from defect ionization. The theoretical performance limits for solar thermochemical hydrogen within the charged defect mechanism are assessed by considering hypothetical materials described by a variation of the CeO2 model parameters within a plausible range.

"Mn3AlN" is Really Mn4N.

We investigate the synthesis of antiperovskite "Mn3AlN" using the published synthesis procedure, as well as several new reaction pathways. In each case, only a combination of antiperovskite Mn4N and Mn5Al8 or precursors is obtained. The identity of the obtained antiperovskite phase is unambiguously determined to be Mn4N via synchrotron powder X-ray diffraction (SPXRD), X-ray absorption spectroscopy (XAS), and magnetometry. The experimental results are further supported by thermochemical calculations informed by density functional theory (DFT), which find Mn3AlN to be metastable versus decomposition into Mn and AlN. The DFT-based calculations also predict an antiferromagnetic ground state for Mn3AlN. This directly contradicts the previously reported ferromagnetic behavior of "Mn3AlN". Instead, the observed magnetic behavior is consistent with ferrimagnetic Mn4N. We examine the data in the original publication and conclude that the compound reported to be Mn3AlN is in fact Mn4N.

Formation of 6H-Ba3Ce0.75Mn2.25O9 during Thermochemical Reduction of 12R-Ba4CeMn3O12: Identification of a Polytype in the Ba(Ce,Mn)O3 Family.

The resurgence of interest in a hydrogen economy and the development of hydrogen-related technologies has initiated numerous research and development efforts aimed at making the generation, storage, and transportation of hydrogen more efficient and affordable. Solar thermochemical hydrogen production (STCH) is a process that potentially exhibits numerous benefits such as high reaction efficiencies, tunable thermodynamics, and continued performance over extended cycling. Although CeO2 has been the de facto standard STCH material for many years, more recently 12R-Ba4CeMn3O12 (BCM) has demonstrated enhanced hydrogen production at intermediate H2/H2O conditions compared to CeO2, making it a contender for large-scale hydrogen production. However, the thermo-reduction stability of 12R-BCM dictates the oxygen partial pressure (pO2) and temperature conditions optimal for cycling. In this study, we identify the formation of a 6H-BCM polytype at high temperature and reducing conditions, experimentally and computationally, as a mechanism and pathway for 12R-BCM decomposition. 12R-BCM was synthesized with high purity and then controllably reduced using thermogravimetric analysis (TGA). Synchrotron X-ray diffraction (XRD) data is used to identify the formation of a 6H-Ba3Ce0.75Mn2.25O9 (6H-BCM) polytype that is formed at 1350 °C under strongly reducing pO2. Density functional theory (DFT) total energy and defect calculations show a window of thermodynamic stability for the 6H-polytype consistent with the XRD results. These data provide the first evidence of the 6H-BCM polytype and could provide a mechanistic explanation for the superior water-splitting behaviors of 12R-BCM.

Ternary nitride semiconductors in the rocksalt crystal structure.

Inorganic nitrides with wurtzite crystal structures are well-known semiconductors used in optical and electronic devices. In contrast, rocksalt-structured nitrides are known for their superconducting and refractory properties. Breaking this dichotomy, here we report ternary nitride semiconductors with rocksalt crystal structures, remarkable electronic properties, and the general chemical formula Mgx TM 1-xN (TM = Ti, Zr, Hf, Nb). Our experiments show that these materials form over a broad metal composition range, and that Mg-rich compositions are nondegenerate semiconductors with visible-range optical absorption onsets (1.8 to 2.1 eV) and up to 100 cm2 V-1⋅s-1 electron mobility for MgZrN2 grown on MgO substrates. Complementary ab initio calculations reveal that these materials have disorder-tunable optical absorption, large dielectric constants, and electronic bandgaps that are relatively insensitive to disorder. These ternary Mgx TM 1-xN semiconductors are also structurally compatible both with binary TMN superconductors and main-group nitride semiconductors along certain crystallographic orientations. Overall, these results highlight Mgx TM 1-xN as a class of materials combining the semiconducting properties of main-group wurtzite nitrides and rocksalt structure of superconducting transition-metal nitrides.

Metal chalcogenides for neuromorphic computing: emerging materials and mechanisms.

The approaching end of Moore's law scaling has significantly accelerated multiple fields of research including neuromorphic-, quantum-, and photonic computing, each of which possesses unique benefits unobtained through conventional binary computers. One of the most compelling arguments for neuromorphic computing systems is power consumption, noting that computations made in the human brain are approximately 106times more efficient than conventional CMOS logic. This review article focuses on the materials science and physical mechanisms found in metal chalcogenides that are currently being explored for use in neuromorphic applications. We begin by reviewing the key biological signal generation and transduction mechanisms within neuronal components of mammalian brains and subsequently compare with observed experimental measurements in chalcogenides. With robustness and energy efficiency in mind, we will focus on short-range mechanisms such as structural phase changes and correlated electron systems that can be driven by low-energy stimuli, such as temperature or electric field. We aim to highlight fundamental materials research and existing gaps that need to be overcome to enable further integration or advancement of metal chalcogenides for neuromorphic systems.

Computational discovery of stable and metastable ternary oxynitrides.

Materials design from first principles enables exploration of uncharted chemical spaces. Extensive computational searches have been performed for mixed-cation ternary compounds, but mixed-anion systems are gaining increased interest as well. Central to computational discovery is the crystal structure prediction, where the trade-off between reliance on prototype structures and size limitations of unconstrained sampling has to be navigated. We approach this challenge by letting two complementary structure sampling approaches compete. We use the kinetically limited minimization approach for high-throughput unconstrained crystal structure prediction in smaller cells up to 21 atoms. On the other hand, ternary-and, more generally, multinary-systems often assume structures formed by atomic ordering on a lattice derived from a binary parent structure. Thus, we additionally sample atomic configurations on prototype lattices with cells up to 56 atoms. Using this approach, we searched 65 different charge-balanced oxide-nitride stoichiometries, including six known systems as the control sample. The convex hull analysis is performed both for the thermodynamic limit and for the case of synthesis with activated nitrogen sources. We identified 34 phases that are either on the convex hull or within a viable energy window for potentially metastable phases. We further performed structure sampling for "missing" binary nitrides whose energies are needed for the convex hull analysis. Among these, we discovered metastable Ce3N4 as a nitride analog of the tetravalent cerium oxide, which becomes stable under slightly activated nitrogen condition ΔµN > +0.07 eV. Given the outsize role of CeO2 in research and application, Ce3N4 is a potentially important discovery.

Combinatorial Synthesis of Magnesium Tin Nitride Semiconductors.

Nitride materials feature strong chemical bonding character that leads to unique crystal structures, but many ternary nitride chemical spaces remain experimentally unexplored. The search for previously undiscovered ternary nitrides is also an opportunity to explore unique materials properties, such as transitions between cation-ordered and -disordered structures, as well as to identify candidate materials for optoelectronic applications. Here, we present a comprehensive experimental study of MgSnN2, an emerging II-IV-N2 compound, for the first time mapping phase composition and crystal structure, and examining its optoelectronic properties computationally and experimentally. We demonstrate combinatorial cosputtering of cation-disordered, wurtzite-type MgSnN2 across a range of cation compositions and temperatures, as well as the unexpected formation of a secondary, rocksalt-type phase of MgSnN2 at Mg-rich compositions and low temperatures. A computational structure search shows that the rocksalt-type phase is substantially metastable (>70 meV/atom) compared to the wurtzite-type ground state. Spectroscopic ellipsometry reveals optical absorption onsets around 2 eV, consistent with band gap tuning via cation disorder. Finally, we demonstrate epitaxial growth of a mixed wurtzite-rocksalt MgSnN2 on GaN, highlighting an opportunity for polymorphic control via epitaxy. Collectively, these findings lay the groundwork for further exploration of MgSnN2 as a model ternary nitride, with controlled polymorphism, and for device applications, enabled by control of optoelectronic properties via cation ordering.

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.

Redox-Mediated Stabilization in Zinc Molybdenum Nitrides.

We report on the theoretical prediction and experimental realization of new ternary zinc molybdenum nitride compounds. We used theory to identify previously unknown ternary compounds in the Zn-Mo-N systems, Zn3MoN4 and ZnMoN2, and to analyze their bonding environment. Experiments show that Zn-Mo-N alloys can form in broad composition range from Zn3MoN4 to ZnMoN2 in the wurtzite-derived structure, accommodating very large off-stoichiometry. Interestingly, the measured wurtzite-derived structure of the alloys is metastable for the ZnMoN2 stoichiometry, in contrast to the Zn3MoN4 stoichiometry, where ordered wurtzite is predicted to be the ground state. The formation of Zn3MoN4-ZnMoN2 alloy with wurtzite-derived crystal structure is enabled by the concomitant ability of Mo to change oxidation state from +VI in Zn3MoN4 to +IV in ZnMoN2, and the capability of Zn to contribute to the bonding states of both compounds, an effect that we define as "redox-mediated stabilization". The stabilization of Mo in both the +VI and +IV oxidation states is due to the intermediate electronegativity of Zn, which enables significant polar covalent bonding in both Zn3MoN4 and ZnMoN2 compounds. The smooth change in the Mo oxidation state between Zn3MoN4 and ZnMoN2 stoichiometries leads to a continuous change in optoelectronic properties-from resistive and semitransparent Zn3MoN4 to conductive and absorptive ZnMoN2. The reported redox-mediated stabilization in zinc molybdenum nitrides suggests there might be many undiscovered ternary compounds with one metal having an intermediate electronegativity, enabling significant covalent bonding, and another metal capable of accommodating multiple oxidation states, enabling stoichiometric flexibility.

Communication: The electronic entropy of charged defect formation and its impact on thermochemical redox cycles.

The ideal material for solar thermochemical water splitting, which has yet to be discovered, must satisfy stringent conditions for the free energy of reduction, including, in particular, a sufficiently large positive contribution from the solid-state entropy. By inverting the commonly used relationship between defect formation energy and defect concentration, it is shown here that charged defect formation causes a large electronic entropy contribution manifesting itself as the temperature dependence of the Fermi level. This result is a general feature of charged defect formation and motivates new materials design principles for solar thermochemical hydrogen production.

Synthesis of a mixed-valent tin nitride and considerations of its possible crystal structures.

Recent advances in theoretical structure prediction methods and high-throughput computational techniques are revolutionizing experimental discovery of the thermodynamically stable inorganic materials. Metastable materials represent a new frontier for these studies, since even simple binary non-ground state compounds of common elements may be awaiting discovery. However, there are significant research challenges related to non-equilibrium thin film synthesis and crystal structure predictions, such as small strained crystals in the experimental samples and energy minimization based theoretical algorithms. Here, we report on experimental synthesis and characterization, as well as theoretical first-principles calculations of a previously unreported mixed-valent binary tin nitride. Thin film experiments indicate that this novel material is N-deficient SnN with tin in the mixed ii/iv valence state and a small low-symmetry unit cell. Theoretical calculations suggest that the most likely crystal structure has the space group 2 (SG2) related to the distorted delafossite (SG166), which is nearly 0.1 eV/atom above the ground state SnN polymorph. This observation is rationalized by the structural similarity of the SnN distorted delafossite to the chemically related Sn3N4 spinel compound, which provides a fresh scientific insight into the reasons for growth of polymorphs of metastable materials. In addition to reporting on the discovery of the simple binary SnN compound, this paper illustrates a possible way of combining a wide range of advanced characterization techniques with the first-principle property calculation methods, to elucidate the most likely crystal structure of the previously unreported metastable materials.

Effect of extended strain fields on point defect phonon scattering in thermoelectric materials.

The design of thermoelectric materials often involves the integration of point defects (alloying) as a route to reduce the lattice thermal conductivity. Classically, the point defect scattering strength follows from simple considerations such as mass contrast and the presence of induced strain fields (e.g. radius contrast, coordination changes). While the mass contrast can be easily calculated, the associated strain fields induced by defect chemistry are not readily predicted and are poorly understood. In this work, we use classical and first principles calculations to provide insight into the strain field component of phonon scattering from isoelectronic point defects. Our results also integrate experimental measurements on bulk samples of SnSe and associated alloys with S, Te, Ge, Sr and Ba. These efforts highlight that the strength and extent of the resulting strain field depends strongly on defect chemistry. Strain fields can have a profound impact on the local structure. For example, in alloys containing Ba, the strain fields have significant spatial extent (1 nm in diameter) and produce large shifts in the atomic equilibrium positions (up to 0.5 Å). Such chemical complexity suggests that computational assessment of point defects for thermal conductivity depression should be hindered. However, in this work, we present and verify several computational descriptors that correlate well with the experimentally measured strain fields. Furthermore, these descriptors are conceptually transparent and computationally inexpensive, allowing computation to provide a pivotal role in the screening of effective alloys. The further development of point defect engineering could complement or replace nanostructuring when optimizing the thermal conductivity, offering the benefits of thermodynamic stability, and providing more clearly defined defect chemistry.

Assessing capability of semiconductors to split water using ionization potentials and electron affinities only.

We show in this article that the position of semiconductor band edges relative to the water reduction and oxidation levels can be reliably predicted from the ionization potentials (IP) and electron affinities (AE) only. Using a set of 17 materials, including transition metal compounds, we show that accurate surface dependent IPs and EAs of semiconductors can be computed by combining density functional theory and many-body GW calculations. From the extensive comparison of calculated IPs and EAs with available experimental data, both from photoemission and electrochemical measurements, we show that it is possible to sort candidate materials solely from IPs and EAs thereby eliminating explicit treatment of semiconductor/water interfaces. We find that at pH values corresponding to the point of zero charge there is on average a 0.5 eV shift of IPs and EAs closer to the vacuum due to the dipoles formed at material/water interfaces.

Charge self-regulation upon changing the oxidation state of transition metals in insulators.

Transition-metal atoms embedded in an ionic or semiconducting crystal can exist in various oxidation states that have distinct signatures in X-ray photoemission spectroscopy and 'ionic radii' which vary with the oxidation state of the atom. These oxidation states are often tacitly associated with a physical ionization of the transition-metal atoms--that is, a literal transfer of charge to or from the atoms. Physical models have been founded on this charge-transfer paradigm, but first-principles quantum mechanical calculations show only negligible changes in the local transition-metal charge as the oxidation state is altered. Here we explain this peculiar tendency of transition-metal atoms to maintain a constant local charge under external perturbations in terms of an inherent, homeostasis-like negative feedback. We show that signatures of oxidation states and multivalence--such as X-ray photoemission core-level shifts, ionic radii and variations in local magnetization--that have often been interpreted as literal charge transfer are instead a consequence of the negative-feedback charge regulation.

Local Ordering, Distortion, and Redox Activity in (La0.75Sr0.25)(Mn0.25Fe0.25Co0.25Al0.25)O3 Investigated by a Computational Workflow for Compositionally Complex Perovskite Oxides.

Mixing multiple cations can result in a significant configurational entropy, offer a new compositional space with vast tunability, and introduce new computational challenges. For applications such as the two-step solar thermochemical hydrogen (STCH) generation techniques, we demonstrate that using density functional theory (DFT) combined with Metropolis Monte Carlo method (DFT-MC) can efficiently sample the possible cation configurations in compositionally complex perovskite oxide (CCPO) materials, with (La0.75Sr0.25)(Mn0.25Fe0.25Co0.25Al0.25)O3 as an example. In the presence of oxygen vacancies (VO), DFT-MC simulations reveal a significant increase of the local site preference of the cations (short-range ordering), compared to a more random mixing without VO. Co is found to be the redox-active element and the VO is the preferentially generated next to Co due to the stretched Co-O bonds. A clear definition of the vacancy formation energy (Evf) is proposed for CCPO in an ensemble of structures evolved in parallel from independent DFT-MC paths. By combining the distribution of Evf with VO interactions into a statistical model, the oxygen nonstoichiometry (δ), under the STCH thermal reduction and oxidation conditions, is predicted and compared with the experiments. Similar to the experiments, the predicted δ can be used to extract the enthalpy and entropy of reduction using the van't Hoff method, providing direct comparisons with the experimental results. This procedure provides a full predictive workflow for using DFT-MC to obtain possible local ordering or fully random structures, understand the redox activity of each element, and predict the thermodynamic properties of CCPOs, for computational screening and design of these CCPO materials at STCH conditions.

Theoretical prediction and experimental realization of new stable inorganic materials using the inverse design approach.

Discovery of new materials is important for all fields of chemistry. Yet, existing compilations of all known ternary inorganic solids still miss many possible combinations. Here, we present an example of accelerated discovery of the missing materials using the inverse design approach, which couples predictive first-principles theoretical calculations with combinatorial and traditional experimental synthesis and characterization. The compounds in focus belong to the equiatomic (1:1:1) ABX family of ternary materials with 18 valence electrons per formula unit. Of the 45 possible V-IX-IV compounds, 29 are missing. Theoretical screening of their thermodynamic stability revealed eight new stable 1:1:1 compounds, including TaCoSn. Experimental synthesis of TaCoSn, the first ternary in the Ta-Co-Sn system, confirmed its predicted zincblende-derived crystal structure. These results demonstrate how discovery of new materials can be accelerated by the combination of high-throughput theoretical and experimental methods. Despite being made of three metallic elements, TaCoSn is predicted and explained to be a semiconductor. The band gap of this material is difficult to measure experimentally, probably due to a high concentration of interstitial cobalt defects.

Defect graph neural networks for materials discovery in high-temperature clean-energy applications.

We present a graph neural network approach that fully automates the prediction of defect formation enthalpies for any crystallographic site from the ideal crystal structure, without the need to create defected atomic structure models as input. Here we used density functional theory reference data for vacancy defects in oxides, to train a defect graph neural network (dGNN) model that replaces the density functional theory supercell relaxations otherwise required for each symmetrically unique crystal site. Interfaced with thermodynamic calculations of reduction entropies and associated free energies, the dGNN model is applied to the screening of oxides in the Materials Project database, connecting the zero-kelvin defect enthalpies to high-temperature process conditions relevant for solar thermochemical hydrogen production and other energy applications. The dGNN approach is applicable to arbitrary structures with an accuracy limited principally by the amount and diversity of the training data, and it is generalizable to other defect types and advanced graph convolution architectures. It will help to tackle future materials discovery problems in clean energy and beyond.

Two-dimensional polaronic behavior in the binary oxides m-HfO2 and m-ZrO2.

We demonstrate that the three-dimensional (3D) binary monoclinic oxides HfO2 and ZrO2 exhibit quasi-2D polaron localization and conductivity, which results from a small difference in the coordination of two oxygen sublattices in these materials. The transition between a 2D large polaron into a zero-dimensional small polaron state requires overcoming a small energetic barrier. These results demonstrate how a small asymmetry in the lattice structure can determine the qualitative character of polaron localization and significantly broaden the realm of quasi-2D polaron systems.

Interplay between Composition, Electronic Structure, Disorder, and Doping due to Dual Sublattice Mixing in Nonequilibrium Synthesis of ZnSnN2 :O.

The opportunity for enhanced functional properties in semiconductor solid solutions has attracted vast scientific interest for a variety of novel applications. However, the functional versatility originating from the additional degrees of freedom due to atomic composition and ordering comes along with new challenges in characterization and modeling. Developing predictive synthesis-structure-property relationships is prerequisite for effective materials design strategies. Here, a first-principles based model for property prediction in such complex semiconductor materials is presented. This framework incorporates nonequilibrium synthesis, dopants and defects, and the change of the electronic structure with composition and short range order. This approach is applied to ZnSnN2 (ZTN) which has attracted recent interest for photovoltaics. The unintentional oxygen incorporation and its correlation with the cation stoichiometry leads to the formation of a solid solution with dual sublattice mixing. A nonmonotonic doping behavior as a function of the composition is uncovered. The degenerate doping of near-stoichiometric ZTN, which is detrimental for potential applications, can be lowered into the 1017 cm-3 range in highly off-stoichiometric material, in quantitative agreement with experiments.

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.