Source author record

Gerbrand Ceder

Gerbrand Ceder appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

29works
11topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

29 published item(s)

preprint2026arXiv

Li+/H+ exchange in solid-state oxide Li-ion conductors

Understanding the moisture stability of oxide Li-ion conductors is important for their practical applications in solid-state batteries. Unlike sulfide or halide conductors, oxide conductors generally better resist degradation when in contact with water, but can still undergo topotactic \ch{Li+}/\ch{H+} exchange (LHX). Here, we combine density functional theory (DFT) calculations with a machine-learning interatomic potential model to investigate the thermodynamic driving force of the LHX reaction for two representative oxide Li-ion conductor families: garnets and NASICONs. Li-stuffed garnets exhibit a strong driving force for proton exchange due to their high Li chemical potential. In contrast, NASICONs demonstrate a higher resistance against proton exchange due to the lower Li chemical potential and the lower O-H bond covalency for polyanion-bonded oxygens. Our findings reveal a critical trade-off: Li stuffing enhances conductivity but increases moisture susceptibility. This study underscores the importance of designing Li-ion conductors that possess both high conductivity and high stability in practical environments.

preprint2026arXiv

Mechanisms of alkali ionic transport in amorphous oxyhalides solid state conductors

Amorphous oxyhalides have attracted significant attention due to their relatively high ionic conductivity ($>$1 mS cm$^{-1}$), excellent chemical stability, mechanical softness, and facile synthesis routes via standard solid-state reactions. These materials exhibit an ionic conductivity that is almost independent of the underlying chemistry, in stark contrast to what occurs in crystalline conductors. In this work, we employ an accurately fine-tuned machine learning interatomic potential to construct large-scale molecular dynamics trajectories encompassing hundreds of nanoseconds to obtain statistically converged transport properties. We find that the amorphous state consists of chain fragments of metal-anion tetrahedra of various lenght. By analyzing the residence time of alkali cations migrating around tetrahedrally-coordinated trivalent metal ions, we find that oxygen anions on the metal-anion tetrahedra limit alkali diffusion. By computing the full Einstein expression of the ionic conductivity, we demonstrate that the alkali transference number of these materials is strongly influenced by distinct-particles correlations, while at the same time they are characterized by an alkali Haven ratio close to one, implying that ionic transport is largely dictated by uncorrelated self-diffusion. Finally, by extending this analysis to chemical compositions $AMX_{2.5}\textsf{O}_{0.75}$, spanning different alkaline ($A$ = Li, Na, K), metallic ($M$ = Al, Ga, In), and halogen ($X$ = Cl, Br, I) species, we clarify why the diffusion properties of these materials remain largely insensitive to variations in atomic chemistry.

preprint2023arXiv

The cluster decomposition of the configurational energy of multicomponent alloys

Lattice models parameterized using first-principles calculations constitute an effective framework to simulate the thermodynamic behavior of physical systems. The cluster expansion method is a flexible lattice-based method used extensively in the study of multicomponent alloys. Yet despite its prevalent use, a well-defined understanding of expansion terms has remained elusive. In this letter, we introduce the cluster decomposition as a unique and basis-agnostic decomposition of any general function of the atomic configuration in a crystal. We demonstrate that cluster expansions constructed from arbitrary orthonormal basis sets are all representations of the same cluster decomposition. We show how the norms of expansion coefficients associated with the same crystallographic orbit are invariant to changes between orthonormal bases. Based on its uniqueness and orthogonality properties, we identify the cluster decomposition as an invariant ANOVA decomposition. We leverage these results to illustrate how functional analysis of variance and sensitivity analysis can be used to directly interpret interactions among species and gain insight into computed thermodynamic properties. The work we present in this letter opens new directions for parameter estimation, interpretation, and use of applied lattice models on well-established mathematical and statistical grounds.

preprint2022arXiv

An $\ell_0\ell_2$-norm regularized regression model for construction of robust cluster expansions in multicomponent systems

