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.

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.

Ab Initio Modeling of MultiWall: A General Algorithm First Applied to Carbon Nanotubes.

A general, versatile and automated computational algorithm to design any type of multiwall nanotubes of any chiralities is presented for the first time. It can be applied to rolling up surfaces obtained from cubic, hexagonal, and orthorhombic lattices. Full exploitation of the helical symmetry permits a drastic reduction of the computational cost and therefore opens to the study of realistic systems. As a test case, the structural, electronic, mechanical, and transport properties of multiwall carbon nanotubes (MWCNT) are calculated using a density functional theory approach, and results are compared with those of the corresponding layered (graphene-like) precursors. The interaction between layers has a general minimum for the inter-wall distance of ≈3.4 Å, in good agreement with experimental and computed optimal distances in graphene sheets. The metallic armchair and semiconductor zigzag MWCNT are almost isoenergetic and their stability increases as the number of walls increases. The vibrational fingerprint provides a reliable tool to identify the chirality and the thickness of the nanostructures. Finally, some promising thermoelectric features of the semiconductor MWCNT are reproduced and discussed.

Topology of the Electron Density and of Its Laplacian from Periodic LCAO Calculations on f-Electron Materials: The Case of Cesium Uranyl Chloride.

The chemistry of f-electrons in lanthanide and actinide materials is yet to be fully rationalized. Quantum-mechanical simulations can provide useful complementary insight to that obtained from experiments. The quantum theory of atoms in molecules and crystals (QTAIMAC), through thorough topological analysis of the electron density (often complemented by that of its Laplacian) constitutes a general and robust theoretical framework to analyze chemical bonding features from a computed wave function. Here, we present the extension of the Topond module (previously limited to work in terms of s-, p- and d-type basis functions only) of the Crystal program to f- and g-type basis functions within the linear combination of atomic orbitals (LCAO) approach. This allows for an effective QTAIMAC analysis of chemical bonding of lanthanide and actinide materials. The new implemented algorithms are applied to the analysis of the spatial distribution of the electron density and its Laplacian of the cesium uranyl chloride, Cs2UO2Cl4, crystal. Discrepancies between the present theoretical description of chemical bonding and that obtained from a previously reconstructed electron density by experimental X-ray diffraction are illustrated and discussed.

A quantum-mechanical investigation of oxygen vacancies and copper doping in the orthorhombic CaSnO3 perovskite.

In this study we explore the implications of oxygen vacancy formation and of copper doping in the orthorhombic CaSnO3 perovskite, by means of density functional theory, focusing on energetic and electronic properties. In particular, the electronic charge distribution is analyzed by Mulliken, Hirshfeld-I, Bader and Wannier approaches. Calculations are performed at the PBE and the PBE0 level (for doping with Cu, only PBE0), with both spin-restricted and spin-unrestricted formulations; unrestricted calculations are used for spin-polarized cases and for the naturally open-shell cases (Cu doping). An oxygen vacancy is found to have the tendency to reduce Sn neighbors by giving rise to an energy band within the energy band-gap of the pristine system, close to the valence band. At variance with what happens in the CaTiO3 perovskite (also investigated here), an oxygen vacancy in the CaSnO3 perovskite is found to lose two valence electrons and thus to be positively charged so that no F-center is formed. Regarding Cu doping, when one Sn atom is substituted by a Cu one, the most stable configuration corresponds to having the Cu atom as a first neighbor to the vacancy. These findings shed some light on the catalytic and phosphorous host properties of this perovskite.

Calibration of 57Fe Mössbauer constants by first principles.

Electron density and eigenvalues of the 3 × 3 matrix of the electric field gradients at the (57)Fe nuclei positions have been evaluated with the periodic ab initio CRYSTAL code for a wide range of crystalline compounds, adopting different computational approaches (Hartree-Fock, gradient corrected and hybrids functionals). The robust calibration procedure, involving experimental isomer shifts and quadrupolar splittings, yields reliable Mössbauer parameters, i.e. the isomer shift constant (α) and the nuclear quadrupolar moment (Q). Dependence of the results on the Hamiltonian is explored and well suited localized basis sets for periodic calculations are provided.

Electron density analysis of large (molecular and periodic) systems: A parallel implementation.

