Path ahead: Tackling the Challenge of Computationally Estimating Lithium Diffusion in Cathode Materials.

In the roadmap toward designing new and improved materials for Lithium ion batteries, the ability to estimate the diffusion coefficient of Li atoms in electrodes, and eventually solid-state electrolytes, is key. Nevertheless, as of today, accurate prediction through computational tools remains challenging. Its experimental measurement does not appear to be much easier. In this work, we devise a computational protocol for the determination of the Li-migration energy barrier and diffusion coefficient, focusing on a common cathode material such as LiNiO2, which represents a prototype of the widely adopted NMC (LiNi1-x-y Mn x Co y O2) class of materials. Different methodologies are exploited, combining ab initio metadynamics, path sampling, and density functional theory. Furthermore, we propose a novel, fast, and simple 1D approximation for the estimation of the effective frequency. The outlined computational protocol aims to be generally applicable to Lithium diffusion in other materials and components for batteries, including anodes and solid electrolytes.

A computational study of the negative LiIn modified anode and its interaction with ß-Li3PS4 solid-electrolyte for battery applications.

All-solid-state lithium batteries (ASSLBs) have sparked interest due to their far superior energy density compared to current commercial material, but the heightened reactivity of the negative Li electrode can compromise the long-term cyclability of the cell, calling for the introduction of passivating layers or alloy anodes. In this article, we aim to explain the outstanding stability of LiIn alloy-based anodes over extended cycling by comparing its bulk and interface properties to Li-metal. Using density functional theory, we conducted an in-depth analysis of the LiIn surfaces' formation and subsequent structural stability in interfaces with the solid electrolyte ß-Li3PS4. Several LiIn facets are shown to possess sufficient structural stability, with the (110) surface being the most stable. The stable interfaces established with the ß-Li3PS4(100) surface featured favorable adhesion energy, low strain energy, and little reconstruction. By comparing these interface properties with the bulk properties of Li-metal and LiIn, we highlighted the influence of the cohesion energy, Fermi energy level, and band position of the two materials in the long-term stability of their anodes under battery conditions.

Understanding Ionic Diffusion Mechanisms in Li2S Coatings for Solid-State Batteries: Development of a Tailored Reactive Force Field for Multiscale Simulations.

In order to investigate Li2S as a potential protective coating for lithium anode batteries using superionic electrolytes, we need to describe reactions and transport for systems at scales of >10,000 atoms for time scales beyond nanoseconds, which is most impractical for quantum mechanics (QM) calculations. To overcome this issue, here, we first report the development of the reactive analytical force field (ReaxFF) based on density functional theory (DFT) calculations on model systems at the PBE0/TZVP and M062X/TZVP levels. Then, we carry out reactive molecular dynamics simulations (RMD) for up to 20 ns to investigate the diffusion mechanisms in bulk Li2S as a function of vacancy density, determining the activation barrier for diffusion and conductivity. We show that RMD predictions for diffusion and conductivity are comparable to experiments, while results on model systems are consistent with and validated by short (10-100 ps) ab initio molecular dynamics (AIMD). This new ReaxFF for Li2S systems enables practical RMD on spatial scales of 10-100 nm (10,000 to 10 million atoms) for the time scales of 20 ns required to investigate predictively the interfaces between electrodes and electrolytes, electrodes and coatings, and coatings and electrolytes during the charging and discharging processes.

Computational Understanding of Delithiation, Overlithiation, and Transport Properties in Disordered Cubic Rock-Salt Type Li2TiS3.

Lithium-titanium-sulfur cathodes have gained attention because of their unique properties and have been studied for their application in lithium-ion batteries. They offer different advantages such as lower cost, higher safety, and higher energy density with respect to commonly adopted transition metal oxides. Moreover, this family of compounds is free from critical raw materials such as cobalt and nickel. For cathode materials, a crucial aspect is evaluating the evolution and behavior of the structure and properties during the cycling process, which means simulating the system under lithium extraction and insertion. Structural optimization, electronic band structures, density of states, and Raman spectra were simulated, looking for fingerprints and peculiar aspects related to the delithiation and overlithiation process. Lithium transport properties were also investigated through the nudged elastic band methodology. This allowed us to evaluate the diffusion coefficient of lithium, which is a crucial parameter for cathode performance evaluation.

Stability and Formation of the Li3PS4/Li, Li3PS4/Li2S, and Li2S/Li Interfaces: A Theoretical Study.