We introduce the $\ell_0\ell_2$-norm regularization and hierarchy constraints into linear regression for the construction of cluster expansion to describe configurational disorder in materials. The approach is implemented through mixed integer quadratic programming (MIQP). The $\ell_2$-norm regularization is used to suppress intrinsic data noise, while $\ell_0$-norm is used to penalize the number of non-zero elements in the solution. The hierarchy relation between clusters imposes relevant physics and is naturally included by the MIQP paradigm. As such, sparseness and cluster hierarchy can be well optimized to obtain a robust, converged, and effective cluster interactions with improved physical meaning. We demonstrate the effectiveness of $\ell_0\ell_2$-norm regularization in two high-component disordered rocksalt cathode material systems, where we compare the cross-validation and convergence speed, reproduction of phase diagrams, voltage profiles, and Li-occupancy energies with those of the conventional $\ell_1$-norm regularized cluster expansion model.

preprint2022arXiv

Dataset of gold nanoparticle sizes and morphologies extracted from literature-mined microscopy images

The factors controlling the size and morphology of nanoparticles have so far been poorly understood. Data-driven techniques are an exciting avenue to explore this field through the identification of trends and correlations in data. However, for these techniques to be utilized, large datasets annotated with the structural attributes of nanoparticles are required. While experimental SEM/TEM images collected from controlled experiments are reliable sources of this information, large-scale collection of these images across a variety of experimental conditions is expensive and infeasible. Published scientific literature, which provides a vast source of high-quality figures including SEM/TEM images, can provide a large amount of data at a lower cost if effectively mined. In this work, we develop an automated pipeline to retrieve and analyse microscopy images from gold nanoparticle literature and provide a dataset of 4361 SEM/TEM images of gold nanoparticles along with automatically extracted size and morphology information. The dataset can be queried to obtain information about the physical attributes of gold nanoparticles and their statistical distributions.

preprint2022arXiv

Machine-learning rationalization and prediction of solid-state synthesis conditions

There currently exist no quantitative methods to determine the appropriate conditions for solid-state synthesis. This not only hinders the experimental realization of novel materials but also complicates the interpretation and understanding of solid-state reaction mechanisms. Here, we demonstrate a machine-learning approach that predicts synthesis conditions using large solid-state synthesis datasets text-mined from scientific journal articles. Using feature importance ranking analysis, we discovered that optimal heating temperatures have strong correlations with the stability of precursor materials quantified using melting points and formation energies ($ΔG_f$, $ΔH_f$). In contrast, features derived from the thermodynamics of synthesis-related reactions did not directly correlate to the chosen heating temperatures. This correlation between optimal solid-state heating temperature and precursor stability extends Tamman's rule from intermetallics to oxide systems, suggesting the importance of reaction kinetics in determining synthesis conditions. Heating times are shown to be strongly correlated with the chosen experimental procedures and instrument setups, which may be indicative of human bias in the dataset. Using these predictive features, we constructed machine-learning models with good performance and general applicability to predict the conditions required to synthesize diverse chemical systems. Codes and data used in this work can be found at: https://github.com/CederGroupHub/s4.

preprint2022arXiv

Text-mined dataset of gold nanoparticle synthesis procedures, morphologies, and size entities

Gold nanoparticles are highly desired for a range of technological applications due to their tunable properties, which are dictated by the size and shape of the constituent particles. Many heuristic methods for controlling the morphological characteristics of gold nanoparticles are well known. However, the underlying mechanisms controlling their size and shape remain poorly understood, partly due to the immense range of possible combinations of synthesis parameters. Data-driven methods can offer insight to help guide understanding of these underlying mechanisms, so long as sufficient synthesis data are available. To facilitate data mining in this direction, we have constructed and made publicly available a dataset of codified gold nanoparticle synthesis protocols and outcomes extracted directly from the nanoparticle materials science literature using natural language processing and text-mining techniques. This dataset contains 5,154 data records, each representing a single gold nanoparticle synthesis article, filtered from a database of 4,973,165 publications. Each record contains codified synthesis protocols and extracted morphological information from a total of 7,608 experimental and 12,519 characterization paragraphs.

preprint2022arXiv

ULSA: Unified Language of Synthesis Actions for Representation of Synthesis Protocols

