(Keynote) An Automatic Rapid Screening Tool for Fast Ionic Conductors

Abstract

Identifying new materials that combine high ionic conductivity with structural and electrochemical stability so far remains a slow trial and error search process. To rationally accelerate materials design and exploit the opportunities in the materials genome a dependable rapid screening of materials is required that can pre-select structures that merit higher level computational as well as experimental characterization. Here we report on the progress of our “softBV” bond-valence site energy-based automated pathway analysis utilizing our new Bond Valence Pathway Analyzer (BVPA) – with fast bond valence site energy calculations to quickly obtain suggested candidate materials for fast ionic conduction. BVPA provides rapid and simplified visualization, in order to bridge the gap between experimentalists and simulation [1,2]. Calculation of ion transport pathways can be done extremely quickly on the order of seconds or minutes on desktop PCs providing a speedup factor of 3 to 5 orders of magnitude compared to DFT-based NEB methods. Combined with a graphical user interface our software suite (that can be downloaded from [2] and is free for academic use) should enable experimentalists to quickly identify candidate solid electrolyte materials. We also aim to integrate the pre-screening into an automated workflow for subsequent DFT characterization [3].Results will be benchmarked against both experimental and DFT NEB migration barriers. Besides the migration barriers the approach now also comprises an AI-based dopant predictor utilizing bond-valence-based crystal chemical descriptors to assist experimentalists in exploring favorable substitutional doping strategies. We will also compare the predictability of absolute room temperature conductivities from static energy landscape analysis, bond-valence based empirical MD simulations and ab initio molecular dynamics (AIMD) simulations. While for small fast-ion conductor structures at sufficiently high temperatures AIMD appears to be the gold standard, the less reliable but computationally empirical approaches have an advantage in modelling complex disordered interfaces at low temperatures over longer periods. This eliminates the hazards involved in extrapolations down to room temperature properties for the frequent cases of order-disorder phase transitions at intermediate temperatures.As an example we will discuss lithium and sodium compounds containing multiple anions, in particular the combination of thiophosphate and halide anions or various MS4 polyanions. Based on computational screening using our bond valence site approach and DFT studies several thiophosphate halides along the A3PS4-LiX (Cl, Br, I; A = Li, Na) tie line [4] and the Ax(MS4)y(M’S4)z phase space [5] have been explored and their properties discussed based on BVSE pathway models and molecular dynamics simulations in combination with experimental (X-ray and neutron) diffraction, solid state NMR and electrochemical characterisation. MD simulations e.g. show that Li5(PS4)Cl2 is found to undergo an order-disorder phase transition and thus should, contrasting to earlier predictions, not be a fast Li+ ion conductor. The newly predicted thermodynamically stable cubic solid electrolyte Li15(PS4)4Cl3 was successfully prepared and characterized. Though its conductivity does not reach a superionic level, it demonstrates that the computational approach can successfully predict a completely new classes of solid electrolytes and can predict its optimization by doping. The simplicity of the approach also facilitates the study of homogeneity ranges as exemplified for the solid solution systems Li4-xPS4Ix (0<x<0.67) and Na9+x(MS4)3-x(SnS4)x (M = P, Sb; x ≈ 2).