EwaldSolidSolution: A High-Throughput Application to Quickly Sample Stable Site Arrangements for Ionic Solid Solutions.

Data-driven materials design of ionic solid solutions often requires sampling (meta)stable site arrangements among the massive number of possibilities, which has been hampered by the lack of relevant methods. Herein, we develop a quick high-throughput sampling application for site arrangements of ionic solid solutions. Given the Ewald Coulombic energies for an initial site arrangement, EwaldSolidSolution updates the modified parts of the energy with varying sites only, which can be exhaustively estimated by using massively parallel processing. Given two representative examples of solid electrolytes, Li10GeP2S12 and Na3Zr2Si2PO12, EwaldSolidSolution successfully calculates the Ewald Coulombic energies of 211,266,225 (235,702,467) site arrangements for Li10GeP2S12 (Na3Zr2Si2PO12) with 216 (160) ion sites per unit cell in 1223.2 (1187.9) seconds: 0.0057898 (0.0050397) milliseconds per site arrangement. The computational cost is enormously saved in comparison with an existing application, which estimates the energy of a site arrangement on the second timescale. The positive correlations between the Ewald Coulombic energies and those estimated by density functional theory calculations show that (meta)stable samples are easily revealed by our computationally inexpensive algorithm. We also reveal that the different-valence nearest-neighbor pairs are distinctively formed in the low-energy site arrangements. EwaldSolidSolution will boost the materials design of ionic solid solutions by attracting broad interest.

On-Surface Synthesis of Porphyrin-Complex Multi-Block Co-Oligomers by Defluorinative Coupling.

On-surface chemical reaction has become a very powerful technique to synthesize nanostructures by linking small molecules in the bottom-up approach. Given the fact that most reactants are simultaneously activated at certain temperatures, a sequential reaction in a controlled way has remained challenging. Here, we present an on-surface synthesis of multi-block co-oligomers from trifluoromethyl (CF3 ) substituted porphyrin metal complexes. The oligomerization on Au(111) is demonstrated with a combination of bond-resolved scanning probe microscopy and density functional theory (DFT) calculations. Even after the first oligomerization of single monomer unit, the termini of the oligomer keep the CF3 group, which can be used as a reactant for further coupling in a sequential order. Consequently, copper, cobalt, and palladium complexes of bisanthracene-fused porphyrin oligomers were linked with each other in a designed order.

Boron position-dependent surface reconstruction and electronic states of boron-doped diamond(111) surfaces: an ab initio study.

Boron-doped diamond (BDD) has attracted much attention in semi-/superconductor physics and electrochemistry, where the surface structures and electronic states play crucial roles. Herein, we systematically examine the structural and electronic properties of the unterminated and H-terminated diamond(111) surfaces by using density functional theory calculations, and the effect of the boron position on them. The surface energy increases compared to that of the undoped case when the boron is located at a deeper position in the diamond bulk, which indicates that boron near the surface can facilitate the surface stability of the BDD in addition to the H-termination. Moreover, the surface energy and projected density of state analyses suggest that the boron can enhance the graphitization of the pristine (ideal) unterminated (111) surface thanks to the alternative sp2-sp3 arrangement on that surface. Finally, we found that surface electronic states depend on the boron's position, i.e., the Fermi energy (EF) is located around the mid-gap position when the boron lies near the surface, instead of showing a p-type semiconductor behavior where the EF lies closer to the valence band maximum.

MgO-Template Synthesis of Extremely High Capacity Hard Carbon for Na-Ion Battery.

Extremely high capacity hard carbon for Na-ion battery, delivering 478 mAh g-1 , is successfully synthesized by heating a freeze-dried mixture of magnesium gluconate and glucose by a MgO-template technique. Influences of synthetic conditions and nano-structures on electrochemical Na storage properties in the hard carbon are systematically studied to maximize the reversible capacity. Nano-sized MgO particles are formed in a carbon matrix prepared by pre-treatment of the mixture at 600 °C. Through acid leaching of MgO and carbonization at 1500 °C, resultant hard carbon demonstrates an extraordinarily large reversible capacity of 478 mAh g-1 with a high Coulombic efficiency of 88 % at the first cycle.

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.

In Situ Spectroscopic Study on the Surface Hydroxylation of Diamond Electrodes.

Carbon-based materials are regarded as an environmentally benign alternative to the conventional metal electrode used in electrochemistry from the viewpoint of sustainable chemistry. Among various carbon electrode materials, boron-doped diamond (BDD) exhibits superior electrochemical properties. However, it is still uncertain how surface chemical species of BDD influence the electrochemical performance, because of the difficulty in characterizing the surface species. Here, we have developed in situ spectroscopic measurement systems on BDD electrodes, i.e., in situ attenuated total reflection infrared spectroscopy (ATR-IR) and electrochemical X-ray photoelectron spectroscopy (EC-XPS). ATR-IR studies at a controlled electrode potential confirmed selective surface hydroxylation. EC-XPS studies confirmed deprotonation of C-OH groups at the BDD/electrolyte interface. These findings should be important not only for better understanding of BDD's fundamentals but also for a variety of applications.