Applying AI power to predict syntheses of novel materials requires high-quality, large-scale datasets. Extraction of synthesis information from scientific publications is still challenging, especially for extracting synthesis actions, because of the lack of a comprehensive labeled dataset using a solid, robust, and well-established ontology for describing synthesis procedures. In this work, we propose the first Unified Language of Synthesis Actions (ULSA) for describing ceramics synthesis procedures. We created a dataset of 3,040 synthesis procedures annotated by domain experts according to the proposed ULSA scheme. To demonstrate the capabilities of ULSA, we built a neural network-based model to map arbitrary ceramics synthesis paragraphs into ULSA and used it to construct synthesis flowcharts for synthesis procedures. Analysis for the flowcharts showed that (a) ULSA covers essential vocabulary used by researchers when describing synthesis procedures and (b) it can capture important features of synthesis protocols. This work is an important step towards creating a synthesis ontology and a solid foundation for autonomous robotic synthesis.

preprint2021arXiv

Sparse expansions of multicomponent oxide configuration energy using coherency & redundancy

Compressed sensing has become a widely accepted paradigm to construct high dimensional cluster expansion models used for statistical mechanical studies of atomic configuration in complex multicomponent crystalline materials. However, strict sampling requirements necessary to obtain minimal coherence measurements for compressed sensing to guarantee accurate estimation of model parameters are difficult and in some cases impossible to satisfy due to the inability of physical systems to access certain configurations. Nevertheless, the dependence of energy on atomic configuration can still be adequately learned without these strict requirements by using compressed sensing by way of coherent measurements using redundant function sets known as frames. We develop a particular frame constructed from the union of all occupancy-based cluster expansion basis sets. We illustrate how using this highly redundant frame yields sparse expansions of the configuration energy of complex oxide materials that are competitive and often surpass the prediction accuracy and sparsity of models obtained from standard cluster expansions.

preprint2020arXiv

A critical examination of compound stability predictions from machine-learned formation energies

Machine learning has emerged as a novel tool for the efficient prediction of materials properties, and claims have been made that machine-learned models for the formation energy of compounds can approach the accuracy of Density Functional Theory (DFT). The models tested in this work include five recently published compositional models, a baseline model using stoichiometry alone, and a structural model. By testing seven machine learning models for formation energy on stability predictions using the Materials Project database of DFT calculations for 85,014 unique chemical compositions, we show that while formation energies can indeed be predicted well, all compositional models perform poorly on predicting the stability of compounds, making them considerably less useful than DFT for the discovery and design of new solids. Most critically, in sparse chemical spaces where few stoichiometries have stable compounds, only the structural model is capable of efficiently detecting which materials are stable. The non-incremental improvement of structural models compared with compositional models is noteworthy and encourages the use of structural models for materials discovery, with the constraint that for any new composition, the ground-state structure is not known a priori. This work demonstrates that accurate predictions of formation energy do not imply accurate predictions of stability, emphasizing the importance of assessing model performance on stability predictions, for which we provide a set of publicly available tests.

preprint2020arXiv

Modeling of Electro-deposition and Mechanical Stability at Li Metal/Solid Electrolyte Interface during Plating in Solid-State Batteries

Interfacial deposition stability between Li metal and a solid electrolyte (SE) is important in preventing interfacial contact loss, mechanical fracture, and dendrite growth in Li-metal solid-state batteries (SSB). In this work, we investigate the deposition and mechanical stability at the Li metal/SE interface and its consequences (such as SE fracture and contact loss). A wide range of contributing factors are investigated, such as charge and mass transfer kinetics, the plasticity of Li metal and fracture of the SE, and the applied stack pressure. We quantify the effect of the ionic conductivity of the SE, the exchange current density of the interfacial charge-transfer reaction and SE surface roughness on the Li deposition stability at the Li metal/SE interface. We also propose a mechanical stability window for the applied stack pressure that can prevent both contact loss and SE fracture, which can be extended to other metal-electrode (such as Sodium) SSB systems.

preprint2020arXiv

Similarity of Precursors in Solid-state Synthesis as Text-Mined from Scientific Literature

Collecting and analyzing the vast amount of information available in the solid-state chemistry literature may accelerate our understanding of materials synthesis. However, one major problem is the difficulty of identifying which materials from a synthesis paragraph are precursors or are target materials. In this study, we developed a two-step Chemical Named Entity Recognition (CNER) model to identify precursors and targets, based on information from the context around material entities. Using the extracted data, we conducted a meta-analysis to study the similarities and differences between precursors in the context of solid-state synthesis. To quantify precursor similarity, we built a substitution model to calculate the viability of substituting one precursor with another while retaining the target. From a hierarchical clustering of the precursors, we demonstrate that "chemical similarity" of precursors can be extracted from text data. Quantifying the similarity of precursors helps provide a foundation for suggesting candidate reactants in a predictive synthesis model.

