The lithiation process and Li diffusion in amorphous SiO2 and Si from first-principles

Silicon is considered the next-generation, high-capacity anode for Li-ion energy storage applications, however, despite significant effort, there are still uncertainties regarding the bulk Si and surface SiO2 structural and chemical evolution as it undergoes lithiation and amorphization. In this paper, we present first-principles calculations of the evolution of the amorphous Si anode, including its oxide surface layer, as a function of Li concentration. We benchmark our methodology by comparing the results for the Si bulk to existing experimental evidence of local structure evolution, ionic diffusivity as well as electrochemical activity. Recognizing the important role of the surface Si oxide (either native or artificially grown), we undertake the same calculations for amorphous SiO2, analyzing its potential impact on the activity of Si anode materials. Derived voltage curves for the amorphous phases compare well to experimental results, highlighting that SiO2 lithiates at approximately 0.7 V higher than Si in the low Li concentration regime, which provides an important electrochemical fingerprint. The combined evidence suggests that i) the inherent diffusivity of amorphous Si is high (in the order 10−9cm2s−1 - 10−7cm2s−1), ii) SiO2 is thermodynamically driven to lithiate, such that Li–O local environments are increasingly favored as compared to Si–O bonding, iii) the ionic diffusivity of Li in LiySiO2 is initially two orders of magnitude lower than that of LiySi at low Li concentrations but increases rapidly with increasing Li content and iv) the final lithiation product of SiO2 is Li2O and highly lithiated silicides. Hence, this work suggests that - excluding explicit interactions with the electrolyte - the SiO2 surface layer presents a kinetic impediment for the lithiation of Si and a sink for Li inventory, resulting in non-reversible capacity loss through strong local Li–O bond formation

Crystallographic map: A general lattice and basis formalism enabling efficient and discretized exploration of crystallographic phase space

Three-dimensional lattices are fundamental to solid-state physics. The description of a lattice with an atomic basis constitutes the necessary information to predict solid phase properties and evolution. Here, we present an algorithm for systematically exploring crystallographic phase space. Coupled with ab initio techniques, such as density functional theory, this algorithm offers an approach for exploring and tuning materials behavior, with a broad range of potential applications: particularly martensitic phase transformations and materials stability

A method to computationally screen for tunable properties of crystalline alloys

Conventionally, high-throughput computational materials searches start from an input set of bulk compounds extracted from material databases, and this set is screened for candidate materials for specific applications. In contrast, many functional materials, and especially semiconductors, are heavily engineered alloys or solid solutions of multiple compounds rather than a single bulk compound. To improve our ability to design functional materials, in this work we propose a framework and open-source code to automatically construct possible "alloy pairs" and "alloy systems" and detect "alloy members" from a set of existing, experimental or calculated ordered compounds, without requiring any additional metadata beyond their crystal structure. We provide analysis tools to estimate stability across each alloy. As a demonstration, we apply this framework to all inorganic materials in the Materials Project database to create a new database of over 600,000 unique alloy pair entries that can then be used in materials discovery studies to search for materials with tunable properties. This new database has been incorporated into the Materials Project website and linked with corresponding material identifiers for any user to query and explore. Using an example of screening for p-type transparent conducting materials, we demonstrate how using this methodology reveals candidate material systems that might otherwise have been excluded by a traditional screening. This work lays a foundation from which materials databases can go beyond stoichiometric compounds, and approach a more realistic description of compositionally tunable materials

Unsupervised word embeddings capture latent knowledge from materials science literature.

The overwhelming majority of scientific knowledge is published as text, which is difficult to analyse by either traditional statistical analysis or modern machine learning methods. By contrast, the main source of machine-interpretable data for the materials research community has come from structured property databases1,2, which encompass only a small fraction of the knowledge present in the research literature. Beyond property values, publications contain valuable knowledge regarding the connections and relationships between data items as interpreted by the authors. To improve the identification and use of this knowledge, several studies have focused on the retrieval of information from scientific literature using supervised natural language processing3-10, which requires large hand-labelled datasets for training. Here we show that materials science knowledge present in the published literature can be efficiently encoded as information-dense word embeddings11-13 (vector representations of words) without human labelling or supervision. Without any explicit insertion of chemical knowledge, these embeddings capture complex materials science concepts such as the underlying structure of the periodic table and structure-property relationships in materials. Furthermore, we demonstrate that an unsupervised method can recommend materials for functional applications several years before their discovery. This suggests that latent knowledge regarding future discoveries is to a large extent embedded in past publications. Our findings highlight the possibility of extracting knowledge and relationships from the massive body of scientific literature in a collective manner, and point towards a generalized approach to the mining of scientific literature