A parallel implementation is presented of a series of algorithms for the evaluation of several one-electron properties of large molecular and periodic (of any dimensionality) systems. The electron charge and momentum densities of the system, the electrostatic potential, X-ray structure factors, directional Compton profiles can be effectively evaluated at low computational cost along with a full topological analysis of the electron charge density (ECD) of the system according to Bader's quantum theory of atoms in molecules. The speedup of the parallelization of the different algorithms is presented. The search of all symmetry-irreducible critical points of the ECD of the crystallized crambin protein and the evaluation of all the corresponding bond paths, for instance, would require about 32 days if run in serial mode and reduces to less than 2 days when run in parallel mode over 32 processors.

Diffraction of helium on MgO(100) surface calculated from first-principles.

In this work we simulate the diffraction peak intensities of He beams scattered on the MgO(100) surface from first principles. It turns out that diffraction peak intensities are extremely sensitive to the quality of the potential describing the He-MgO surface interaction. Achieving the required accuracy in first principles calculations is very challenging indeed. The present work describes a first principles protocol able to achieve very high accuracy for reasonable computational cost. This method is based on periodic local second-order Møller-Plesset perturbation theory where systematic corrections for basis set truncation and for high-order electronic correlation are introduced using coupled cluster calculations on finite model systems mimicking the target system. For the He-MgO system the requirements with respect to the level of theory are very high; it turns out that contributions from connected quadruple excitations are non-negligible. Here we demonstrate that using this protocol, it is possible to reach the accuracy in the He-MgO potential that is required to predict the observed He diffraction peak intensities.

Born-Oppenheimer Molecular Dynamics with a Linear Combination of Atomic Orbitals and Hybrid Functionals for Condensed Matter Simulations Made Possible. Theory and Performance for the Microcanonical and Canonical Ensembles.

The implementation of an original Born-Oppenheimer molecular dynamics module is presented, which is able to perform simulations of large and complex condensed phase systems for sufficiently long time scales at the level of density functional theory with hybrid functionals, in the microcanonical (NVE) and canonical (NVT) ensembles. The algorithm is fully integrated in the Crystal code, a program for quantum mechanical simulations of materials, whose peculiarity stems from the use of atom-centered basis functions within a linear combination of atomic orbitals to describe the wave function. The corresponding efficiency in the evaluation of the exact Fock exchange series has led to the implementation of a rich variety of hybrid density functionals at a low computational cost. In addition, the molecular dynamics implementation benefits also from the effective MPI parallelization of the code, suited to exploit high-performance computing resources available on current generation supercomputer architectures. Furthermore, the information contained in the trajectory of the dynamics is extracted through a series of postprocessing algorithms that provide the radial distribution function, the diffusion coefficient and the vibrational density of states. In this work, we present a detailed description of the theoretical framework and the algorithmic implementation, followed by a critical evaluation of the accuracy and parallel performance (e.g., strong and weak scaling) of this approach, when ice and liquid water simulations are performed in the microcanonical and canonical ensembles.

Crucial Role of Ni Point Defects and Sb Doping for Tailoring the Thermoelectric Properties of ZrNiSn Half-Heusler Alloy: An Ab Initio Study.

In the wide group of thermoelectric compounds, the half-Heusler ZrNiSn alloy is one of the most promising materials thanks to its thermal stability and narrow band gap, which open it to the possibility of mid-temperature applications. A large variety of defects and doping can be introduced in the ZrNiSn crystalline structure, thus allowing researchers to tune the electronic band structure and enhance the thermoelectric performance. Within this picture, theoretical studies of the electronic properties of perfect and defective ZrNiSn structures can help with the comprehension of the relation between the topology of defects and the thermoelectric features. In this work, a half-Heusler ZrNiSn alloy is studied using different defective models by means of an accurate Density Functional Theory supercell approach. In particular, we decided to model the most common defects related to Ni, which are certainly present in the experimental samples, i.e., interstitial and antisite Ni and a substitutional defect consisting of the replacement of Sn with Sb atoms using concentrations of 3% and 6%. First of all, a comprehensive characterization of the one-electron properties is performed in order to gain deeper insight into the relationship between structural, topological and electronic properties. Then, the effects of the modeled defects on the band structure are analyzed, with particular attention paid to the region between the valence and the conduction bands, where the defective models introduce in-gap states with respect to the perfect ZrNiSn crystal. Finally, the electronic transport properties of perfect and defective structures are computed using semi-classical approximation in the framework of the Boltzmann transport theory as implemented in the Crystal code. The dependence obtained of the Seebeck coefficient and the power factor on the temperature and the carrier concentration shows reasonable agreement with respect to the experimental counterpart, allowing possible rationalization of the effect of the modeled defects on the thermoelectric performance of the synthesized samples. As a general conclusion, defect-free ZrNiSn crystal appears to be the best candidate for thermoelectric applications when compared to interstitial and antisite Ni defective models, and substitutional defects of Sn with Sb atoms (using concentrations of 3% and 6%) do not appreciably improve electronic transport properties.