Reaction energy benchmarks of hydrocarbon combustion by Gaussian basis and plane wave basis approaches.

The reaction energies of 275 elementary reactions from the hydrocarbon combustion model GRI-Mech 3.0 were evaluated by electronic structure calculations using both localized Gaussian basis and plane wave basis sets. In the Gaussian basis calculations, the d-polarization function on C, N, and O elements reduces the mean absolute deviation (MAD) from the experimental value by 53%, a significant improvement in computational accuracy. In the plane wave basis calculation using different exchange-correlation (XC) functionals, the MAD values were 0.316-0.426 eV when non-hybrid type XC functionals such as RPBE, PBE, PW91, revPBE, and PBEsol were used. On the other hand, hybrid functionals like B3LYP and HSE06 reduced the MAD values significantly down to 0.182 and 0.233 eV, respectively. The B3LYP results have 49% less MAD compared to the PBE results. These demonstrated the strong advantage of the hybrid functional for calculating gas-phase reaction energies. The present comprehensive benchmarks will be crucial for future microkinetics as well as machine learning studies on the catalytic reactions. © 2019 Wiley Periodicals, Inc.

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.

On principal features of organic electrolyte molecules in lithium ion battery performance.

The lithium (Li) complexes of organic electrolyte solvents are theoretically investigated using the long-range correction for density functional theory in order to figure out the cause for the high performance of cyclic carbonate electrolytes in lithium ion batteries (LIBs). Calculations of the Li complexes with ethylene carbonate solvent molecules prove that ten ligand molecules should be incorporated to obtain near-degenerate four- and five-coordination optimum structures and dramatically improved orbital energies. The geometry optimizations of the Li complexes with thirteen types of organic solvent molecules give four-coordination neutral and five-coordination cation complexes for many solvent molecules. The five-coordination Li complexes are considered to use Berry pseudorotation to approach the electrodes from the Li atom. The calculated Koopmans, vertical and adiabatic ionization potentials and electron affinities show that near-degeneracy and structural deformation effects play significant roles in the electronic states of the Li complexes. Mulliken charge and dipole moment analyses indicate that the Li complexes of cyclic carbonates construct a deep electric double layer near electrodes due to the electron-donating ability of the ligand molecules. Molecular orbital analyses also explain that the Li complexes of cyclic carbonates easily construct a solid electrolyte interface, which contributes to Li ion conductance, by localizing the accepted electron to one ligand molecule. In conclusion, the Li complexes of cyclic carbonates have three main features: preference of five-coordination structures, high electron-donating ability of ligand molecules, and localization of the accepted electron to one ligand molecule.

Electrochemical properties of fluorinated boron-doped diamond electrodes via fluorine-containing plasma treatment.

Since a boron-doped diamond (BDD) exhibits excellent electrode properties such as wide potential window, low back-ground current, and high physical and chemical durability, it has been studied as an electrode material for various electrochemical applications. The electrochemical behavior of BDD depends on the surface termination, which can be easily converted by chemical reactions. Fluorine termination has attracted interest because it exhibits unique surface properties such as high hydrophobicity and a low coefficient of friction, and the electrochemical properties also drastically change. However, so far, it has not been elucidated why fluorinated BDD exhibits specific electrochemical properties. In this article, fluorine-terminated BDD was fabricated by a fluorine-containing plasma treatment, and the electrochemical properties were systematically investigated. Together with experiments, we have calculated the interfacial structures and electronic states of hydrogenated, oxygenated, and fluorinated BDD electrodes. As a result, fluorinated BDD showed lower electrochemical reactivity than hydrogenated and oxygenated BDD. Especially, electron transfer between anionic compound and fluorinated BDD was significantly suppressed. Considered together with theoretical calculations, this reactivity could be attributed to the larger interfacial band bending in fluorinated BDD and electrostatic interactions between BDD and redox species.

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.

Modifications in coordination structure of Mg[TFSA]2-based supporting salts for high-voltage magnesium rechargeable batteries.

