Search for: All records

This work presents a comprehensive computational study showing how aliovalent doping, crystal structure, and oxygen vacancy interactions impact the oxygen vacancy conductivity of lanthanum strontium ferrite (LSF) as a function of temperature in air. First, density functional theory (DFT) calculations were performed to obtain the oxygen vacancy migration barriers and understand the oxidation state changes on neighboring Fe atoms during oxygen vacancy migration. The oxygen migration barrier energy and the corresponding diffusion coefficient were then combined with previously determined mobile oxygen vacancy concentrations to predict the overall oxygen vacancy conductivity and compare it with experimentally measured values. More importantly, the impact of phase changes, the La/Sr ratio, and the oxygen non-stoichiometry on the mobile oxygen vacancy concentration, diffusivity, and conductivity were analyzed. It was found that stabilizing rhombohedral LSF or cubic SFO (through doping or other means), such that oxygen-vacancy-ordering-induced phase transitions are prevented, leads to high oxygen conductivity under solid oxide fuel cell operating conditions. 

Accurate characterization of chemical strain is required to study a broad range of chemical–mechanical coupling phenomena. One of the most studied mechano-chemically active oxides, nonstoichiometric ceria (CeO 2−δ ), has only been described by a scalar chemical strain assuming isotropic deformation. However, combined density functional theory (DFT) calculations and elastic dipole tensor theory reveal that both the short-range bond distortions surrounding an oxygen-vacancy and the long-range chemical strain are anisotropic in cubic CeO 2−δ . The origin of this anisotropy is the charge disproportionation between the four cerium atoms around each oxygen-vacancy (two become Ce 3+ and two become Ce 4+ ) when a neutral oxygen-vacancy is formed. Around the oxygen-vacancy, six of the Ce 3+ –O bonds elongate, one of the Ce 3+ –O bond shorten, and all seven of the Ce 4+ –O bonds shorten. Further, the average and maximum chemical strain values obtained through tensor analysis successfully bound the various experimental data. Lastly, the anisotropic, oxygen-vacancy-elastic-dipole induced chemical strain is polarizable, which provides a physical model for the giant electrostriction recently discovered in doped and non-doped CeO 2−δ . Together, this work highlights the need to consider anisotropic tensors when calculating the chemical strain induced by dilute point defects in all materials, regardless of their symmetry. 

The solid electrolyte interphase (SEI) forms on electrode surfaces from decomposition of the electrolyte. However, there is almost no atomistic detail of SEI formation on Li metal anode, a major obstacle in understanding the highly complex battery electrochemistry sufficiently to design high performance batteries. Herein, a realistic atomistic model (39 000 atoms) for the SEI formation at the interface between the Li metal anode and ionic liquid electrolyte using reactive molecular dynamics simulations is provided. A ≈10 nm thick SEI composed of a dense ordered inorganic layer near the Li‐metal anode and a porous organic layer near the electrolyte is found. These results provide new insights into a deeper understanding of the complex SEI that should be useful in developing a new generation of highly efficient batteries.

 

Both aliovalent doping and the charge state of multivalent lattice ions determine the oxygen non-stoichiometry ( δ ) of mixed ionic and electronic conductors (MIECs). Unfortunately, it has been challenging for both modeling and experiments to determine the multivalent ion charge states in MIECs. Here, the Fe charge state distribution was determined for various compositions and phases of the MIEC La 1−x Sr x FeO 3−δ (LSF) using the spin-polarized density functional theory (DFT)-predicted magnetic moments on Fe. It was found that electron occupancy and crystal-field-splitting-induced differences between the Fe 3d-orbitals of the square pyramidally coordinated, oxygen-vacancy-adjacent Fe atoms and the octahedrally-coordinated, oxygen-vacancy-distant-Fe atoms determined whether the excess electrons produced during oxygen vacancy formation remained localized at the first nearest neighbor Fe atoms (resulting in small oxygen vacancy polarons, as in LaFeO 3 ) or were distributed to the second-nearest-neighbor Fe atoms (resulting in large oxygen vacancy polarons, as in SrFeO 3 ). The progressively larger polaron size and anisotropic shape changes with increasing Sr resulted in increasing oxygen vacancy interactions, as indicated by an increase in the oxygen vacancy formation energy above a critical δ threshold. This was consistent with experimental results showing that Sr-rich LSF and highly oxygen deficient compositions are prone to oxygen-vacancy-ordering-induced phase transformations, while Sr-poor and oxygen-rich LSF compositions are not. Since oxygen vacancy induced phase transformations cause a decrease in the mobile oxygen vacancy site fraction ( X ), both δ and X were predicted as a function of temperature and oxygen partial pressure, for multiple LSF compositions and phases using a combined thermodynamics and DFT approach.