preprint2016arXiv

Constructing and proving the ground state of a generalized Ising model by the cluster tree optimization algorithm

Generalized Ising models, also known as cluster expansions, are an important tool in many areas of condensed-matter physics and materials science, as they are often used in the study of lattice thermodynamics, solid-solid phase transitions, magnetic and thermal properties of solids, and fluid mechanics. However, the problem of finding the global ground state of generalized Ising model has remained unresolved, with only a limited number of results for simple systems known. We propose a method to efficiently find the periodic ground state of a generalized Ising model of arbitrary complexity by a new algorithm which we term cluster tree optimization. Importantly, we are able to show that even in the case of an aperiodic ground state, our algorithm produces a sequence of states with energy converging to the true ground state energy, with a provable bound on error. Compared to the current state-of-the-art polytope method, this algorithm eliminates the necessity of introducing an exponential number of variables to counter frustration, and thus significantly improves tractability. We believe that the cluster tree algorithm offers an intuitive and efficient approach to finding and proving ground states of generalized Ising Hamiltonians of arbitrary complexity, which will help validate assumptions regarding local vs. global optimality in lattice models, as well as offer insights into the low-energy behavior of highly frustrated systems.

preprint2016arXiv

Finding and proving the exact ground state of a generalized Ising model by convex optimization and MAX-SAT

Lattice models, also known as generalized Ising models or cluster expansions, are widely used in many areas of science and are routinely applied to alloy thermodynamics, solid-solid phase transitions, magnetic and thermal properties of solids, and fluid mechanics, among others. However, the problem of finding the true global ground state of a lattice model, which is essential for all of the aforementioned applications, has remained unresolved, with only a limited number of results for highly simplified systems known. In this article, we present the first general algorithm to find the exact ground states of complex lattice models and to prove their global optimality, resolving this fundamental problem in condensed matter and materials theory. We transform the infinite-discrete-optimization problem into a pair of combinatorial optimization (MAX-SAT) and non-smooth convex optimization (MAX-MIN) problems, which provide upper and lower bounds on the ground state energy respectively. By systematically converging these bounds to each other, we find and prove the exact ground state of realistic Hamiltonians whose solutions are completely intractable via traditional methods. Considering that currently such Hamiltonians are solved using simulated annealing and genetic algorithms that are often unable to find the true global energy minimum, and never able to prove the optimality of their result, our work opens the door to resolving long-standing uncertainties in lattice models of physical phenomena.

preprint2016arXiv

Impact of Intermediate Sites on Bulk Diffusion Barriers: Mg Intercalation in Mg$_2$Mo$_3$O$_8$

The ongoing search for high voltage positive electrode materials for Mg batteries has been primarily hampered by poor Mg mobility in bulk oxide frameworks. Motivated by the presence of Mo$_3$ clusters that can facilitate charge redistribution and the presence of Mg in a non-preferred (tetrahedral) coordination environment, we have investigated the Mg (de)intercalation behavior in layered-Mg$_2$Mo$_3$O$_8$, a potential positive electrode. While no electrochemical activity is observed, chemical demagnesiation of Mg$_2$Mo$_3$O$_8$ is successful but leads to amorphization. Subsequent first-principles calculations predict a strong thermodynamic driving force for structure decomposition at low Mg concentrations and high activation barriers for bulk Mg diffusion, in agreement with experimental observations. Further analysis of the Mg diffusion pathway reveals an O--Mg--O dumbbell intermediate site that creates a high Mg$^{2+}$ migration barrier, indicating the influence of transition states on setting the magnitude of migration barriers.

preprint2016arXiv

Role of structural H$_2$O in intercalation electrodes: the case of Mg in nano-crystalline Xerogel-V$_2$O$_5$

