Stabilization of lithium metal in concentrated electrolytes: effects of electrode potential and solid electrolyte interphase formation.

Lithium (Li) metal negative electrodes have attracted wide attention for high-energy-density batteries. However, their low coulombic efficiency (CE) due to parasitic electrolyte reduction has been an alarming concern. Concentrated electrolytes are one of the promising concepts that can stabilize the Li metal/electrolyte interface, thus increasing the CE; however, its mechanism has remained controversial. In this work, we used a combination of LiN(SO2F)2 (LiFSI) and weakly solvating 1,2-diethoxyethane (DEE) as a model electrolyte to study how its liquid structure changes upon increasing salt concentration and how it is linked to the Li plating/stripping CE. Based on previous works, we focused on the Li electrode potential (ELi with reference to the redox potential of ferrocene) and solid-electrolyte-interphase (SEI) formation. Although ELi shows a different trend with DEE compared to conventional 1,2-dimethoxyethane (DME), which is accounted for by different ion-pair states of Li+ and FSI-, the ELi-CE plots overlap for both electrolytes, suggesting that ELi is one of the dominant factors of the CE. On the other hand, the extensive ion pairing results in the upward shift of the FSI- reduction potential, as demonstrated both experimentally and theoretically, which promotes the FSI--derived inorganic SEI. Both ELi and SEI contribute to increasing the Li plating/stripping CE.

Atomic and Electronic Structure in MgO-SiO2.

Understanding disordered structure is difficult due to insufficient information in experimental data. Here, we overcome this issue by using a combination of diffraction and simulation to investigate oxygen packing and network topology in glassy (g-) and liquid (l-) MgO-SiO2 based on a comparison with the crystalline topology. We find that packing of oxygen atoms in Mg2SiO4 is larger than that in MgSiO3, and that of the glasses is larger than that of the liquids. Moreover, topological analysis suggests that topological similarity between crystalline (c)- and g-(l-) Mg2SiO4 is the signature of low glass-forming ability (GFA), and high GFA g-(l-) MgSiO3 shows a unique glass topology, which is different from c-MgSiO3. We also find that the lowest unoccupied molecular orbital (LUMO) is a free electron-like state at a void site of magnesium atom arising from decreased oxygen coordination, which is far away from crystalline oxides in which LUMO is occupied by oxygen's 3s orbital state in g- and l-MgO-SiO2, suggesting that electronic structure does not play an important role to determine GFA. We finally concluded the GFA of MgO-SiO2 binary is dominated by the atomic structure in terms of network topology.

Thermodynamic aspect of sulfur, polysulfide anion and lithium polysulfide: plausible reaction path during discharge of lithium-sulfur battery.

The elucidation of elemental redox reactions of sulfur is important for improving the performance of lithium-sulfur batteries. The energies of stable structures of Sn, Sn˙-, Sn2-, [LiSn]- and Li2Sn (n = 1-8) were calculated at the CCSD(T)/cc-pVTZ//MP3/cc-pVDZ level. The heats of reduction reactions of S8 and Li2Sn with Li in the solid phase were estimated from the calculated energies and sublimation energies. The estimated heats of the redox reactions show that there are several redox reactions with nearly identical heats of reaction, suggesting that several reactions can proceed simultaneously at the same discharge voltage, although the discharging process was often explained by stepwise reduction reactions. The reduction reaction for the formation of Li2Sn (n = 2-6 and 8) from S8 normalized as a one electron reaction is more exothermic than that for the formation of Li2S directly from S8, while the reduction reactions for the formation of Li2S from Li2Sn are slightly less exothermic than that for the formation of Li2S directly from S8. If the reduction reactions with large exotherm occur first, these results suggest that the reduction reactions forming Li2Sn (n = 2-6 and 8) from S8 occur first, then Li2S is formed, and therefore, a two-step discharge-curve is observed.

Determining Factor on the Polarization Behavior of Magnesium Deposition for Magnesium Battery Anode.