To achieve a sustainable-energy society in the future, next-generation highly efficient energy storage technologies, particularly those based on multivalent metal negative electrodes, are urgently required to be developed. Magnesium rechargeable batteries (MRBs) are promising options owing to the many advantageous chemical and electrochemical properties of magnesium. However, the substantially low working voltage of sulfur-based positive electrodes may hinder MRBs in becoming alternatives to current Li-ion batteries. We proposed halide-free noncorrosive ionic liquid-based electrolytes incorporating Mg[TFSA]2 for high-voltage MRB applications. Upon the complexation of Mg[TFSA]2 with tetraglyme (G4) and strict control of the liquid states, the electrolytes achieved excellent anodic stability up to 4.1 V vs. Mg2+/Mg even at 100 °C. The modest electrochemical activities for magnesium deposition/dissolution in the [Mg(G4)][TFSA]2/ionic liquid electrolyte can be improved by certain modifications to the coordination state of [TFSA]-. Dialkyl sulfone was found to be effective in changing the coordination state of [TFSA]- from associated to isolated (free). This coordination change successfully promoted magnesium deposition/dissolution reactions, particularly in the coexistence of ether ligand. By contrast, the coordination of Mg2+ by strongly donating agents such as dimethyl sulfoxide and alkylimidazole led to the complexes inactive electrochemically, suggesting that interaction between Mg2+ and coordination agents predominates the fundamental electrochemical activity. We also demonstrated that an enhancement in the electrochemical activity of electrolytes contributed to improvements in the cycling ability of magnesium batteries with 2.5 V-class MgMn2O4 positive electrodes.

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.

A general representation scheme for crystalline solids based on Voronoi-tessellation real feature values and atomic property data.

Increasing attention has been paid to materials informatics approaches that promise efficient and fast discovery and optimization of functional inorganic materials. Technical breakthrough is urgently requested to advance this field and efforts have been made in the development of materials descriptors to encode or represent characteristics of crystalline solids, such as chemical composition, crystal structure, electronic structure, etc. We propose a general representation scheme for crystalline solids that lifts restrictions on atom ordering, cell periodicity, and system cell size based on structural descriptors of directly binned Voronoi-tessellation real feature values and atomic/chemical descriptors based on the electronegativity of elements in the crystal. Comparison was made vs. radial distribution function (RDF) feature vector, in terms of predictive accuracy on density functional theory (DFT) material properties: cohesive energy (CE), density (d), electronic band gap (BG), and decomposition energy (Ed). It was confirmed that the proposed feature vector from Voronoi real value binning generally outperforms the RDF-based one for the prediction of aforementioned properties. Together with electronegativity-based features, Voronoi-tessellation features from a given crystal structure that are derived from second-nearest neighbor information contribute significantly towards prediction.

Surface Hydrogenation of Boron-Doped Diamond Electrodes by Cathodic Reduction.

Boron-doped diamond (BDD) has attracted much attention as a promising electrode material especially for electrochemical sensing systems, because it has excellent properties such as a wide potential window and low background current. It is known that the electrochemical properties of BDD electrodes are very sensitive to the surface termination such as to whether it is hydrogen- or oxygen-terminated. Pretreating BDD electrodes by cathodic reduction (CR) to hydrogenate the surface has been widely used to achieve high sensitivity. However, little is known about the effects of the CR treatment conditions on surface hydrogenation. In this Article, we report on a systematic study of CR treatments that can achieve effective surface hydrogenation. As a result, we found that the surface hydrogenation could be improved by applying a more negative potential in a lower pH solution. This is because hydrogen atoms generated from protons in the CR treatment contribute to the surface hydrogenation. After CR treatments, BDD surface could be hydrogenated not completely but sufficiently to achieve high sensitivity for electrochemical sensing. In addition, we confirmed that hydrogenation with high repeatability could be achieved.

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.

Catalytic Proton Dynamics at the Water/Solid Interface of Ceria-Supported Pt Clusters.

Wet conditions in heterogeneous catalysis can substantially improve the rate of surface reactions by assisting the diffusion of reaction intermediates between surface reaction sites. The atomistic mechanisms underpinning this accelerated mass transfer are, however, concealed by the complexity of the dynamic water/solid interface. Here we employ ab initio molecular dynamics simulations to disclose the fast diffusion of protons and hydroxide species along the interface between water and ceria, a catalytically important, highly reducible oxide. Up to 20% of the interfacial water molecules are shown to dissociate at room temperature via proton transfer to surface O atoms, leading to partial surface hydroxylation and to a local increase of hydroxide species in the surface solvation layer. A water-mediated Grotthus-like mechanism is shown to activate the fast and long-range proton diffusion at the water/oxide interface. We demonstrate the catalytic importance of this dynamic process for water dissociation at ceria-supported Pt nanoparticles, where the solvent accelerates the spillover of ad-species between oxide and metal sites.

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.

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 method to calculate redox potentials relative to the normal hydrogen electrode in nonaqueous solution by using density functional theory-based molecular dynamics.

We demonstrate the redox potential calculation relative to the normal hydrogen electrode (NHE) in nonaqueous solution using a density functional theory-based molecular dynamics (DFT-MD) simulation. The calculation of the NHE in nonaqueous solution consists of two processes: the first step is the equilibrated simulation for a proton in nonaqueous solution to determine the space for inserting a proton in solution, and the second step is the thermodynamic integration method to calculate the solvation energy of the proton in the nonaqueous solution. In this work, we apply the method for a cation and an anion, i.e., copper(ii)/copper(i) and iodine/iodide in acetonitrile solution, and show that the errors in the calculated redox potential from experiments are within 0.21 V.