Solid electrolytes have shown superior behavior and many advantages over liquid electrolytes, including simplicity in battery design. However, some chemical and structural instability problems arise when solid electrolytes form a direct interface with the negative Li-metal electrode. In particular, it was recognized that the interface between the ß-Li3PS4 crystal and lithium anode is quite unstable and tends to promote structural defects that inhibit the correct functioning of the device. As a possible way out of this problem, we propose a material, Li2S, as a passivating coating for the Li/ß-Li3PS4 interface. We investigated the mutual affinity between Li/Li2S and Li2S/ß-Li3PS4 interfaces by DFT methods and investigated the structural stability through the adhesion energy and mechanical stress. Furthermore, a topological analysis of the electron density identified preferential paths for the migration of Li ions.

From symmetry breaking in the bulk to phase transitions at the surface: a quantum-mechanical exploration of Li6PS5Cl argyrodite superionic conductor.

Lithium superionic conductor electrolytes may enable the safe use of metallic lithium anodes in all-solid-state batteries. The key to a successful application is a high Li conductivity in the electrolyte material, to be achieved through the maintenance of intimate contact with the electrodes and the knowledge of the chemical nature of that contact. In this manuscript, we tackle this issue by a theoretical ab initio approach. Focusing on the Li6PS5Cl, a thiophosphate with high ionic conductivity, we carry on thorough modeling of the surfaces together with the prediction of the thermal and elastic behaviour. Our investigation leads to some new findings: the bulk structure, as reported in the literature, appears to be metastable, with spontaneous symmetry breaking. Moreover, the relevant stoichiometric surfaces identified for stable and metastable crystal structures are not up-down symmetry related and they expose from one side Li2S and LiCl. Surface reconstructions can be interpreted as local phase transitions. We also predict entirely ab initio the morphology of crystallites, charge, and electrostatic potential at surfaces, together with the effect of temperature on structural properties and the elastic behaviour of this material. Such findings may constitute the relevant groundwork for a better understanding of ionic transport in Li-ion conductors at the electrolyte/anode and electrolyte/cathode interfaces.

Computational Characterization of ß-Li3PS4 Solid Electrolyte: From Bulk and Surfaces to Nanocrystals.

The all-solid-state lithium-ion battery is a new class of batteries being developed following today's demand for renewable energy storage, especially for electric cars. The key component of such batteries is the solid-state electrolyte, a technology that promises increased safety and energy density with respect to the traditional liquid electrolytes. In this view, ß-Li3PS4 is emerging as a good solid-state electrolyte candidate due to its stability and ionic conductivity. Despite the number of recent studies on this material, there is still much to understand about its atomic structure, and in particular its surface, a topic that becomes of key relevance for ionic diffusion and chemical stability in grain borders and contact with the other device components. In this study, we performed a density functional study of the structural and electronic properties of ß-Li3PS4 surfaces. Starting from the bulk, we first verified that the thermodynamically stable structure featured slight distortion to the structure. Then, the surfaces were cut along different crystallographic planes and compared with each other. The (100) surface is confirmed as the most stable at T = 298 K, closely followed by (011), (010), and (210). Finally, from the computed surface energies, the Wulff nanocrystals were obtained and it was verified that the growth along the (100) and (011) directions reasonably reproduces the shape of the experimentally observed nanocrystal. With this study, we demonstrate that there are other surfaces besides (100) that are stable and can form interfaces with other components of the battery as well as facilitate the Li-migration according to their porous structures.

Disordered Rock-Salt Type Li2TiS3 as Novel Cathode for LIBs: A Computational Point of View.

The development of high-energy cathode materials for lithium-ion batteries with low content of critical raw materials, such as cobalt and nickel, plays a key role in the progress of lithium-ion batteries technology. In recent works, a novel and promising family of lithium-rich sulfides has received attention. Among the possible structures and arrangement, cubic disordered Li2TiS3 has shown interesting properties, also for the formulation of new cell for all-solid-state batteries. In this work, a computational approach based on DFT hybrid Hamiltonian, localized basis functions and the use of the periodic CRYSTAL code, has been set up. The main goal of the present study is to determine accurate structural, electronic, and spectroscopic properties for this class of materials. Li2TiS3 precursors as Li2S, TiS2, and TiS3 alongside other formulations and structures such as LiTiS2 and monoclinic Li2TiS3 have been selected as benchmark systems and used to build up a consistent and robust predictive scheme. Raman spectra, XRD patterns, electronic band structures, and density of states have been simulated and compared to available literature data. Disordered rock-salt type Li2TiS3 structures have been derived via a solid solution method as implemented into the CRYSTAL code. Representative structures were extensively characterized through the calculations of their electronic and vibrational properties. Furthermore, the correlation between structure and Raman fingerprint was established.