To clarify the origin of the polarization of magnesium deposition/dissolution reactions, we combined electrochemical measurement, operando soft X-ray absorption spectroscopy (operando SXAS), Raman, and density functional theory (DFT) techniques to three different electrolytes: magnesium bis(trifluoromethanesulfonyl)amide (Mg(TFSA)2)/triglyme, magnesium borohydride (Mg(BH4)2)/tetrahydrofuran (THF), and Mg(TFSA)2/2-methyltetrahydrofuran (2-MeTHF). Cyclic voltammetry revealed that magnesium deposition/dissolution reactions occur in Mg(TFSA)2/triglyme and Mg(BH4)2/THF, while the reactions do not occur in Mg(TFSA)2/2-MeTHF. Raman spectroscopy shows that the [TFSA]- in the Mg(TFSA)2/triglyme electrolyte largely does not coordinate to the magnesium ions, while all of the [TFSA]- in Mg(TFSA)2/2-MeTHF and [BH4]- in Mg(BH4)2/THF coordinate to the magnesium ions. In operando SXAS measurements, the intermediate, such as the Mg+ ion, was not observed at potentials above the magnesium deposition potential, and the local structure distortion around the magnesium ions increases in all of the electrolytes at the magnesium electrode|electrolyte interface during the cathodic polarization. Our DFT calculation and X-ray photoelectron spectroscopy results indicate that the [TFSA]-, strongly bound to the magnesium ion in the Mg(TFSA)2/2-MeTHF electrolyte, undergoes reduction decomposition easily, instead of deposition of magnesium metal, which makes the electrolyte inactive electrochemically. In the Mg(BH4)2/THF electrolyte, because the [BH4]- coordinated to the magnesium ions is stable even under the potential of the magnesium deposition, the magnesium deposition is not inhibited by the decomposition of [BH4]-. Conversely, because [TFSA]- is weakly bound to the magnesium ion in Mg(TFSA)2/triglyme, the reduction decomposition occurs relatively slowly, which allows the magnesium deposition in the electrolyte.

Li-ion transport at the interface between a graphite anode and Li2CO3 solid electrolyte interphase: ab initio molecular dynamics study.

Understanding and the control of Li-ion (Li+) transport across the interface between the anode and solid electrolyte interphase (SEI) film or electrolyte is a key issue in battery electrochemistry and interface science. In this study, we investigated the structural, electronic and free energy properties of Li+ migration between a Li-intercalated graphite anode LiCx and Li2CO3 SEI film, by using ab initio molecular dynamics and free energy calculations. We compared three types of graphite edges: H-, OH- and mixed (H, OH, COOH)-terminations, and three cases of transferred Li-ions: Li+ constructing the SEI, excess Li+ and excess Li0 (excess Li+ + e- in anode). After validation of our calculations with Li2CO3 and LiCx bulk systems, we sampled the interfacial structures under thermodynamic equilibrium and demonstrated that the OH- and mixed-terminations had larger binding energies. The calculated free energy profiles of Li+ intercalation from the Li2CO3 SEI to LiC24 showed barriers larger than 1.2 eV irrespective of the terminations and Li+ cases. We also clarified that the charges of Li ions did not change much upon the intercalation. Based on these results and the calculated Li chemical potential, we constructed the probable free energy profile of Li+ between the anode and cathode under charging and discharging. This profile model suggest a possible electric field approximation for the charging stage, and the resultant free energy profiles with such fields gave a ca. 0.5 eV barrier under charging, which was consistent with the experimental values. The present picture will give a crucial insight into Li-ion transport at the battery interfaces.

Machine learning prediction of coordination energies for alkali group elements in battery electrolyte solvents.

We combined a data science-driven method with quantum chemistry calculations, and applied it to the battery electrolyte problem. We performed quantum chemistry calculations on the coordination energy (Ecoord) of five alkali metal ions (Li, Na, K, Rb, and Cs) to electrolyte solvent, which is intimately related to ion transfer at the electrolyte/electrode interface. Three regression methods, namely, multiple linear regression (MLR), least absolute shrinkage and selection operator (LASSO), and exhaustive search with linear regression (ES-LiR), were employed to find the relationship between Ecoord and descriptors. Descriptors include both ion and solvent properties, such as the radius of metal ions or the atomic charge of solvent molecules. Our results clearly indicate that the ionic radius and atomic charge of the oxygen atom that is connected to the metal ion are the most important descriptors. Good prediction accuracy for Ecoord of 0.127 eV was obtained using ES-LiR, meaning that we can predict Ecoord for any alkali ion without performing quantum chemistry calculations for ion-solvent pairs. Further improvement in the prediction accuracy was made by applying the exhaustive search with Gaussian process, which yields 0.016 eV for the prediction accuracy of Ecoord.

