Abstract:Proton-conducting perovskites such as Y-doped BaZrO 3 (BYZ) are promising candidates as electrolytes for a proton ceramic fuel cell (PCFC) that might permit much lower temperatures (from 400 to 600 degrees C). However, these materials lead to relatively poor total conductivity ( approximately 10 (-4) S/cm) because of extremely high grain boundary resistance. In order to provide the basis for improving these materials, we developed the ReaxFF reactive force field to enable molecular dynamics (MD) simulations of proton diffusion in the bulk phase and across grain boundaries of BYZ. This allows us to elucidate the atomistic structural details underlying the origin of this poor grain boundary conductivity and how it is related to the orientation of the grains. The parameters in ReaxFF were based entirely on the results of quantum mechanics (QM) calculations for systems related to BYZ. We apply here the ReaxFF to describe the proton diffusion in crystalline BYZ and across grain boundaries in BYZ. The results are in excellent agreement with experiment, validating the use of ReaxFF for studying the transport properties of these membranes. Having atomistic structures for the grain boundaries from simulations that explain the overall effect of the grain boundaries on diffusion opens the door to in silico optimization of these materials. That is, we can now use theory and simulation to examine the effect of alloying on both the interfacial structures and on the overall diffusion. As an example, these calculations suggest that the reduced diffusion of protons across the grain boundary results from the increased average distances between oxygen atoms in the interface, which necessarily leads to larger barriers for proton hopping. Assuming that this is the critical issue in grain boundary diffusion, the performance of BYZ for multigranular systems might be improved using additives that would tend to precipitate to the grain boundary and which would tend to pull the oxygens atoms together. Possibilities might be to use a small amount of larger trivalent ions, such as La or Lu or of tetravalent ions such as Hf or Th. Since ReaxFF can also be used to describe the chemical processes on the anode and cathode and the migration of ions across the electrode-membrane interface, ReaxFF opens the door to the possibility of atomistic first principles predictions on models of a complete fuel cell.

ReaxFF reactive force field for molecular dynamics simulations of liquid Cu and Zr metals

Development of a ReaxFF reactive force field for aqueous chloride and copper chloride

Development and Validation of a General-Purpose ReaxFF Reactive Force Field for Earth Material Modeling.

Phase Field Simulation for Non-Isothermal Solidification of Multicomponent Alloys Coupled with Thermodynamics Database

Development of suitable interatomic potentials for simulation of liquid and amorphous Cu–Zr alloys

Development and application of ReaxFF methodology for understanding the chemical dynamics of metal carbonates in aqueous solutions

Development of ReaxFF Reactive Force Field for Aqueous Iron–Sulfur Clusters with Applications to Stability and Reactivity in Water

Reaxff(Sio) Reactive Force Field For Silicon And Silicon Oxide Systems

An Efficient Reactive Force Field without Explicit Coordination Dependence for Studying Caustic Aluminum Chemistry

Interatomic force fields for zirconium based on the embedded atom method and the tabulated Gaussian Approximation Potential

Development of a Ni–Al reactive force field for Ni-based superalloy: revealing electrostatic effects on mechanical deformation

Optimization and Validation of a Deep Learning CuZr Atomistic Potential: Robust Applications for Crystalline and Amorphous Phases with near-DFT Accuracy

ReaxFF reactive force field for the Y-doped BaZrO3 proton conductor with applications to diffusion rates for multigranular systems

An investigation of the combined size and rate effects on the mechanical responses of FCC metals

Melting Point and Crystal Growth Kinetics of Metals and Metal Oxides Using Reactive Force Fields: The Case of Aluminum and Alumina

Molecular Dynamics Simulations of Melting Iron Nanoparticles with/without Defects Using a Reaxff Reactive Force Field

Parameterization of a reactive force field using a Monte Carlo algorithm

Analysis of semi-empirical interatomic potentials appropriate for simulation of crystalline and liquid Al and Cu

Optimization of Reactive Force Field Simulation: Refactor, Parallelization, and Vectorization for Interactions

Development and validation of a ReaxFF reactive force field for Fe/Al/Ni alloys: molecular dynamics study of elastic constants, diffusion, and segregation

ReaxFF Reactive Force Field Development for Cu/Si Systems and application to Copper Cluster Formation During Cu Diffusion Inside Silicon