Combined DFT-D3 Computational and Experimental Studies on g-C3N4: New Insight into Structure, Optical, and Vibrational Properties.

Graphitic carbon nitride (g-C3N4) has emerged as one of the most promising solar-light-activated polymeric metal-free semiconductor photocatalysts due to its thermal physicochemical stability but also its characteristics of environmentally friendly and sustainable material. Despite the challenging properties of g-C3N4, its photocatalytic performance is still limited by the low surface area, together with the fast charge recombination phenomena. Hence, many efforts have been focused on overcoming these drawbacks by controlling and improving the synthesis methods. With regard to this, many structures including strands of linearly condensed melamine monomers, which are interconnected by hydrogen bonds, or highly condensed systems, have been proposed. Nevertheless, complete and consistent knowledge of the pristine material has not yet been achieved. Thus, to shed light on the nature of polymerised carbon nitride structures, which are obtained from the well-known direct heating of melamine under mild conditions, we combined the results obtained from XRD analysis, SEM and AFM microscopies, and UV-visible and FTIR spectroscopies with the data from the Density Functional Theory method (DFT). An indirect band gap and the vibrational peaks have been calculated without uncertainty, thus highlighting a mixture of highly condensed g-C3N4 domains embedded in a less condensed "melon-like" framework.

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.

CRYSCOR: a program for the post-Hartree-Fock treatment of periodic systems.

Cryscor is a periodic post-Hartree-Fock program based on local functions in direct space, i.e., Wannier functions and projected atomic orbitals. It uses atom centered Gaussians as basis functions. The Hartree-Fock reference, as well as symmetry information, is provided by the Crystal program. Cryscor presently features an efficient and parallel implementation of periodic local second order Møller-Plesset perturbation theory (MP2), which allows us to study 1D-, 2D- and 3D-periodic systems beyond 1000 basis functions per unit cell. Apart from the correlation energy also the MP2 density matrix, and from that the Compton profile, are available. Very recently, a new module for calculating excitonic band gaps at the uncorrelated Configuration-Interaction-Singles (CIS) level has been added. Other advancements include new extrapolation techniques for calculating surface adsorption on semi-infinite solids. In this paper the diverse features and recent advances of the present Cryscor version are illustrated by exemplary applications to various systems: the adsorption of an argon monolayer on the MgO (100) surface, the rolling energy of a boron nitride nanoscroll, the relative stability of different aluminosilicates, the inclusion energy of methane in methane-ice-clathrates, and the effect of electron correlation on charge and momentum density of α-quartz. Furthermore, we present some first tentative CIS results for excitonic band gaps of simple 3D-crystals, and their dependence on the diffuseness of the basis set.

Modeling of BN-Doped Carbon Nanotube as High-Performance Thermoelectric Materials.

Ternary BNC nanotubes were modeled and characterized through a periodic density functional theory approach with the aim of investigating the influence on the structural, electronic, mechanical, and transport properties of the quantity and pattern of doping. The main energy band gap is easily tunable as a function of the BN percentage, the mechanical stability is generally preserved, and an interesting piezoelectric character emerges in the BNC structures. Moreover, C@(BN)1-xCx double-wall presents promising values of the thermoelectric coefficients due to the combined lowering of the thermal conductivity and increase of charge carriers. Computed results are in qualitative agreement with the little experimental evidence and therefore can provide insights on an atomic scale of the real samples and direct the synthesis towards increasingly performing hybrid nanomaterials.

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.