Prediction and optimization of epoxy adhesive strength from a small dataset through active learning.

Machine learning is emerging as a powerful tool for the discovery of novel high-performance functional materials. However, experimental datasets in the polymer-science field are typically limited and they are expensive to build. Their size (< 100 samples) limits the development of chemical intuition from experimentalists, as it constrains the use of machine-learning algorithms for extracting relevant information. We tackle this issue to predict and optimize adhesive materials by combining laboratory experimental design, an active learning pipeline and Bayesian optimization. We start from an initial dataset of 32 adhesive samples that were prepared from various molecular-weight bisphenol A-based epoxy resins and polyetheramine curing agents, mixing ratios and curing temperatures, and our data-driven method allows us to propose an optimal preparation of an adhesive material with a very high adhesive joint strength measured at 35.8 ± 1.1 MPa after three active learning cycles (five proposed preparations per cycle). A Gradient boosting machine learning model was used for the successive prediction of the adhesive joint strength in the active learning pipeline, and the model achieved a respectable accuracy with a coefficient of determination, root mean square error and mean absolute error of 0.85, 4.0 MPa and 3.0 MPa, respectively. This study demonstrates the important impact of active learning to accelerate the design and development of tailored highly functional materials from very small datasets.

Combined Theoretical and Experimental Studies of Sodium Battery Materials.

Owing to developments in theoretical chemistry and computer power, the combination of calculations and experiments is now standard practice in understanding and developing new materials for battery systems. Here, we briefly review our recent combined studies based on density functional theory and molecular dynamics calculations for electrode and electrolyte materials for sodium-ion batteries. These findings represent case studies of successful combinations of experimental and theoretical methods.

Liquid electrolyte informatics using an exhaustive search with linear regression.

Exploring new liquid electrolyte materials is a fundamental target for developing new high-performance lithium-ion batteries. In contrast to solid materials, disordered liquid solution properties have been less studied by data-driven information techniques. Here, we examined the estimation accuracy and efficiency of three information techniques, multiple linear regression (MLR), least absolute shrinkage and selection operator (LASSO), and exhaustive search with linear regression (ES-LiR), by using coordination energy and melting point as test liquid properties. We then confirmed that ES-LiR gives the most accurate estimation among the techniques. We also found that ES-LiR can provide the relationship between the "prediction accuracy" and "calculation cost" of the properties via a weight diagram of descriptors. This technique makes it possible to choose the balance of the "accuracy" and "cost" when the search of a huge amount of new materials was carried out.

First-Principles Study of Electron Injection and Defects at the TiO2/CH3NH3PbI3 Interface of Perovskite Solar Cells.

We investigated electron injection rates and vacancy defect properties by performing first-principles calculations on the interface of an anatase-TiO2(001) and a tetragonal CH3NH3PbI3(110) (MAPbI3(110)). We found that the coupling matrix element between the lowest unoccupied molecular orbital of MAPbI3 and the TiO2 conduction band (CB) minimum is negligibly small, the indication being that electron-injection times for low-energy excited states are quite long (more than several tens of picoseconds). We also found that higher-lying CB states coupled more strongly; injection was expected to take place on a femtosecond time scale. Furthermore, we found that vacancy defects in the TiO2 layer produced undesired defect levels that caused hole traps and recombination centers. Whereas most of the vacancy defects in the MAPbI3 layer produced no additional states in the MAPbI3 gap, a Pb vacancy (VPb) at the interface created an energy level below the MAPbI3 CB edge and had a lower energy of formation than the VPb defect in bulk because of the interaction with the TiO2 surface.

Unusual Passivation Ability of Superconcentrated Electrolytes toward Hard Carbon Negative Electrodes in Sodium-Ion Batteries.