Co-intercalation is a potential approach to influence the voltage and mobility with which cations insert in electrodes for energy storage devices. Combining a robust thermodynamic model with first-principles calculations, we present a detailed investigation revealing the important role of H$_2$O during ion intercalation in nano-materials. We examine the scenario of Mg$^{2+}$ and H$_2$O co-intercalation in nano-crystalline Xerogel-V$_2$O$_5$, a potential cathode material to achieve energy density greater than Li-ion batteries. Water co-intercalation in cathode materials could broadly impact an electrochemical system by influencing its voltages or causing passivation at the anode. The analysis of the stable phases of Mg-Xerogel V$_2$O$_5$ and voltages at different electrolytic conditions reveals a range of concentrations for Mg in the Xerogel and H$_2$O in the electrolyte where there is no thermodynamic driving force for H$_2$O to shuttle with Mg during electrochemical cycling. Also, we demonstrate that H$_2$O shuttling with the Mg$^{2+}$ ions in wet electrolytes yields higher voltages than in dry electrolytes. The thermodynamic framework used to study water and Mg$^{2+}$ co-intercalation in this work opens the door for studying the general phenomenon of solvent co-intercalation observed in other complex solvent-electrode pairs used in the Li- and Na-ion chemical spaces.

preprint2016arXiv

Synthesis of a mixed-valent tin nitride and considerations of its possible crystal structures

Recent advances in theoretical structure prediction methods and high-throughput computational techniques are revolutionizing experimental discovery of the thermodynamically stable inorganic materials. Metastable materials represent a new frontier for studies, since even simple binary non ground state compounds of common elements may be awaiting discovery. However, there are significant research challenges related to non-equilibrium thin film synthesis and crystal structure predictions, such as small strained crystals in the experimental samples and energy minimization based theoretical algorithms. Here we report on experimental synthesis and characterization, as well as theoretical first-principles calculations of a previously unreported mixed-valent binary tin nitride. Thin film experiments indicate that this novel material is N-deficient SnN with tin in the mixed II/IV valence state and a small low-symmetry unit cell. Theoretical calculations suggest that the most likely crystal structure has the space group 2 (SG2) related to the distorted delafossite (SG166), which is nearly 0.1 eV/atom above the ground state SnN polymorph. This observation is rationalized by the structural similarity of the SnN distorted delafossite to the chemically related Sn3N4 spinel compound, which provides a fresh scientific insight into the reasons for growth of polymorphs of the metastable material. In addition to reporting on the discovery of the simple binary SnN compound, this paper illustrates a possible way of combining a wide range of advanced characterization techniques with the first-principle property calculation methods, to elucidate the most likely crystal structure of the previously unreported metastable materials.

preprint2015arXiv

Calibrating transition metal energy levels and oxygen bands in first principles calculations: accurate prediction of redox potentials and charge transfer in lithium transition metal oxides

Transition metal (TM) oxides play an increasingly important role in technology today including applications such as catalysis, solar energy harvesting, and energy storage. In many of these applications, the details of their electronic structure near the Fermi level are critically important for their properties. We propose a first-principles based computational methodology for the accurate prediction of oxygen charge transfer in TM oxides and lithium TM (Li-TM) oxides. To obtain accurate electronic structures, the Heyd-Scuseria-Ernzerhof (HSE06) hybrid functional is adopted and the amount of exact Hartree-Fock exchange (mixing parameter) is adjusted to reproduce reference band gaps. We show that the HSE06 functional with optimal mixing parameter yields not only improved electronic densities of states but also better energetics (Li-intercalation voltages) for LiCoO2 and LiNiO2 as compared to GGA, GGA+U and standard HSE06. We find that the optimal mixing parameters for TM oxides are system-specific and correlate with the covalency (ionicity) of the TM species. Strong covalent (ionic) nature of TM-O bonding leads to lower (higher) optimal mixing parameters. We find that optimized HSE06 functionals predict stronger hybridization of the Co 3d and O 2p orbitals than GGA, resulting in a greater contribution from oxygen states to charge compensation upon delithiation in LiCoO2. We also find that the band gaps of Li-TM oxides increase linearly with the mixing parameter, enabling the straightforward determination of optimal mixing parameters based on GGA (α = 0.0) and HSE06 (α = 0.25) calculations. Our results also show that G0W0@GGA+U band gaps of TM oxides (MO, M = Mn, Co, Ni) and LiCoO2 agree well with experimental references, suggesting that G0W0 calculations can be used as a reference for the calibration of the mixing parameter in case no experimental band gap has been reported.

