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.

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.

Experimental and computational study of the role of defects and secondary phases on the thermoelectric properties of TiNi1+xSn (0 ≤x≤ 0.12) half Heusler compounds.

The half Heusler TiNiSn compound is a model system for understanding the relationship among structural, electronic, microstructural and thermoelectric properties. However, the role of defects that deviate from the ideal crystal structure is far from being fully described. In this work, TiNi1+xSn alloys (x= 0, 0.03, 0.06, 0.12) were synthesized by arc melting elemental metals and annealed to achieve equilibrium conditions. Experimental values of the Seebeck coefficient and electrical resistivity, obtained from this work and from the literature, scale with the measured carrier concentration, due to different amounts of secondary phases and interstitial nickel. Density functional theory calculations showed that the presence of both interstitial Ni defects and composition conserving defects narrows the band gap with respect to the defect free structure, affecting the transport properties. Accordingly, results of experimental investigations have been explained confirming that interstitial Ni defects, as well as secondary phases, promote a metallic behavior, raising the electrical conductivity and lowering the absolute values of the Seebeck coefficient.

Roadmap on thermoelectricity.

The increasing energy demand and the ever more pressing need for clean technologies of energy conversion pose one of the most urgent and complicated issues of our age. Thermoelectricity, namely the direct conversion of waste heat into electricity, is a promising technique based on a long-standing physical phenomenon, which still has not fully developed its potential, mainly due to the low efficiency of the process. In order to improve the thermoelectric performance, a huge effort is being made by physicists, materials scientists and engineers, with the primary aims of better understanding the fundamental issues ruling the improvement of the thermoelectric figure of merit, and finally building the most efficient thermoelectric devices. In this Roadmap an overview is given about the most recent experimental and computational results obtained within the Italian research community on the optimization of composition and morphology of some thermoelectric materials, as well as on the design of thermoelectric and hybrid thermoelectric/photovoltaic devices.

First-principles evaluation of the second harmonic generation response of reference organic and inorganic crystals.

Using the CRYSTAL17 package at the coupled-perturbed Kohn-Sham (CPKS) level, periodic boundary conditions first-principles calculations are enacted to predict the second harmonic generation second-order nonlinear optical (NLO) susceptibility, χ(2), values of six historical NLO crystals. This selection allowed the comparison between state-of-the-art calculations and experiment. Several computational aspects are tackled to define conditions where the results are converged with respect to the range of lattice summations, to the number of k-points in the first Brillouin zone, to the order of the multipole expansions for evaluating the long-range part of the electrostatic interactions, as well as to the atomic basis set size. A valence triple zeta basis set supplemented with polarization functions has been selected. Then, χ(2) calculations have been performed using a range of exchange-correlation functionals (XCFs). Results show the large impact of the amount of Hartree-Fock (HF) exchange on the amplitude but also on the sign on the χ(2) tensor components. To a given extent, these amplitude effects are consistent with results on molecules, but the sign reversal effects and the non-monotonic behavior of the χ(2) tensor components as a function of the amount of HF exchange are scarcely found for molecules. Then, using the recommended range-separated hybrid XCFs, the CPKS scheme leads to good agreement with experimental data for potassium dihydrogenophosphate, urea, and χZXX (2) of LiNbO3. The agreement is more questionable for χZZZ (2) of LiNbO3 whereas it remains poor for ammonium dihydrogenophosphate and 2-methyl-4-nitroaniline, with large underestimations by about a factor of 3, opening a path to further fine-tuning of the ranges of inclusion of HF exchange.

CRYSTAL23: A Program for Computational Solid State Physics and Chemistry.

The Crystal program for quantum-mechanical simulations of materials has been bridging the realm of molecular quantum chemistry to the realm of solid state physics for many years, since its first public version released back in 1988. This peculiarity stems from the use of atom-centered basis functions within a linear combination of atomic orbitals (LCAO) approach and from the corresponding efficiency in the evaluation of the exact Fock exchange series. In particular, this has led to the implementation of a rich variety of hybrid density functional approximations since 1998. Nowadays, it is acknowledged by a broad community of solid state chemists and physicists that the inclusion of a fraction of Fock exchange in the exchange-correlation potential of the density functional theory is key to a better description of many properties of materials (electronic, magnetic, mechanical, spintronic, lattice-dynamical, etc.). Here, the main developments made to the program in the last five years (i.e., since the previous release, Crystal17) are presented and some of their most noteworthy applications reviewed.

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.

Reaction barriers on non-conducting surfaces beyond periodic local MP2: Diffusion of hydrogen on α-Al2O3(0001) as a test case.