The passivation of negative electrodes is key to achieving prolonged charge-discharge cycling with Na-ion batteries. Here, we report the unusual passivation ability of superconcentrated Na-salt electrolytes. For example, a 50 mol % sodium bis(fluorosulfonyl)amide (NaFSA)/succinonitrile (SN) electrolyte enables highly reversible Na+ insertion into a hard carbon negative electrode without any electrolyte additive, functional binder, or electrode pretreatment. Importantly, an anion-derived passivation film is formed via preferential reduction of the anion upon charging, which can effectively suppress further electrolyte reduction. As a structural characteristic of the electrolyte, most anions are coordinated to multiple Na+ cations at high concentration, which shifts the lowest unoccupied molecular orbitals of the anions downward, resulting in preferential anion reduction. The present work provides a new understanding of the passivation mechanism with respect to the coordination state of the anion.

Cation Mixing Properties toward Co Diffusion at the LiCoO2 Cathode/Sulfide Electrolyte Interface in a Solid-State Battery.

All-solid-state Li-ion batteries (ASS-LIBs) are expected to be the next-generation battery, however, their large interfacial resistance hinders their widespread application. To understand and resolve the possible causes of this resistance, we examined mutual diffusion properties of the cation elements at LiCoO2 (LCO) cathode/ß-Li3PS4 (LPS) solid electrolyte interface as a representative system as well as the effect of a LiNbO3 buffer layer by first-principles calculations. Evaluating energies of exchanging ions between the cathode and the electrolyte, we found that the mixing of Co and P is energetically preferable to the unmixed states at the LCO/LPS interface. We also demonstrated that the interposition of the buffer layer suppresses such mixing because the exchange of Co and Nb is energetically unfavorable. Detailed analyses of the defect levels and the exchange energies by using the individual bulk crystals as well as the interfaces suggest that the lower interfacial states in the energy gap can make a major contribution to the stabilization of the Co ↔ P exchange, although the anion bonding preference of Co and P as well as the electrostatic interactions may have effects as well. Finally, the Co ↔ P exchanges induce interfacial Li sites with low chemical potentials, which enhance the growth of the Li depletion layer. These atomistic understandings can be meaningful for the development of ASS-LIBs with smaller interfacial resistances.

The Solvation Structure of Lithium Ions in an Ether Based Electrolyte Solution from First-Principles Molecular Dynamics.

The solvation and desolvation of the Li ion play a crucial role in the electrolytes of Li based secondary batteries, and their understanding at the microscopic level is of great importance. Oligoether (glyme) based electrolytes have attracted much attention as electrolytes used in Li based secondary batteries, such as Li-ion, Li-S, and Li-O2 batteries. However, the solvation structure of the Li ion in glyme based electrolytes has not been fully clarified yet. We present a computational study on the solvation structure of lithium ions in the mixture of triglyme and lithium bis(trifluoromethylsulfonyl)-amide (LiTFSA) by means of molecular orbital and molecular dynamics calculations based on density functional theory. We found that, in the electrolyte solution composed of the equimolar mixture of triglyme and LiTFSA, lithium ions are solvated mainly by crown-ether-like curled triglyme molecules and in direct contact with an TFSA anion. We also found the aggregate formed with Li ion and TFSA anions and/or triglyme molecule(s) is equally stable, which has not been reported in the previous classical molecular dynamics simulations, suggesting that in reality a small fraction of Li ions form aggregates and they might have a significant impact on the Li ion transport. Our results demonstrate the importance of performing electronic structure based molecular dynamics of electrolyte solution to clarify the detailed solvation structure of the Li ion.

Superconcentrated electrolytes for a high-voltage lithium-ion battery.

Finding a viable electrolyte for next-generation 5 V-class lithium-ion batteries is of primary importance. A long-standing obstacle has been metal-ion dissolution at high voltages. The LiPF6 salt in conventional electrolytes is chemically unstable, which accelerates transition metal dissolution of the electrode material, yet beneficially suppresses oxidative dissolution of the aluminium current collector; replacing LiPF6 with more stable lithium salts may diminish transition metal dissolution but unfortunately encounters severe aluminium oxidation. Here we report an electrolyte design that can solve this dilemma. By mixing a stable lithium salt LiN(SO2F)2 with dimethyl carbonate solvent at extremely high concentrations, we obtain an unusual liquid showing a three-dimensional network of anions and solvent molecules that coordinate strongly to Li(+) ions. This simple formulation of superconcentrated LiN(SO2F)2/dimethyl carbonate electrolyte inhibits the dissolution of both aluminium and transition metal at around 5 V, and realizes a high-voltage LiNi0.5Mn1.5O4/graphite battery that exhibits excellent cycling durability, high rate capability and enhanced safety.