preprint2015arXiv

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.

preprint2015arXiv

The intercalation phase diagram of Mg in V$_2$O$_5$ from first principles

We have investigated Mg intercalation into orthorhombic V$_2$O$_5$, one of only three cathodes known to reversibly intercalate Mg ions. By calculating the ground state Mg$_x$V$_2$O$_5$ configurations and by developing a cluster expansion for the configurational disorder in $δ$-V$_2$O$_5$, a full temperature-composition phase diagram is derived. Our calculations indicate an equilibrium phase separating behavior between fully demagnesiated $α$-V$_2$O$_5$ and fully magnesiated $δ$-V$_2$O$_5$, but also motivate the existence of potentially metastable solid solution transformation paths in both phases. We find significantly better mobility for Mg in the $δ$ polymorph suggesting that better performance can be achieved by cycling Mg in the $δ$ phase.

preprint2015arXiv

Understanding the Initial Stages of Reversible Mg Deposition and Stripping in Inorganic Non-Aqueous Electrolytes

Multi-valent (MV) battery architectures based on pairing a Mg metal anode with a high-voltage ($\sim$ 3 V) intercalation cathode offer a realistic design pathway toward significantly surpassing the energy storage performance of traditional Li-ion based batteries, but there are currently only few electrolyte systems that support reversible Mg deposition. Using both static first-principles calculations and $ab\; initio$ molecular dynamics, we perform a comprehensive adsorption study of several salt and solvent species at the interface of Mg metal with an electrolyte of Mg$^{2+}$ and Cl$^-$ dissolved in liquid tetrahydrofuran (THF). Our findings not only provide a picture of the stable species at the interface, but also explain how this system can support reversible Mg deposition and as such we provide insights in how to design other electrolytes for Mg plating and stripping. The active depositing species are identified to be (MgCl)$^+$ monomers coordinated by THF, which exhibit preferential adsorption on Mg compared to possible passivating species (such as THF solvent or neutral MgCl$_2$ complexes). Upon deposition, the energy to desolvate these adsorbed complexes and facilitate charge-transfer is shown to be small ($\sim$ 61 $-$ 46.2 kJ mol$^{-1}$ to remove 3 THF from the strongest adsorbing complex), and the stable orientations of the adsorbed but desolvated (MgCl)$^+$ complexes appear favorable for charge-transfer. Finally, observations of Mg-Cl dissociation at the Mg surface at very low THF coordinations (0 and 1) suggest that deleterious Cl incorporation in the anode may occur upon plating. In the stripping process, this is beneficial by further facilitating the Mg removal reaction.

preprint2013arXiv

High-Throughput Screening of Perovskite Alloys for Piezoelectric Performance and Formability

We screen a large chemical space of perovskite alloys for systems with the right properties to accommodate a morphotropic phase boundary (MPB) in their composition-temperature phase diagram, a crucial feature for high piezoelectric performance. We start from alloy end-points previously identified in a high-throughput computational search. An interpolation scheme is used to estimate the relative energies between different perovskite distortions for alloy compositions with a minimum of computational effort. Suggested alloys are further screened for thermodynamic stability. The screening identifies alloy systems already known to host a MPB, and suggests a few new ones that may be promising candidates for future experiments. Our method of investigation may be extended to other perovskite systems, e.g., (oxy-)nitrides, and provides a useful methodology for any application of high-throughput screening of isovalent alloy systems.

preprint2011arXiv

Lithium Diffusion in Graphitic Carbon

Graphitic carbon is currently considered the state-of-the-art material for the negative electrode in lithium-ion cells, mainly due to its high reversibility and low operating potential. However, carbon anodes exhibit mediocre charge/discharge rate performance, which contributes to severe transport-induced surface-structural damage upon prolonged cycling, and limits the lifetime of the cell. Lithium bulk diffusion in graphitic carbon is not yet completely understood, partly due to the complexity of measuring bulk transport properties in finite-sized, non-isotropic particles. To solve this problem for graphite, we use the Devanathan-Stachurski electrochemical methodology combined with ab-initio computations to deconvolute, and quantify the mechanism of lithium-ion diffusion in highly oriented pyrolytic graphite (HOPG). The results reveal inherent high lithium-ion diffusivity in the direction parallel to the graphene plane (ca. 10^-7 - 10^-6 cm2 s-1), as compared to sluggish lithium-ion transport along grain boundaries (ca. 10^-11 cm^2 s^-1), indicating the possibility of rational design of carbonaceous materials and composite electrodes with very high rate capability.