The quest for "chemical accuracy" is becoming more and more demanded in the field of structure and kinetics of molecules at solid surfaces. In this paper, as an example, we focus on the barrier for hydrogen diffusion on a α-Al2O3(0001) surface, aiming for a couple cluster singles, doubles, and perturbative triples [CCSD(T)]-level benchmark. We employ the density functional theory (DFT) optimized minimum and transition state structures reported by Heiden, Usvyat, and Saalfrank [J. Phys. Chem. C 123, 6675 (2019)]. The barrier is first evaluated at the periodic Hartree-Fock and local Møller-Plesset second-order perturbation (MP2) level of theory. The possible sources of errors are then analyzed, which includes basis set incompleteness error, frozen core, density fitting, local approximation errors, as well as the MP2 method error. Using periodic and embedded fragment models, corrections to these errors are evaluated. In particular, two corrections are found to be non-negligible (both from the chemical accuracy perspective and at the scale of the barrier value of 0.72 eV): the correction to the frozen core-approximation of 0.06 eV and the CCSD(T) correction of 0.07 eV. Our correlated wave function results are compared to barriers obtained from DFT. Among the tested DFT functionals, the best performing for this barrier is B3LYP-D3.

Electronic Excitations in Crystalline Solids through the Maximum Overlap Method.

The maximum overlap method (MOM) has emerged from molecular quantum chemistry as a convenient practical procedure for studying excited states. Unlike the Aufbau principle, during self-consistent field (SCF) iterations, the MOM forces orbital occupation to be maximally similar to that of a reference state. Although still within a single-particle framework, this approach allows for the evaluation of excitation energies (Δ-SCF) and geometry optimization of electronic configurations other than the ground state. In this work, we present an extension of the MOM to periodic crystalline solids, within the framework of an atom-centered Gaussian basis set. In order to obtain a realistic concentration of excited electrons, we allow excitation in only one-or a few-points of the Brillouin zone, leading to a fractional occupation of crystalline Kohn-Sham states. Since periodic SCF solution techniques involve an iteration between direct and reciprocal spaces, only totally symmetric excitations are allowed in our treatment, in order to preserve the translational symmetry: vertical Γ-point excitations or collective excitations in a sphere around Γ. Other types of excitations are accessible through folding of the Brillouin zone subsequent to the creation of a supercell. The features and performance of the method are presented through its application to prototypical solids such as bulk silicon, diamond, and lithium fluoride and comparing the results with the available experimental data. The demonstrative application to nickel oxide and solid CuI(piperazine)-a luminescent copper halide compound-highlights the promising potential of the MOM in solid-state quantum chemistry.

Fragment-Based Restricted Active Space Configuration Interaction with Second-Order Corrections Embedded in Periodic Hartree-Fock Wave Function.

We present a computational scheme for restricted-active-space configuration interaction (RASCI) calculations combined with second-order perturbation theory (RASCI-PT2) on a fragment of a periodic system embedded in the periodic Hartree-Fock (HF) wave function. This method allows one to calculate the electronic structure of localized strongly correlated features in crystals and surfaces. The scheme was implemented via an interface between the Cryscor and Q-Chem codes. To evaluate the performance of the embedding method, we explored dissociation of a fluorine atom from a lithium fluoride surface and partially fluorinated graphane layer. The results show that RASCI and RASCI-PT2 embedded in periodic HF are able to produce well-behaved potential energy surfaces and accurate dissociation energies.

Exfoliation Energy of Layered Materials by DFT-D: Beware of Dispersion!

In this work, we have computed the exfoliation energy (the energy required to separate a single layer from the bulk structure), the interlayer distance, and the thermodynamic state functions for representative layered inorganic minerals such as Brucite, Portlandite, and Kaolinite, while leaving the more classical 2D transition-metal dichalcogenides (like MoS2) for future work. Such materials are interesting for several applications in the field of adsorption and in prebiotic chemistry. Their peculiar features are directly controlled by the exfoliation energy. In materials without cations/anions linking different layers, the interactions keeping the layers together are of weak nature, mainly dispersion London interactions and hydrogen bonds, somehow challenging to deal with computationally. We used Hartree-Fock (HF) and density functional theory (DFT) approaches focusing on the role of dispersion forces using the popular and widespread Grimme's pairwise dispersion schemes (-D2 and -D3) and, as a reference method, the periodic MP2 approach based on localized orbitals (LMP2). The results are highly dependent on the choice of the scheme adopted to account for dispersion interactions. D2 and D3 schemes combined with either HF or DFT lead to overestimated exfoliation energies (about 2.5 and 1.7 times higher than LMP2 data for Brucite/Portlandite and Kaolinite) and underestimated interlayer distances (by about 3.5% for Brucite/Portlandite). The reason is that D2 and D3 corrections are based on neutral atomic parameters for each chemical element which, instead, behave as cations in the considered layered material (Mg, Ca, and Al), causing overattractive interaction between layers. More sophisticated dispersion corrections methods, like those based on nonlocal vdW functionals, many body dispersion model, and exchange-hole dipole moment not relying on atom-typing, are, in principle, better suited to describe the London interaction of ionic species. Nonetheless, we demonstrate that good results can be achieved also within the simpler D2 and D3 schemes, in agreement with previous literature suggestions, by adopting the dispersion coefficients of the preceding noble gas for the ionic species, leading to energetics in good agreement with LMP2 and structures closer to the experiments.