Gold–Thiolate Nanocluster Dynamics and Intercluster Reactions Enabled by a Machine Learned Interatomic Potential

Monolayer protected metal clusters comprise a rich class of molecular systems and are promising candidate materials for a variety of applications. While a growing number of protected nanoclusters have been synthesized and characterized in crystalline forms, their dynamical behavior in solution, including prenucleation cluster formation, is not well understood due to limitations both in characterization and first-principles modeling techniques. Recent advancements in machine-learned interatomic potentials are rapidly enabling the study of complex interactions such as dynamical behavior and reactivity on the nanoscale. Here, we develop an Au-S-C-H atomic cluster expansion (ACE) interatomic potential for efficient and accurate molecular dynamics simulations of thiolate-protected gold nanoclusters (Aun(SCH3)m). Trained on more than 30,000 density functional theory calculations of gold nanoclusters, the interatomic potential exhibits ab initio level accuracy in energies and forces and replicates nanocluster dynamics including thermal vibration and chiral inversion. Long dynamics simulations (up to 0.1 μs time scale) reveal a mechanism explaining the thermal instability of neutral Au25(SR)18 clusters. Specifically, we observe multiple stages of isomerization of the Au25(SR)18 cluster, including a chiral isomer. Additionally, we simulate coalescence of two Au25(SR)18 clusters and observe series of clusters where the formation mechanisms are critically mediated by ligand exchange in the form of [Au-S]n rings

High-throughput optical absorption spectra for inorganic semiconductors

An optical absorption spectrum constitutes one of the most fundamental material characteristics, with relevant applications ranging from material identification to energy harvesting and optoelectronics. However, the database of both experimental and computational spectra is currently lacking. In this study, we designed a computational workflow for the optical absorption spectrum and integrated the simulated spectra into the Materials Project. Using density-functional theory, we computed the frequency-dependent dielectric function and the corresponding absorption coefficient for more than 1000 solid compounds of varying crystal structure and chemistry. The computed spectra show excellent agreement, as quantified by a high value of the Pearson correlation, with experimental results when applying the band gap correction from the HSE functional. The demonstrated calculated accuracy in the spectra suggests that the workflow can be applied in screening studies for materials with specific optical properties

Elucidating the Structure of the Magnesium Aluminum Chloride Complex electrolyte for Magnesium-ion batteries

We present a rigorous analysis of the Magnesium Aluminum Chloro Complex (MACC) in tetrahydrofuran (THF), one of the few electrolytes that can reversibly plate and strip Mg. We use \emph{ab initio} calculations and classical molecular dynamics simulations to interrogate the MACC electrolyte composition with the goal of addressing two urgent questions that have puzzled battery researchers: \emph{i}) the functional species of the electrolyte, and \emph{ii}) the complex equilibria regulating the MACC speciation after prolonged electrochemical cycling, a process termed as conditioning, and after prolonged inactivity, a process called aging. A general computational strategy to untangle the complex structure of electrolytes, ionic liquids and other liquid media is presented. The analysis of formation energies and grand-potential phase diagrams of Mg-Al-Cl-THF suggests that the MACC electrolyte bears a simple chemical structure with few simple constituents, namely the electro-active species MgCl$^+$ and AlCl$_4^-$ in equilibrium with MgCl$_2$ and AlCl$_3$. Knowledge of the stable species of the MACC electrolyte allows us to determine the most important equilibria occurring during electrochemical cycling. We observe that Al deposition is always preferred to Mg deposition, explaining why freshly synthesized MACC cannot operate and needs to undergo preparatory conditioning. Similarly, we suggest that aluminum displacement and depletion from the solution upon electrolyte resting (along with continuous MgCl$_2$ regeneration) represents one of the causes of electrolyte aging. Finally, we compute the NMR shifts from shielding tensors of selected molecules and ions providing fingerprints to guide future experimental investigations

Conformational Entropy as a Means to Control the Behavior of Poly(diketoenamine) Vitrimers In and Out of Equilibrium.

Control of equilibrium and non-equilibrium thermomechanical behavior of poly(diketoenamine) vitrimers is shown by incorporating linear polymer segments varying in molecular weight (MW) and conformational degrees of freedom into the dynamic covalent network. While increasing MW of linear segments yields a lower storage modulus at the rubbery plateau after softening above the glass transition (Tg ), both Tg and the characteristic time of stress relaxation are independently governed by the conformational entropy of the embodied linear segments. Activation energies for bond exchange in the solid state are lower for networks incorporating flexible chains; the network topology freezing temperature decreases with increasing MW of flexible linear segments but increases with increasing MW of stiff segments. Vitrimer reconfigurability is therefore influenced not only by the energetics of bond exchange for a given network density, but also the entropy of polymer chains within the network