Effect of Internal Donors on Raman and IR Spectroscopic Fingerprints of MgCl2/TiCl4 Nanoclusters Determined by Machine Learning and DFT.

To go deep into the origin of MgCl2 supported Ziegler-Natta catalysis we need to fully understand the structure and properties of precatalytic nanoclusters MgCl2/TiCl4 in presence of Lewis bases as internal donors (ID). In this work MgCl2/TiCl4 nanoplatelets derived by machine learning and DFT calculations have been used to model the interaction with ethyl-benzoate EB as ID, with available exposed sites of binary TixCly/MgCl2 systems. The influence of vicinal Ti2Cl8 and coadsorbed TiCl4 on energetic, structural and spectroscopic behaviour of EB has been considered. The adsorption of homogeneous-like TiCl4EB and TiCl4(EB)2 at the various surface sites have been also simulated. B3LYP-D2 and M06 functionals combined with TZVP quality basis set have been adopted for calculations. The adducts have been characterized by computing IR and Raman spectra that have been found to provide specific fingerprints useful to identify surface species; IR spectra have been successfully compared to available experimental data.

The NV-N+ charged pair in diamond: a quantum-mechanical investigation.

The NV-N+ charged pair in diamond has been investigated by using a Gaussian-type basis set, the B3LYP functional, the supercell scheme and the CRYSTAL code. It turns out that: (i) when the distance between the two defects is larger than 6-7 Å, the properties of the double defect are the superposition of the properties of the individual defects. (ii) The energy required for the reaction NV0 + Ns→ NV- + N+ is roughly -1.3 eV at about 12 Å, irrespective of the basis set and functional adopted, and remains negative at any larger distance. (iii) These results support the observation of a charge transfer mechanism through a Ns→ NV0 donation occurring in the ground state, through a tunnelling process, without irradiation. (iv) The IR spectrum of the two subunits is characterized by specific peaks, that might be used as fingerprints. (v) Calculation of electrostatic interaction permitted an estimate of the effective charge of the defects.

Nitrogen substitutional defects in silicon. A quantum mechanical investigation of the structural, electronic and vibrational properties.

The vibrational infrared (IR) and Raman spectra of seven substitutional defects in bulk silicon are computed, by using the quantum mechanical CRYSTAL code, the supercell scheme, an all electron Gaussian type basis set and the B3LYP functional. The relative stability of various spin states has been evaluated, the geometry optimized, the electronic structure analyzed. The IR and Raman intensities have been evaluated analitically. In all cases the IR spectrum is dominated by a single N peak (or by two or three peaks with very close wavenumbers), whose intensity is at least 20 times larger than the one of any other peak. These peaks fall in the 645-712 cm-1 interval, and a shift of few cm-1 is observed from case to case. The Raman spectrum of all defects is dominated by an extremely intense peak at about 530 cm-1, resulting from the (weak) perturbation of the peak of pristine silicon.

Binding of Nucleic Acid Components to the Serpentinite-Hosted Hydrothermal Mineral Brucite.

The adsorption of nucleic acid components onto the serpentinite-hosted hydrothermal mineral brucite has been investigated experimentally by determining the equilibrium adsorption isotherms in aqueous solution. Thermodynamic characterization of the adsorption data has been performed using the extended triple-layer model (ETLM) to establish a model for the stoichiometry and equilibrium constants of surface complexes. Infrared characterization of the molecule-mineral complexes has helped gain insight into the molecular functional groups directly interacting with the mineral surface. Quantum mechanical calculations have been carried out to identify the possible complexes formed on surfaces by nucleic acid components and their binding configurations on mineral surfaces, both in the presence of water molecules and in water-free conditions. The results indicate that brucite favors adsorption of nucleotides with respect to nucleosides and nucleobases from dilute aqueous environments. The surface of this mineral is able to induce well-defined orientations of the molecules through specific molecule-mineral interactions. This result suggests plausible roles of the mineral brucite in assisting prebiotic molecular self-organization. Furthermore, the detection of the infrared spectroscopic features of such building blocks of life adsorbed on brucite at very low degrees of coverage provides important support to life detection investigations.