Decomposition of the fluoroethylene carbonate additive and the glue effect of lithium fluoride products for the solid electrolyte interphase: an ab initio study.

Additives in the electrolyte solution of lithium-ion batteries (LIBs) have a large impact on the performance of the solid electrolyte interphase (SEI) that forms on the anode and is a key to the stability and durability of LIBs. We theoretically investigated effects of fluoroethylene carbonate (FEC), a representative additive, that has recently attracted considerable attention for the enhancement of cycling stability of silicon electrodes and the improvement of reversibility of sodium-ion batteries. First, we intensively examined the reductive decompositions by ring-opening, hydrogen fluoride (HF) elimination to form a vinylene carbonate (VC) additive and intermolecular chemical reactions of FEC in the ethylene carbonate (EC) electrolyte, by using density functional theory (DFT) based molecular dynamics and the blue-moon ensemble technique for the free energy profile. The results show that the most plausible product of the FEC reductive decomposition is lithium fluoride (LiF), and that the reactivity of FEC to anion radicals is found to be inert compared to the VC additive. We also investigated the effects of the generated LiF on the SEI by using two model systems; (1) LiF molecules distributed in a model aggregate of organic SEI film components (SFCs) and (2) a LiF aggregate interfaced with the SFC aggregate. DFT calculations of the former system show that F atoms form strong bindings with the Li atoms of multiple organic SFC molecules and play as a joint connecting them. In the latter interface system, the LiF aggregate adsorbs the organic SFCs through the F-Li bindings. These results suggest that LiF moieties play the role of glue in the organic SFC within the SEI film. We also examined the interface structure between a LiF aggregate and a lithiated silicon anode, and found that they are strongly bound. This strong binding is likely to be related to the effectiveness of the FEC additive in the electrolyte for the silicon anode.

Sodium-Ion Intercalation Mechanism in MXene Nanosheets.

MXene, a family of layered compounds consisting of nanosheets, is emerging as an electrode material for various electrochemical energy storage devices including supercapacitors, lithium-ion batteries, and sodium-ion batteries. However, the mechanism of its electrochemical reaction is not yet fully understood. Herein, using solid-state (23)Na magic angle spinning NMR and density functional theory calculation, we reveal that MXene Ti3C2Tx in a nonaqueous Na(+) electrolyte exhibits reversible Na(+) intercalation/deintercalation into the interlayer space. Detailed analyses demonstrate that Ti3C2Tx undergoes expansion of the interlayer distance during the first sodiation, whereby desolvated Na(+) is intercalated/deintercalated reversibly. The interlayer distance is maintained during the whole sodiation/desodiation process due to the pillaring effect of trapped Na(+) and the swelling effect of penetrated solvent molecules between the Ti3C2Tx sheets. Since Na(+) intercalation/deintercalation during the electrochemical reaction is not accompanied by any substantial structural change, Ti3C2Tx shows good capacity retention over 100 cycles as well as excellent rate capability.

Surface Properties of CH3NH3PbI3 for Perovskite Solar Cells.

