1
|
Kanhaiya K, Nathanson M, In 't Veld PJ, Zhu C, Nikiforov I, Tadmor EB, Choi YK, Im W, Mishra RK, Heinz H. Accurate Force Fields for Atomistic Simulations of Oxides, Hydroxides, and Organic Hybrid Materials up to the Micrometer Scale. J Chem Theory Comput 2023; 19:8293-8322. [PMID: 37962992 DOI: 10.1021/acs.jctc.3c00750] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2023]
Abstract
The simulation of metals, oxides, and hydroxides can accelerate the design of therapeutics, alloys, catalysts, cement-based materials, ceramics, bioinspired composites, and glasses. Here we introduce the INTERFACE force field (IFF) and surface models for α-Al2O3, α-Cr2O3, α-Fe2O3, NiO, CaO, MgO, β-Ca(OH)2, β-Mg(OH)2, and β-Ni(OH)2. The force field parameters are nonbonded, including atomic charges for Coulomb interactions, Lennard-Jones (LJ) potentials for van der Waals interactions with 12-6 and 9-6 options, and harmonic bond stretching for hydroxide ions. The models outperform DFT calculations and earlier atomistic models (Pedone, ReaxFF, UFF, CLAYFF) up to 2 orders of magnitude in reliability, compatibility, and interpretability due to a quantitative representation of chemical bonding consistent with other compounds across the periodic table and curated experimental data for validation. The IFF models exhibit average deviations of 0.2% in lattice parameters, <10% in surface energies (to the extent known), and 6% in bulk moduli relative to experiments. The parameters and models can be used with existing parameters for solvents, inorganic compounds, organic compounds, biomolecules, and polymers in IFF, CHARMM, CVFF, AMBER, OPLS-AA, PCFF, and COMPASS, to simulate bulk oxides, hydroxides, electrolyte interfaces, and multiphase, biological, and organic hybrid materials at length scales from atoms to micrometers. The nonbonded character of the models also enables the analysis of mixed oxides, glasses, and certain chemical reactions, and well-performing nonbonded models for silica phases, SiO2, are introduced. Automated model building is available in the CHARMM-GUI Nanomaterial Modeler. We illustrate applications of the models to predict the structure of mixed oxides, and energy barriers of ion migration, as well as binding energies of water and organic molecules in outstanding agreement with experimental data and calculations at the CCSD(T) level. Examples of model building for hydrated, pH-sensitive oxide surfaces to simulate solid-electrolyte interfaces are discussed.
Collapse
Affiliation(s)
- Krishan Kanhaiya
- Department of Chemical and Biological Engineering, University of Colorado at Boulder, Boulder, Colorado 80309, United States
| | - Michael Nathanson
- Department of Chemical and Biological Engineering, University of Colorado at Boulder, Boulder, Colorado 80309, United States
| | - Pieter J In 't Veld
- BASF SE, Molecular Modeling & Drug Discovery, Carl Bosch Str. 38, 67056 Ludwigshafen, Germany
| | - Cheng Zhu
- Department of Chemical and Biological Engineering, University of Colorado at Boulder, Boulder, Colorado 80309, United States
| | - Ilia Nikiforov
- Department of Aerospace Engineering and Mechanics, University of Minnesota, Minneapolis, Minnesota 55455, United States
| | - Ellad B Tadmor
- Department of Aerospace Engineering and Mechanics, University of Minnesota, Minneapolis, Minnesota 55455, United States
| | - Yeol Kyo Choi
- Department of Biological Sciences, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| | - Wonpil Im
- Department of Biological Sciences, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| | - Ratan K Mishra
- BASF SE, Molecular Modeling & Drug Discovery, Carl Bosch Str. 38, 67056 Ludwigshafen, Germany
| | - Hendrik Heinz
- Department of Chemical and Biological Engineering, University of Colorado at Boulder, Boulder, Colorado 80309, United States
| |
Collapse
|
3
|
Microstructure Evolution in a Solid Oxide Fuel Cell Stack Quantified with Interfacial Free Energy. ENERGIES 2021. [DOI: 10.3390/en14123476] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
Degradation of electrode microstructure is one of the key factors affecting long term performance of Solid Oxide Fuel Cell systems. Evolution of a multiphase system can be described quantitatively by the change in its interfacial energy. In this paper, we discuss free energy of a microstructure to showcase the anisotropy of its evolution during a long-term performance experiment involving an SOFC stack. Ginzburg Landau type functional is used to compute the free energy, using diffuse phase distributions based on Focused Ion Beam Scanning Electron Microscopy images of samples taken from nine different sites within the stack. It is shown that the rate of microstructure evolution differs depending on the position within the stack, similar to phase anisotropy. However, the computed spatial relation does not correlate with the observed distribution of temperature.
Collapse
|
4
|
Valadez Huerta G, Raabe G. Genetic Parameterization of Interfacial Force Fields Based on Classical Bulk Force Fields and Ab Initio Data: Application to the Methanol-ZnO Interfaces. J Chem Inf Model 2020; 60:6033-6043. [DOI: 10.1021/acs.jcim.0c01093] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Gerardo Valadez Huerta
- Institut für Thermodynamik, Technische Universität Braunschweig, Hans-Sommer-Straße 5, D-38106 Braunschweig, Germany
| | - Gabriele Raabe
- Institut für Thermodynamik, Technische Universität Braunschweig, Hans-Sommer-Straße 5, D-38106 Braunschweig, Germany
| |
Collapse
|
5
|
Kalra P, Garg RK, Kumar A. Parametric sensitivity analysis for a natural gas fueled high temperature tubular solid oxide fuel cell. Heliyon 2020; 6:e04450. [PMID: 32760822 PMCID: PMC7393395 DOI: 10.1016/j.heliyon.2020.e04450] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2019] [Revised: 03/06/2020] [Accepted: 07/10/2020] [Indexed: 11/26/2022] Open
Abstract
For more than the last two decades, there has been research going on to develop advanced energy technologies involving minimum environmental pollution, to replace the conventional fossil energy systems. Solid oxide fuel cells (SOFCs) are one of the most promising, eco-friendly and efficient means for the generation of electricity to meet the future energy requirements. This research work focuses on the parametric sensitivity analysis for natural gas fueled high temperature tubular SOFC. Firstly, for the tubular SOFC, a one-dimensional radially symmetrical model has been developed and solved using the finite-difference method. Then, the effect of the variation of important operational and design parameters on its performance has been analyzed. The parameters typically include composition, inlet fuel temperature, pressure, length of SOFC tube and thicknesses of its components. The composition is expressed as steam to methane ratio and it has been observed that the voltage and power density developed by the SOFC diminishes as the ratio increases. Further, a change in the inlet fuel pressure of the tubular SOFC has no pronounced influence on the voltage and power density developed. On the other hand, with an increase in the inlet fuel temperature, a small improvement in these performance characteristics is exhibited. The axial length of the tubular SOFC does conspicuously influence its performance characteristics but solely at current densities greater than 4000A/m2. An increase in the thickness of its components results in a reduction in its voltage and power density developed. The largest decline in these performance characteristics with the increase in thickness is observed for electrolyte followed by cathode and anode respectively.
Collapse
Affiliation(s)
- Pankaj Kalra
- Department of Chemical Engineering, SBSSTC, Ferozepur, 152004, Punjab, India
| | - Rajeev Kumar Garg
- Department of Chemical Engineering, SBSSTC, Ferozepur, 152004, Punjab, India
| | - Ajay Kumar
- Department of Applied Science & Humanities, SBSSTC, Ferozepur, 152004, Punjab India
| |
Collapse
|