Previous efforts to understand structure‐function relationships in high ionic conductivity materials for solid state batteries have predominantly relied on density functional theory (DFT‐) based ab initio molecular dynamics (MD). Such simulations, however, are computationally demanding and cannot be reasonably applied to large systems containing more than a hundred atoms. Here, an artificial neural network (ANN) is trained to accelerate the calculation of high accuracy atomic forces and energies used during such MD simulations. After carefully training a robust ANN for four and five element systems, nearly identical lithium ion diffusivities are obtained for Li10GeP2S12(LGPS) when benchmarking the ANN‐MD results with DFT‐MD. Applying the ANN‐MD approach, the effect of chlorine doping on the lithium diffusivity is calculated in an LGPS‐like structure and it is found that a dopant concentration of 1.3% maximizes ionic conductivity. The optimal concentration balances the competing consequences of effective atomic radii and dielectric constants on lithium diffusion and agrees with the experimental composition. Performing simulations at the resolution necessary to model experimentally relevant and optimal concentrations would be infeasible with traditional DFT‐MD. Systems that require a large number of simulated atoms can be studied more efficiently while maintaining high accuracy with the proposed ANN‐MD framework.
- PAR ID:
- 10406140
- Date Published:
- Journal Name:
- Frontiers in Materials
- Volume:
- 10
- ISSN:
- 2296-8016
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
more » « less
Mismatched complex oxide thin films and heterostructures have gained significant traction for use as electrolytes in intermediate temperature solid oxide fuel cells, wherein interfaces exhibit variation in ionic conductivity as compared to the bulk. Although misfit dislocations present at interfaces in these structures impact ionic conductivity, the fundamental mechanisms responsible for this effect are not well understood. To this end, a kinetic lattice Monte Carlo (KLMC) model was developed to trace oxygen vacancy diffusion at misfit dislocations in SrTiO3/BaZrO3 heterostructures and elucidate the atomistic mechanisms governing ionic diffusion at oxide interfaces. The KLMC model utilized oxygen vacancy migration energy barriers computed using molecular statics. While some interfaces promote oxygen vacancy diffusion, others impede their transport. Fundamental factors such as interface layer chemistry, misfit dislocation structure, and starting and ending sites of migrating ions play a crucial role in oxygen diffusivity. Molecular dynamics (MD) simulations were further performed to support qualitative trends for oxygen vacancy diffusion. Overall, the agreement between KLMC and MD is quite good, though MD tends to predict slightly higher conductivities, perhaps a reflection of nuanced structural relaxations that are not captured by KLMC. The current framework comprising KLMC modeling integrated with molecular statics offers a powerful tool to perform mechanistic studies focused on ionic transport in thin film oxide electrolytes and facilitate their rational design.
-
more » « less
Abstract Glass for pharmaceutical packaging requires high chemical durability for the safe storage and distribution of newly developed medicines. In borosilicate pharmaceutical glasses which typically contain a mixture of different modifier ions (alkali or alkaline earth), the dependence of the chemical durability on alkaline earth oxide concentrations is not well understood. Here, we have designed a series of borosilicate glasses with systematic substitutions of CaO with MgO while keeping their total concentrations at 13 mol% and a fixed Na2O concentration of 12.7 mol%. We used these glasses to investigate the influence of
R = [MgO]/([MgO] + [CaO]) on the resistance to aqueous corrosion at 80°C for 40 days. It was found that this type of borosilicate glass undergoes both leaching of modifier ions through an ion exchange process and etching of the glass network, leading to dissolution of the glass surface. Based on the concentration analysis of the Si and B species dissolved into the solution phase, the dissolved layer thickness was found to increase from ~100 to ~170 nm asR increases from 0 to 1. The depth profiling analysis of the glasses retrieved from the solution showed that the concentration of modifier ions (Na+, Ca2+, and Mg2+) at the interface between the solution and the corroded glass surface decreased to around 40%–60% of the corresponding bulk concentrations, regardless ofR and the leaching of modifier cations resulted in a silica‐rich layer in the surface. The leaching of Ca2+and Mg2+ions occurred within ~50 and <25 nm, respectively, from the glass surface and this thickness was not a strong function ofR . The leaching of Na+ions varied monotonically; the thickness of the Na+depletion layer increased from ~100 nm atR = 0 to ~200 nm atR = 1. Vibrational spectroscopy analysis suggested that the partial depletion of the ions may have caused some degree of the network re‐arrangement or re‐polymerization in the corroded layer. Overall, these results suggested that for the borosilicate glass, replacing [CaO] with [MgO] deteriorates the chemical durability in aqueous solution. -
null (Ed.)Induced-charge electro-osmotic (ICEO) flow caused by an alternating electric field applied around an infinitely long, ideally polarizable, uncharged circular cylinder in a binary electrolyte with unequal cation and anion diffusion coefficients is analysed. The thin-Debye-layer and weak-field approximations are invoked to compute the time-averaged, or rectified, quadrupolar ICEO flow around the cylinder. The inequality of ionic diffusion coefficients leads to transient ion concentration gradients, or concentration polarization, in the electroneutral bulk electrolyte outside the Debye layer. Consequently, the electric potential in the bulk is non-harmonic. Further, the concentration polarization alters the electro-osmotic slip at the surface of the cylinder and generates body forces in the bulk, both of which affect the rectified ICEO flow. Predictions for the strength of the rectified flow for varying ratio of ionic diffusion coefficients are in reasonable agreement with available experimental data. Our work highlights that an inequality in ionic diffusion coefficients – which all electrolytes possess to some extent – is an important factor in modelling ICEO flows.more » « less
-
Network glass structures are commonly characterized by the network formers and their linkage but modifiers can also play an important role on various features of glass structures. In this work, we investigated the effect of cation field strength (CFS) of common modifier cations with large differences of CFS on the structures of aluminoborosilicate glasses by performing molecular dynamics (MD) simulations with recently developed potentials. It was found that modifier cations with higher CFS such as Mg2+ significantly reduced the fraction of fourfold coordinated boron, suggesting that the cations with higher field strength favor nonbridging oxygen generation in the silicate network and are less effective for charge compensation. The findings from our MD simulations are compared with the results from NMR and Raman spectroscopy studies in the literature as well as those from other MD simulations. Insights of the CFS effect on glass structures and the structural role of Mg2+ ions are gained from these simulations results and related discussionsmore » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.3389/fmats.2023.1123213