Perovskite solar cells (PSCs) have attracted considerable interest because of their high potential for solar energy conversion. Power conversion efficiencies of the PSCs have rapidly increased from 3.8 to over 20% only in the past few years. PSCs have several similarities to dye-sensitized solar cells in their device compositions; mesoporous TiO2 (mp-TiO2) is sensitized by light-absorbing components and placed into a medium containing hole transporting materials (HTMs). On the other hand, the perovskite materials for the light-harvesting, for example, CH3NH3PbI3 (MAPbI3), have a greater advantage for the photovoltaic applications; extremely long photocarrier diffusion lengths (over 1 µm) enable carrier transports without singnificant loss. In this respect, the surface states, that can be possible recombination centers, are also of great importance. Availability of solution processes is another important aspect in terms of low cost fabrication of PSCs. Two-step methods, where PbI2 is first introduced from solution onto a mp-TiO2 film and subsequently transformed into the MAPbI3 by the exposition of a solution containing MAI, suggest that use of such a high PbI2 concentration is crucial to obtain higher performance. The experiments also indicate that the PbI2-rich growth condition modifies TiO2/ or HTM/MAPbI3 interfaces in such a way that the photocarrier transport is improved. Thus, the characteristics of surfaces and interfaces play key roles in the high efficiencies of the PSCs. In this Account, we focus on the structural stability and electronic states of the representative (110), (001), (100), and (101) surfaces of tetragonal MAPbI3, which can be regarded as reasonable model HTM/MAPbI3 interfaces, by use of first-principles calculations. By examining various types of PbIx polyhedron terminations, we found that there are two major phases on all of the four surface facets. They can be classified as vacant- and flat-type terminations, and the former is more stable than the latter under thermodynamically equilibrium conditions. More interestingly, both terminations can coexist especially on the more probable (110) and (001) surfaces. Electronic states, that is, projected density of states, of the stable-vacant and PbI2-rich-flat terminations on the two surfaces are almost the same as that in bulk MAPbI3. These surfaces can contribute to the long carrier lifetime actually observed for the PSCs because they have no midgap surface states. Furthermore, the shallow surface states on the (110) and (001) flat terminations can be efficient intermediates for hole transport to HTMs. Consequently, the formation of the flat terminations under the PbI2-rich condition will be beneficial for the improvement of PSC performance.

First-Principles Study of Ion Diffusion in Perovskite Solar Cell Sensitizers.

Hysteresis in current-voltage curves has been an important issue for conversion efficiency evaluation and development of perovskite solar cells (PSCs). In this study, we explored the ion diffusion effects in tetragonal CH3NH3PbI3 (MAPbI3) and trigonal (NH2)2CHPbI3 (FAPbI3) by first-principles calculations. The calculated activation energies of the anionic and cationic vacancy migrations clearly show that I(-) anions in both MAPbI3 and FAPbI3 can easily diffuse with low barriers of ca. 0.45 eV, comparable to that observed in ion-conducting materials. More interestingly, typical MA(+) cations and larger FA(+) cations both have rather low barriers as well, indicating that the cation molecules can migrate in the perovskite sensitizers when a bias voltage is applied. These results can explain the ion displacement scenario recently proposed by experiments. With the dilute diffusion theory, we discuss that smaller vacancy concentrations (higher crystallinity) and replacement of MA(+) with larger cation molecules will be essential for suppressing hysteresis as well as preventing aging behavior of PSC photosensitizers.

A double-QM/MM method for investigating donor-acceptor electron-transfer reactions in solution.

We developed a double-quantum mechanical/molecular mechanical (d-QM/MM) method for investigation of full outer-sphere electron transfer (ET) processes between a donor and an acceptor (DA) in condensed matter. In the d-QM/MM method, which employs the novel concept of multiple QM regions, one can easily specify the number of electrons, spin states and appropriate exchange-correlation treatment in each QM region, which is especially important in the cases of ET involving transition metal sites. We investigated Fe(2+/3+) self-exchange and Fe(3+) + Ru(2+) → Fe(2+) + Ru(3+) in aqueous solution as model reactions, and demonstrated that the d-QM/MM method gives reasonable accuracy for the redox potential, reorganization free energy and electronic coupling. In particular, the DA distance dependencies of those quantities are clearly shown at the density functional theory hybrid functional level. The present d-QM/MM method allows us to explore the intermediate DA distance region, important for long-range ET phenomena observed in electrochemistry (on the solid-liquid interfaces) and biochemistry, which cannot be dealt by the half-reaction scheme with the conventional QM/MM.

Arylpyrrole oligomers as tunable anion receptors.

A novel type of arylpyrrole oligomer possessing an appropriate electropositive cavity has been designed, prepared and analysed for use as readily accessible receptors for negatively charged guests. Affinities of the receptors for various anions were determined by UV/Vis titration experiments and in depth insights into the host-guest interactions were gained by performing (1)H NMR titration experiments and X-ray crystallographic structure analyses. Experimentally determined association constants were correlated with the calculated maximum electrostatic potentials of the electropositive cavities of the receptors, allowing estimation of the strengths of host-guest associations in similar compounds. The joint contribution of aryl C-H and pyrrole N-H hydrogens was shown to be key to a strong guest association, resulting in the arylpyrrole oligomers being efficient anion receptors.