On the structure of superbasic (MgO)n sites solvated in a faujasite zeolite.

We report the synthesis and characterisation of a HY/MgO zeolite/oxide nanocomposite material with high crystallinity and highly dispersed, highly basic MgO sites. Preparation was optimized in order to preserve sample crystallinity, to avoid the formation of mesoporosity and to minimize the formation of separate Mg-containing phases. These features were checked by means of electron microscopy, X-ray powder diffraction, porosimetry and IR spectroscopy. A highly dispersed material was obtained, comprising nanoclusters of magnesium oxide and hydroxide hosted by the microporous zeolite framework. The location and structure of the Mg-containing clusters have been studied by means of a combination of Rietveld refinement of XRPD data and high quality quantum mechanical simulations. The refinement has shown the presence of magnesium and oxygen atoms in the double six-membered ring cages, consistent with the presence of mononuclear Mg moieties. However, the composition and IR spectroscopy demonstrate that other Mg species must exist, likely located in the zeolite pores. In order to propose candidate structures for these species, several hypothetic periodic models of the material were built by placing (MgO)n clusters in different locations of the zeolite structure, taking into account the material composition and other constraints imposed by the experimental observations. Periodic structures with P1 symmetry were optimized at the B3LYP-D*/DZVP level with the CRYSTAL code and classified according to their stability. Two families of possible sites were identified: highly solvated (MgO)n units in narrow cavities and less coordinated clusters in the supercages. The stability of these clusters appears to be regulated by the ability of Mg2+ and O2- ions to interact with the pore walls and by the formation of Mg-OH species as result of the reaction of Mg-O couples with remaining acidic protons. The reactivity of four representative models with CO2 has been simulated at the B3LYP-D*/TZVP level. CO2 forms very stable linear end-on adducts with low coordinated Mg ions in most cases. Isolated sites give rise to bridge bidentate complexes in agreement with previous spectroscopic observations. The formation of hydrogen-carbonates is observed only on specific sites, through a process having a low adsorption energy because of the high deformation of the adsorption site.

Probing the coordinative unsaturation and local environment of Ti³âº sites in an activated high-yield Ziegler-Natta catalyst.

The typical activation of a fourth generation Ziegler-Natta catalyst TiCl4/MgCl2/phthalate with triethyl aluminum generates Ti(3+) centers that are investigated by multi-frequency continuous wave and pulse EPR methods. Two families of isolated, molecule-like Ti(3+) species have been identified. A comparison of the experimentally derived g tensors and (35,37)Cl hyperfine and nuclear-quadrupole tensors with DFT-computed values suggests that the dominant EPR-active Ti(3+)  species is located on MgCl2(110) surfaces (or equivalent MgCl2 terminations with tetra-coordinated Mg). O2 reactivity tests show that a fraction of these Ti sites is chemically accessible, an important result in view of the search for the true catalyst active site in olefin polymerization.

Periodic and high-temperature disordered conformations of polytetrafluoroethylene chains: an ab initio modeling.

We report here the main results of a successful attempt to predict some macroscopic properties of representative polymers of technological relevance both in regular and disordered forms by using first principle quantum mechanical approaches at microscopic level. Until now, the prediction of the structural and thermal properties of those polymers has been mostly a domain of molecular mechanics methods. To overcome the limits of those classical computational tools whenever physical properties are significantly influenced by stereoelectronic effects (e.g., electron rich substituents), we employed methods rooted in the Density Functional Theory (DFT). A general computational strategy including the proper choice of periodic boundary conditions (PBC), functional, basis set, and model system size, has been tested and validated for saturated polymers such as polyethylene and isotactic/syndiotactic polypropylenes. On the basis of these results, a comprehensive study of poly(tetrafluoroethylene) (PTFE) chains in both regular periodic and disordered conformations has been performed. A statistical approach has been next applied to obtain the thermal concentration of defects and to reproduce the thermal behavior of the investigated polymer. At the end, a very good agreement with experimental X-ray diffraction and IR results has been achieved, definitely reaching a good understanding of the widely studied disorder phenomenon determining the main technological properties of poly(tetrafluoroethylene) (the trade Teflon) and, at the same time, identifying the proper computational tools to investigate perfluoro-compounds and other complex polymeric systems.