1
|
Xu B, Park J, Zhang D, De Santiago HA, Li W, Liu X, Luo J, Lany S, Qi Y. Local Ordering, Distortion, and Redox Activity in (La 0.75Sr 0.25)(Mn 0.25Fe 0.25Co 0.25Al 0.25)O 3 Investigated by a Computational Workflow for Compositionally Complex Perovskite Oxides. CHEMISTRY OF MATERIALS : A PUBLICATION OF THE AMERICAN CHEMICAL SOCIETY 2024; 36:4990-5001. [PMID: 38828189 PMCID: PMC11138530 DOI: 10.1021/acs.chemmater.3c03038] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/30/2023] [Revised: 04/28/2024] [Accepted: 04/30/2024] [Indexed: 06/05/2024]
Abstract
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.
Collapse
Affiliation(s)
- Boyuan Xu
- Department
of Physics, Brown University, Providence, Rhode Island 02912, United States
| | - Jiyun Park
- School
of Engineering, Brown University, Providence, Rhode Island 02912, United States
| | - Dawei Zhang
- Program
in Materials Science and Engineering, University
of California San Diego, La Jolla, California 92093, United States
| | - Héctor A De Santiago
- Department
of Mechanical, Materials and Aerospace Engineering, Benjamin M. Statler
College of Engineering and Mineral Resources, West Virginia University, Morgantown, West Virginia 26506, United States
| | - Wei Li
- Department
of Mechanical, Materials and Aerospace Engineering, Benjamin M. Statler
College of Engineering and Mineral Resources, West Virginia University, Morgantown, West Virginia 26506, United States
| | - Xingbo Liu
- Department
of Mechanical, Materials and Aerospace Engineering, Benjamin M. Statler
College of Engineering and Mineral Resources, West Virginia University, Morgantown, West Virginia 26506, United States
| | - Jian Luo
- Program
in Materials Science and Engineering, University
of California San Diego, La Jolla, California 92093, United States
- Department
of NanoEngineering, University of California
San Diego, La Jolla, California 92093, United States
| | - Stephan Lany
- Materials
Science Center, National Renewable Energy
Laboratory, Golden, Colorado 80401, United States
| | - Yue Qi
- School
of Engineering, Brown University, Providence, Rhode Island 02912, United States
| |
Collapse
|
2
|
Lany S. Chemical Potential Analysis as an Alternative to the van't Hoff Method: Hypothetical Limits of Solar Thermochemical Hydrogen. J Am Chem Soc 2024; 146:14114-14127. [PMID: 38739418 PMCID: PMC11117408 DOI: 10.1021/jacs.4c02688] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/22/2024] [Revised: 04/21/2024] [Accepted: 04/25/2024] [Indexed: 05/14/2024]
Abstract
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.
Collapse
Affiliation(s)
- Stephan Lany
- National Renewable Energy
Laboratory, Golden, Colorado 80401, United States
| |
Collapse
|
3
|
Streibel V, Schönecker JL, Wagner LI, Sirotti E, Munnik F, Kuhl M, Jiang CM, Eichhorn J, Santra S, Sharp ID. Zirconium Oxynitride Thin Films for Photoelectrochemical Water Splitting. ACS APPLIED ENERGY MATERIALS 2024; 7:4004-4015. [PMID: 38756865 PMCID: PMC11094725 DOI: 10.1021/acsaem.4c00303] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 02/04/2024] [Revised: 03/19/2024] [Accepted: 04/18/2024] [Indexed: 05/18/2024]
Abstract
Transition metal oxynitrides are a promising class of functional materials for photoelectrochemical (PEC) applications. Although these compounds are most commonly synthesized via ammonolysis of oxide precursors, such synthetic routes often lead to poorly controlled oxygen-to-nitrogen anion ratios, and the harsh nitridation conditions are incompatible with many substrates, including transparent conductive oxides. Here, we report direct reactive sputter deposition of a family of zirconium oxynitride thin films and the comprehensive characterization of their tunable structural, optical, and functional PEC properties. Systematic increases of the oxygen content in the reactive sputter gas mixture enable access to different crystalline structures within the zirconium oxynitride family. Increasing oxygen contents lead to a transition from metallic to semiconducting to insulating phases. In particular, crystalline Zr2ON2-like films have band gaps in the UV-visible range and are n-type semiconductors. These properties, together with a valence band maximum position located favorably relative to the water oxidation potential, make them viable photoanode candidates. Using chopped linear sweep voltammetry, we indeed confirm that our Zr2ON2 films are PEC-active for the oxygen evolution reaction in alkaline electrolytes. We further show that high-vacuum annealing boosts their PEC performance characteristics. Although the observed photocurrents are low compared to state-of-the-art photoanodes, these dense and planar thin films can offer a valuable platform for studying oxynitride photoelectrodes, as well as for future nanostructuring, band gap engineering, and defect engineering efforts.
Collapse
Affiliation(s)
- Verena Streibel
- Walter
Schottky Institute, Technical University
of Munich, Garching 85748, Germany
- Physics
Department, TUM School of Natural Sciences, Technical University of Munich, Garching 85748, Germany
| | - Johanna L. Schönecker
- Walter
Schottky Institute, Technical University
of Munich, Garching 85748, Germany
- Physics
Department, TUM School of Natural Sciences, Technical University of Munich, Garching 85748, Germany
| | - Laura I. Wagner
- Walter
Schottky Institute, Technical University
of Munich, Garching 85748, Germany
- Physics
Department, TUM School of Natural Sciences, Technical University of Munich, Garching 85748, Germany
| | - Elise Sirotti
- Walter
Schottky Institute, Technical University
of Munich, Garching 85748, Germany
- Physics
Department, TUM School of Natural Sciences, Technical University of Munich, Garching 85748, Germany
| | - Frans Munnik
- Institute
of Ion Beam Physics and Materials Research, Helmholtz-Zentrum Dresden-Rossendorf (HZDR), Dresden 01328, Germany
| | - Matthias Kuhl
- Walter
Schottky Institute, Technical University
of Munich, Garching 85748, Germany
- Physics
Department, TUM School of Natural Sciences, Technical University of Munich, Garching 85748, Germany
| | - Chang-Ming Jiang
- Walter
Schottky Institute, Technical University
of Munich, Garching 85748, Germany
- Physics
Department, TUM School of Natural Sciences, Technical University of Munich, Garching 85748, Germany
| | - Johanna Eichhorn
- Walter
Schottky Institute, Technical University
of Munich, Garching 85748, Germany
- Physics
Department, TUM School of Natural Sciences, Technical University of Munich, Garching 85748, Germany
| | - Saswati Santra
- Walter
Schottky Institute, Technical University
of Munich, Garching 85748, Germany
- Physics
Department, TUM School of Natural Sciences, Technical University of Munich, Garching 85748, Germany
| | - Ian D. Sharp
- Walter
Schottky Institute, Technical University
of Munich, Garching 85748, Germany
- Physics
Department, TUM School of Natural Sciences, Technical University of Munich, Garching 85748, Germany
| |
Collapse
|
4
|
Kistanov AA. Atomic insights into the interaction of N 2, CO 2, NH 3, NO, and NO 2 gas molecules with Zn 2(V, Nb, Ta)N 3 ternary nitride monolayers. Phys Chem Chem Phys 2024; 26:13719-13730. [PMID: 38669029 DOI: 10.1039/d4cp01225a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/09/2024]
Abstract
The search for promising carrier blocking layer materials with high stability, including resistance to surface inhibition by environmental molecules that cause a drop in carrier mobility, is critical for the production of tandem solar cells. Based on density functional theory calculations, the reaction of atmospheric gases, including N2, CO2, NH3, NO, and NO2, with three promising Zn2(V, Nb, Ta)N3 monolayers is discovered. The results suggest the chemical adsorption of NH3 and physical adsorption of NO and NO2. In addition, the Zn2(V, Nb, Ta)N3 monolayers are characterized by a weak bonding with N2 and CO2. Charge redistribution is found at the interface between the monolayers and NH3, NO and NO2 molecules, leading to the formation of a local surface dipole that affects the functionality of the Zn2(V, Nb, Ta)N3 monolayers. The Zn2VN3 monolayer is less reactive with atmospheric gases and thus is the most promising for application in tandem solar cells. Notably, the revealed nontrivial behavior of the Zn2(V, Nb, Ta)N3 monolayers towards N-containing gases makes them promising for application in gas sensing. Specifically, the Zn2TaN3 monolayer is the most promising for application in molecular sensing due to its high reversibility and distinguished interaction with NH3, NO, and NO2 gases.
Collapse
Affiliation(s)
- Andrey A Kistanov
- The Laboratory of Metals and Alloys Under Extreme Impacts, Ufa University of Science and Technology, Ufa 450076, Russia.
| |
Collapse
|
5
|
Witman MD, Goyal A, Ogitsu T, McDaniel AH, Lany S. Defect graph neural networks for materials discovery in high-temperature clean-energy applications. NATURE COMPUTATIONAL SCIENCE 2023; 3:675-686. [PMID: 38177319 DOI: 10.1038/s43588-023-00495-2] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/15/2022] [Accepted: 07/05/2023] [Indexed: 01/06/2024]
Abstract
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.
Collapse
Affiliation(s)
| | - Anuj Goyal
- National Renewable Energy Laboratory, Golden, CO, USA
- Indian Institute of Technology Hyderabad, Kandi, Telangana, India
| | - Tadashi Ogitsu
- Lawrence Livermore National Laboratory, Livermore, CA, USA
| | | | - Stephan Lany
- National Renewable Energy Laboratory, Golden, CO, USA.
| |
Collapse
|
6
|
Affiliation(s)
- Josiah Roberts
- Department of Chemistry, State University of New York at Buffalo, Buffalo, New York 14260-3000, USA
| | - Eva Zurek
- Department of Chemistry, State University of New York at Buffalo, Buffalo, New York 14260-3000, USA
| |
Collapse
|
7
|
Strange NA, Park JE, Goyal A, Bell RT, Trindell JA, Sugar JD, Stone KH, Coker EN, Lany S, Shulda S, Ginley DS. Formation of 6H-Ba 3Ce 0.75Mn 2.25O 9 during Thermochemical Reduction of 12R-Ba 4CeMn 3O 12: Identification of a Polytype in the Ba(Ce,Mn)O 3 Family. Inorg Chem 2022; 61:6128-6137. [PMID: 35404603 DOI: 10.1021/acs.inorgchem.2c00282] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
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.
Collapse
Affiliation(s)
- Nicholas A Strange
- SLAC National Accelerator Laboratory, Menlo Park, California 94025, United States
| | - James Eujin Park
- Sandia National Laboratories, Albuquerque, New Mexico 87185, United States
| | - Anuj Goyal
- National Renewable Energy Laboratory, Golden, Colorado 80401, United States
| | - Robert T Bell
- National Renewable Energy Laboratory, Golden, Colorado 80401, United States
| | - Jamie A Trindell
- Sandia National Laboratories, Livermore, California 94550, United States
| | - Joshua D Sugar
- Sandia National Laboratories, Livermore, California 94550, United States
| | - Kevin H Stone
- SLAC National Accelerator Laboratory, Menlo Park, California 94025, United States
| | - Eric N Coker
- Sandia National Laboratories, Albuquerque, New Mexico 87185, United States
| | - Stephan Lany
- National Renewable Energy Laboratory, Golden, Colorado 80401, United States
| | - Sarah Shulda
- National Renewable Energy Laboratory, Golden, Colorado 80401, United States
| | - David S Ginley
- National Renewable Energy Laboratory, Golden, Colorado 80401, United States
| |
Collapse
|
8
|
Sharan A, Singh N. Intrinsic Valley Polarization in Computationally Discovered Two‐Dimensional Ferrovalley Materials: LaI
2
and PrI
2
Monolayers. ADVANCED THEORY AND SIMULATIONS 2022. [DOI: 10.1002/adts.202100476] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Affiliation(s)
- Abhishek Sharan
- Department of Physics Khalifa University of Science and Technology Abu Dhabi 127788 UAE
| | - Nirpendra Singh
- Department of Physics Khalifa University of Science and Technology Abu Dhabi 127788 UAE
| |
Collapse
|