preprint2010arXiv

First-principles determined charge and orbital interactions in Fe$_3$O$_4$

The interactions between charge and orbitally ordered $d$-electrons are important in many transition metal oxides. We propose an effective energy model for such interactions, parameterized with DFT+U calculations, so that energy contributions of both electronic and lattice origin can be simultaneously accounted for. The model is applied to the low-temperature phase of magnetite, for which we propose a new ground state structure. The effective interactions on the B-lattice of Fe$_3$O$_4$ can be interpreted in terms of electrostatics and short-range Kugel-Khomskii exchange coupling. The frustration between optimal charge and orbital orderings leads to a complex energy landscape whereby the supercell for the charge ordering, orbital ordering and ionic displacements can all be different.

preprint2006arXiv

Configurational electronic entropy and the phase diagram of mixed-valence oxides: the case of Li$_x$FePO$_4$

We demonstrate that configurational electronic entropy, previously neglected, in {\it ab initio} thermodynamics of materials can qualitatively modify the finite-temperature phase stability of mixed-valence oxides. While transformations from low-T ordered or immiscible states are almost always driven by configurational disorder (i.e. random occupation of lattice sites by multiple species), in FePO$_4$--LiFePO$_4$ the formation of a solid solution is almost entirely driven by electronic, rather than ionic configurational entropy. We argue that such an electronic entropic mechanism may be relevant to most other mixed-valence systems.

preprint2005arXiv

Coarse-graining protein energetics in sequence variables

We show that cluster expansions (CE), previously used to model solid-state materials with binary or ternary configurational disorder, can be extended to the protein design problem. We present a generalized CE framework, in which properties such as energy can be unambiguously expanded in the amino-acid sequence space. The CE coarse grains over nonsequence degrees of freedom (e.g., side-chain conformations) and thereby simplifies the problem of designing proteins, or predicting the compatibility of a sequence with a given structure, by many orders of magnitude. The CE is physically transparent, and can be evaluated through linear regression on the energies of training sequences. We show, as example, that good prediction accuracy is obtained with up to pairwise interactions for a coiled-coil backbone, and that triplet interactions are important in the energetics of a more globular zinc-finger backbone.

preprint2005arXiv

The electronic structure and band gap of LiFePO4 and LiMnPO4

Materials with the olivine LixMPO4 structure form an important new class of materials for rechargeable Li batteries. There is significant interest in their electronic properties because of the importance of electronic conductivity in batteries for high rate applications. The density of states of LixMPO4 (x = 0, 1 and M = Fe, Mn) has been determined with the ab initio GGA+U method, appropriate for these correlated electron systems. Computed results are compared with the optical gap of LiFePO4, as measured using UV-Vis-NIR diffuse reflectance spectroscopy. The results obtained from experiment (3.8-4.0 eV) and GGA+U computations (3.7 eV) are in very good agreement. However, standard GGA, without the same level of treatment of electron correlation, is shown to make large errors in predicting the electronic structure. It is argued that olivines are likely to be polaronic conductors with extrinsically determined carrier levels and that their electronic conductivity is therefore not simply related to the band gap.

preprint2005arXiv

The Li intercalation potential of LiMPO4 and LiMSiO4 olivines with M = Fe, Mn, Co, Ni

The Li intercalation potential of LiMPO4 and LiMSiO4 compounds with M = Fe, Mn, Co, and Ni is computed with the GGA+U method. It is found that this approach is considerably more accurate than standard LDA or GGA methods. The calculated potentials for LiFePO4, LiMnPO4 and LiCoOPO4 agree to within 0.1 V with experimental results. The LiNiPO4 potential is predicted to be above 5 V. The potentials of the silicate materials are all found to be rather high, but LiFeSiO4 and LiCoSiO4 have negligible volume change upon Li extraction.