Gaussian Basis Sets for Crystalline Solids: All-Purpose Basis Set Libraries vs System-Specific Optimizations.

It is customary in molecular quantum chemistry to adopt basis set libraries in which the basis sets are classified according to either their size (triple-ζ, quadruple-ζ, ...) and the method/property they are optimal for (correlation-consistent, linear-response, ...) but not according to the chemistry of the system to be studied. In fact the vast majority of molecules is quite homogeneous in terms of density (i.e., atomic distances) and types of bond involved (covalent or dispersive). The situation is not the same for solids, in which the same chemical element can be found having metallic, ionic, covalent, or dispersively bound character in different crystalline forms or compounds, with different packings. This situation calls for a different approach to the choice of basis sets, namely a system-specific optimization of the basis set that requires a practical algorithm that could be used on a routine basis. In this work we develop a basis set optimization method based on an algorithm-similar to the direct inversion in the iterative subspace-that we name BDIIS. The total energy of the system is minimized together with the condition number of the overlap matrix as proposed by VandeVondele et al. [VandeVondele et al. J. Chem. Phys. 2007, 227, 114105]. The details of the method are here presented, and its performance in optimizing valence orbitals is shown. As demonstrative systems we consider simple prototypical solids such as diamond, graphene sodium chloride, and LiH, and we show how basis set optimizations have certain advantages also toward the use of large (quadruple-ζ) basis sets in solids, both at the DFT and Hartree-Fock level.

Coupled Perturbation Theory Approach to Dual Basis Sets for Molecules and Solids. 1. General Theory and Application to Molecules.

We present a new coupled Hartree-Fock(HF)/Kohn-Sham DFT perturbation method that accounts for the effect of enlarging the basis set in electronic structure calculations. In contrast with previous approaches, our dual basis set treatment yields not only a correction for the total energy but also for the orbital eigenvalues and density. The zeroth order solution is obtained from the projection of the small basis set coefficients. Diagonalization of the full Fock matrix in the large basis set is avoided. In this first paper of a series, we develop the theoretical foundations of our approach for molecules, including the coupled-perturbed equations through second order and the energy expressions through fourth order-as our method complies with Wigner's 2n + 1 rule. The first-order perturbation equation turns out to be uncoupled, and odd-order terms in the energy expansion vanish. In calculations on simple molecules, our method recovers over 93% (84%) of the missing DFT(HF) energy when going from the cc-pVDZ to the aug-cc-pVDZ basis, and over about 95% in all cases if an energy extrapolation formula is used. Mulliken charges, the orbital eigenvalue spectrum, and HOMO-LUMO gaps of the large basis are well reproduced. Charge density maps show that the differences between the perturbatively corrected density and the reference nearly vanish through second-order.

Evolutionary Algorithm-Based Crystal Structure Prediction for Copper(I) Fluoride.

Despite numerous experimental studies since 1824, the binary copper(I) fluoride remains unknown. A crystal structure prediction has been carried out for CuF using the USPEX evolutionary algorithm and a dispersion-corrected hybrid density functional method. In total about 5000 hypothetical structures were investigated. The energetics of the predicted structures were also counter-checked with local second-order Møller-Plesset perturbation theory. Herein 39 new hypothetical copper(I) fluoride structures are reported that are lower in energy compared to the previously predicted cinnabar-type structure. Cuprophilic Cu-Cu interactions are present in all the low-energy structures, leading to ordered Cu substructures such as helical or zig-zag-type Cu-Cu motifs. The lowest-energy structure adopts a trigonal crystal structure with space group P31 21. From an electronic point of view, the predicted CuF modification is a semiconductor with an indirect band gap of 2.3 eV.

Thermoelectric Properties of p-Type Cu2O, CuO, and NiO from Hybrid Density Functional Theory.