Charge Density Analysis of Actinide Compounds from the Quantum Theory of Atoms in Molecules and Crystals.

The nature of chemical bonding in actinide compounds (molecular complexes and materials) remains elusive in many respects. A thorough analysis of their electron charge distribution can prove decisive in elucidating bonding trends and oxidation states along the series. However, the accurate determination and robust analysis of the charge density of actinide compounds pose several challenges from both experimental and theoretical perspectives. Significant advances have recently been made on the experimental reconstruction and topological analysis of the charge density of actinide materials [Gianopoulos et al. IUCrJ, 2019, 6, 895]. Here, we discuss complementary advances on the theoretical side, which allow for the accurate determination of the charge density of actinide materials from quantum-mechanical simulations in the bulk. In particular, the extension of the Topond software implementing Bader's quantum theory of atoms in molecules and crystals (QTAIMAC) to f- and g-type basis functions is introduced, which allows for an effective study of lanthanides and actinides in the bulk and in vacuo, on the same grounds. Chemical bonding of the tetraphenyl phosphate uranium hexafluoride cocrystal [PPh4+][UF6-] is investigated, whose experimental charge density is available for comparison. Crystal packing effects on the charge density and chemical bonding are quantified and discussed. The methodology presented here allows reproducing all subtle features of the topology of the Laplacian of the experimental charge density. Such a remarkable qualitative and quantitative agreement represents a strong mutual validation of both approaches-experimental and computational-for charge density analysis of actinide compounds.

DFT and local-MP2 periodic study of the structure and stability of two proton-ordered polymorphs of ice.

The equilibrium geometry and the formation energy of two periodic polymorphs of Ice have been theoretically studied: the former (Ice XI, crystal group Cmc2(1)) is experimentally observed as the most stable structure at low temperature and pressure; the latter (crystal group Pna2(1)) is the simplest proton-ordered model of ordinary ice. With the Crystal code, the problem is solved using Hartree-Fock (HF), pure Kohn-Sham (PW91), or hybrid (B3LYP) one-electron Hamiltonians. The B3LYP results are those in best agreement with the experiment. Using the B3LYP-optimized geometry and starting from the corresponding HF Crystal solution, the energetics of the two polymorphs have been investigated at an ab initio MP2 level using the Cryscor code, based on a local-correlation approach: these calculations have allowed us not only to confirm the excellent B3LYP results as concerns the formation energy and the relative stability of the two structures but also to analyze the role in this respect of the intra- and intermolecular contributions to the correlation energy. Since both Crystal and Cryscor adopt a basis set of localized Gaussian-type functions and since very small energy differences are involved, utter attention has been paid to correcting for the basis set superposition error in the calculation of formation energies.

Periodic density functional theory and local-MP2 study of the librational modes of Ice XI.

Two periodic codes, CRYSTAL and CRYSCOR, are here used to simulate and characterize the librational modes of the nu(R) band of Ice XI: this band has been found experimentally to be the region of the vibrational spectrum of ordinary ice most affected by the transition from the proton-disordered (Ice Ih) to the proton-ordered (Ice XI) phase. With CRYSTAL, the problem is solved using Hartree-Fock (HF), pure Kohn-Sham (PW91) or hybrid (B3LYP) one-electron Hamiltonians: the harmonic approximation is employed to obtain the vibrational spectrum after optimizing the geometry. The B3LYP results are those in best agreement with the experiment. For a given crystalline geometry, CRYSCOR computes the energy per cell in an ab initio HF+MP2 approximation using a local-correlation approach; this technique is employed for recalculating the frequencies of the different modes identified by the B3LYP approach, by fully accounting for long range dispersive interactions. The effect of anharmonicity is evaluated separately for each mode both in the B3LYP and HF+MP2 case. The two approaches accurately reproduce the four-peak structure of the librational band. The harmonic B3LYP nu(R) bandwidth of 70 meV is lowered to 60 meV by anharmonic corrections, and becomes 57 meV in the HF+MP2 anharmonic calculation, in excellent agreement with the experimental IINS data (56-59 meV). The assignment of the librational modes is discussed.