The electronic transport coefficients of three Earth-abundant metal oxides Cu2O, CuO, and NiO were investigated using hybrid density functional theory (DFT). Hybrid DFT methods combined with local Gaussian-type basis sets enabled band structure studies on both non-magnetic and magnetic p-type metal oxides without empirical corrections. The CRYSTAL code was used for obtaining the wavefunction, and the transport properties were calculated with two different methodologies to benchmark their accuracy: a numerical approach as implemented in the BoltzTraP code and an analytical approach recently implemented in CRYSTAL17. Both computational methods produce identical results in good agreement with experimental measurements of the Seebeck coefficient. The predicted electrical conductivities are overestimated, owing likely to the used approximation of a constant electronic relaxation time in the calculations, as explicit electron scattering is neglected and relaxation time is considered only as a free parameter. The obtained results enable us to critically review and complement the available theoretical and experimental literature on the studied p-type thermoelectric metal oxide materials.

On the exfoliation and anisotropic thermal expansion of black phosphorus.

Black phosphorus is a bulk solid allotrope of elemental phosphorus and can be seen as an infinite stack of phosphorene sheets. It is interesting from a technological point of view as well as from an electronic structure perspective due to the importance of electron correlation effects. In a recent paper [M. Schütz, L. Maschio, A. J. Karttunen and D. Usvyat, J. Phys. Chem. Lett., 2017, 8, 1290-1294] a highly accurate exfoliation energy has been computed. Building upon these results we carefully benchmark various dispersion-corrected density functional approximations. The choice of the range-separating function that suppresses London dispersion at short interatomic distances apparently has a substantial influence on the results. Having chosen the suitable functional, we have computed the thermal expansion coefficients of black phosphorous via a quasi-harmonic approximation. The computed coefficients manifest a strong anisotropy between the two in-plane directions. Our calculations, however, do not support the existence of negative thermal expansion in black phosphorus, as reported in some theoretical studies.

Ab initio electronic transport and thermoelectric properties of solids from full and range-separated hybrid functionals.

Within the semiclassical Boltzmann transport theory in the constant relaxation-time approximation, we perform an ab initio study of the transport properties of selected systems, including crystalline solids and nanostructures. A local (Gaussian) basis set is adopted and exploited to analytically evaluate band velocities as well as to access full and range-separated hybrid functionals (such as B3LYP, PBE0, or HSE06) at a moderate computational cost. As a consequence of the analytical derivative, our approach is computationally efficient and does not suffer from problems related to band crossings. We investigate and compare the performance of a variety of hybrid functionals in evaluating Boltzmann conductivity. Demonstrative examples include silicon and aluminum bulk crystals as well as two thermoelectric materials (CoSb3, Bi2Te3). We observe that hybrid functionals other than providing more realistic bandgaps-as expected-lead to larger bandwidths and hence allow for a better estimate of transport properties, also in metallic systems. As a nanostructure prototype, we also investigate conductivity in boron-nitride (BN) substituted graphene, in which nanoribbons (nanoroads) alternate with BN ones.

One-Dimensional Phosphorus Nanostructures: from Nanorings to Nanohelices.

Phosphorus nanorings and nanohelices-speculated to exist over 20 years ago-have been systematically derived from one parent structure and studied with quantum chemical methods. The (P8 P2 )n nanorings have been recently synthetized inside carbon nanotube templates, and our comprehensive analysis of possible structural arrangements strongly supports the possibility to experimentally determine the closely related (P8 P2 )n nanohelices. The nanohelices possess very low stiffness, suggesting interesting mechanical properties with nano-spring-like behavior.

One Step Toward a New Generation of C-MOS Compatible Oxide P-N Junctions: Structure of the LSMO/ZnO Interface Elucidated by an Experimental and Theoretical Synergic Work.

Heterostructures formed by La0.7Sr0.3MnO3/ZnO (LSMO/ZnO) interfaces exhibit extremely interesting electronic properties making them promising candidates for novel oxide p-n junctions, with multifunctional features. In this work, the structure of the interface is studied through a combined experimental/theoretical approach. Heterostructures were grown epitaxially and homogeneously on 4″ silicon wafers, characterized by advanced electron microscopy imaging and spectroscopy and simulated by ab initio density functional theory calculations. The simulation results suggest that the most stable interface configuration is composed of the (001) face of LSMO, with the LaO planes exposed, in contact with the (112̅0) face of ZnO. The ab initio predictions agree well with experimental high-angle annular dark field scanning transmission electron microscopy images and confirm the validity of the suggested structural model. Electron energy loss spectroscopy confirms the atomic sharpness of the interface. From statistical parameter estimation theory, it has been found that the distances between the interfacial planes are displaced from the respective ones of the bulk material. This can be ascribed to the strain induced by the mismatch between the lattices of the two materials employed.