Norm-conserving Pseudopotentials Research Articles

Coherent spin transport in a copper protein.

The coherent electron/spin transport in azurin, a species of copper protein, was calculated based on the Landauer model. The research is motivated by the fast electron transport and spin selectivity/polarization in azurin, which have been reported in relation to the chiral-induced spin selectivity of the peptide structure. The calculated spin polarization of copper proteins was large. This phenomenon was strongly influenced by the spin density of the atoms in the ligand group, whereas the contribution of copper was negligible. The results suggest that spin polarization in copper proteins is enhanced by that of the ligand groups. The predicted spin polarization aligns primarily with the scanning tunneling microscope-based break-junction technique to study the electronic properties of single-molecule junctions. Computational techniques employed in this study are nonequilibrium Green's functions (NEGF) and density functional theory (DFT) based on the Landauer model, implemented using the QuantumATK software (Synopsys Inc.). The Perdew-Burke-Ernzerhof (PBE) exchange-correlation functional was adopted for spin-polarized generalized gradient approximation (SGGA). The valence atomic orbitals were constructed using the wavefunctions of the SIESTA package, which was based on the norm-conserving Troullier-Martins relativistic pseudopotentials for describing core electrons. The mesh used for real-space integration was 150 Ha.

Elastic and Electronically Inelastic Cross Sections for the Scattering of Electrons by Pyrrole.

Integral and differential cross sections for elastic and electronically inelastic electron scattering from the pyrrole molecule are reported. The cross section calculations employed the Schwinger multichannel method with norm-conserving pseudopotentials. The collision dynamics was described according to a model in which up to 209 energetically accessible channels were treated as open. In the elastic channel, calculations carried out in the interval of energies from 0 to 50 eV revealed the presence of four resonances with peaks located at 2.56 eV (π1*), 3.82 eV (π2*), 4.70 eV (σNH*), and between 8.30 and 9.50 eV (σ*) positions which are in good agreement with previous assignments. Moreover, the role of the multichannel coupling effects in obtaining accurate cross sections was evaluated by comparing the present results with theoretical results recently reported in the literature and early measurements performed for elastic electron collisions with furan. Electronic excitation cross sections involving the transitions from ground state to the 13B2, 13A1, 11A2, and 11A1 excited states of pyrrole driven by electron impact are presented for energies from thresholds up to 50 eV and, whenever possible, critically compared with the data available in the literature.

Optoelectronic and mechanical properties of gallium arsenide alloys: Based on density functional theory

First principles calculations based on density functional theory (DFT) were performed to investigate the structural, electronic, optical and mechanical properties of pristine GaAs compound and its alloy; Ga0.75Al0.25As, Ga0.75In0.25As, Ga0.75Sn0.25As, Ga0.75Ti0.25As. WIEN2K and Quantum expresso (QE) codes were adopted for calculations using generalized gradient approximation (GGA) in Perdew-Burke Erzenhoff (PBE) as exchange correlation function for both codes. Full potential linear augmented plane wave (FPLAPW) with the local orbital method was adopted as implement in WIEN2K code. In QE code, norm-conserving pseudopotentials were employed on a plane-wave expansion of the wave functions. Structural and electronic properties were elaborated since their result gives information about the optical and mechanical performance. Electronic band structure and optical parameters were performed using WIEN2K code. Underestimation of band gap observed from DFT calculations were corrected by using Modified Becke and Johnson (mBJ). Mechanical components were determined using QE with thermo_pw package. Lattice constant, volume, bulk modulus and other physical parameters were calculated for structural properties. Discrepancy in these parameters as observed in crystal structure is associated to difference in ionic radius of host and substituted atom. The results of band structure and density of states were calculated for electronic properties. All the studied compounds were semiconductors in nature except Ga0.75Sn0.25As which displayed metallic character. Optical parameters including extinction coefficient, absorption coefficient, refractive index, optical conductivity, optical reflectivity and energy loss function have been computed from the dielectric function at energy range of 0 to 25 eV using the Kramers-Kronig transformations. Calculated elastic function were used to compute the mechanical properties such as anisotropic, brittle characteristics, stiffness and many others. All the results were compared with available theoretical and experimental records.

Dichotomy of the description of the elastic and electronically inelastic scattering of electrons by methane

The interaction between electrons and methane is studied with a particular focus on describing the dynamics of elastic and electronic excitation processes under the influence of the multichannel coupling effects. Elastic and electronically inelastic integral and differential cross sections are reported. These cross sections were calculated with the Schwinger multichannel method implemented with norm-conserving pseudopotentials within the minimal orbital basis for single-configuration interaction approach with up to 181 open channels. While an excellent agreement with previously published data is found in the elastic channel, the comparisons involving electronic excitation cross sections reveal significant discrepancies. Furthermore, the lack of reliable experimental data regarding the electronically inelastic channels hampers a more detailed analysis of the theoretical results. Additionally, total ionization and total cross sections for electron collisions with methane are also presented. By highlighting the current theoretical challenges in dealing with electron-induced electronic excitation of molecules, we seek to contribute to a deeper understanding of this scattering process and, thus, pave the way for future research.

Hydrogen storage: Investigation of the elastic properties of Mg7NbH16 hydride

Structural, mechanical and electronic structure properties of Mg7NbH16 were investigated using the norm-conserving pseudo-potentials and plane waves (PP-PW) within the general gradient approximation (GGA) in the frame of density functional theory (DFT). The computed equilibrium lattice constants are in excellent agreement with the available experimental and theoretical data. The elastic constants Cij of Mg7NbH16 were calculated for the first time. The Mg7NbH16 compound is found to be mechanically stable. The bulk modulus of single crystals has been derived using elastic constants. The results are discussed and compared with the calculated bulk modulus reported in the literature. The polycrystalline elastic moduli (namely: the shear modulus, Young's modulus, Poisson's ratio, Lamé's coefficients, sound velocities and the Debye temperature) were derived from the obtained single-crystal elastic constants. According to the obtained results, Mg7NbH16 can be classified as ductile material. The shear anisotropic factors and the elastic anisotropy are also discussed. A Debye temperature of 619 K was also determined using theoretical elastic constants.

Ab initio study of elastic, electronic, and vibrational properties of SnTe and PbTe.

The elastic, electronic, and vibrational properties of the ground state of the rocksalt SnTe and PbTe are investigated. The deduced elastic constants, namely, shear modulus, Young's modulus, and Poisson's ratio are in very good agreement with the experimental and other theoretical data. The electronic band structure and density of states are obtained with and without considering the spin-orbit coupling. The bandgaps of SnTe and PbTe with (without) spin-orbit coupling are 0.11 (0.05) eV and 0.01 (0.78) eV, respectively. The bandgaps with spin-orbit interactions are nearer to experimental data. The hybrid functionals give higher values of bandgaps for both the SnTe and PbTe. In both compounds, the bandgap increases with volume. The valence bandwidths, however, decrease with increasing volume. The vibrational frequencies are found in reasonable agreement with the experiment. The frequencies increase with pressure. In this work, the ab initio calculations of SnTe and PbTe crystals are carried out applying plane wave pseudopotential method using the QUANTUM ESPRESSO package. The PBE exchange and correlation functional based on GGA is considered. The fully relativistic norm-conserving pseudopotentials for Sn, Pb, and Te are used. The self-consistent field calculations are performed over a dense MP net of 18 × 18 × 18k-points. The energy cut-off of 70 Ryd was found sufficient to achieve convergence of 10-6 Ryd in total energy of the crystals.

First-Principles Calculation of Structural, Electronic, and Optical Properties of Cubic Perovskite CsPbF3

Lead halide perovskites have attracted considerable attention as one of the most promising materials for optoelectronic applications. The structural, electronic, and optical properties of the cubic perovskite CsPbF3 were studied using density functional theory in conjunction with plane waves, norm-conserving pseudopotentials, and Perdew-Berg-Erzenhof flavor of generalized gradient approximation. The obtained structural parameters are a good agreement with the experimentally measured and other’s theoretically predicted values. The obtained electronic band structure revealed that cubic CsPbF3 has a direct fundamental band gap of 2.99 eV at point R. The calculated energy band gaps at the high symmetry points agree with the other available theoretical results. The GW method is adapted to correct the underestimated fundamental energy gap value to 4.05 eV. The contribution of the different bands was analyzed from the total and partial density of states. The electron densities show that Cs and F have strong ionic bonds, whereas Pb and F have strong covalent bonds. The optical properties of CsPbF3 were calculated using the density functional perturbation theory and Kramers-Kronig relations. The wide and direct bandgap nature and the calculated optical properties imply that cubic CsPbF3 can be used in optical and optoelectronic devices for high frequencies visible and low frequencies ultraviolet electromagnetic radiation.

On the excitonic effects of the 1T and 1O[formula omitted] phases of PdS2, PdSe2, and PdSSe monolayers

We report the excitonic effects in the 1T and 1OT phases of PdS2, PdSe2, and PdSSe monolayers. Using fully relativistic norm-conserving non-local PBE pseudopotentials and including the spin–orbit (SOC) coupling, the 1T geometries are indirect band gap semiconductors; conversely, the 1OT structures are semimetallic. The band gaps were corrected using the HSE06 hybrid functional, making all structures semiconductors. Maximally localized Wannier functions were utilized to fit the bands of all six structures. The excitonic band structures were obtained using the Bethe–Salpeter equation (BSE), and the absorption coefficients were estimated at both the BSE and the independent particle approximation (IPA) levels. While the 1OT monolayers are more inclined to absorb light polarized in the y-direction, the 1T structures absorb light equally, independent of the polarization direction. We perceive that the 1T phase is the most promising candidate for photovoltaic devices; in particular, the 1T PdSSe system possesses a maximum power conversion efficiency (PCE) of 29.31% if light trapping techniques are used to enhance absorbance. Amongst the 1OT materials, PdS2 features the highest PCE of 8.38%. Our results constitute the first step to encourage experimental studies towards the fabrication of Pd-based regular and Janus transition metal dichalcogenides photoabsorbers and give hints of how to increase the PCE of such materials.

Negative Thermal Expansion in the Polymorphic Modification of Double Sulfate β-AEu(SO4)2 (A-Rb+, Cs+).

New polymorphic modifications of double sulfates β-AEu(SO4)2 (A-Rb+, Cs+) were obtained by the hydrothermal method, the structure of which differs significantly from the monoclinic modifications obtained earlier by solid-state methods. According to single-crystal diffraction data, it was found that the compounds crystallize in the orthorhombic system, space group Pnna, with parameters β-RbEu(SO4)2: a = 9.4667(4) Å, b = 13.0786(5) Å, c = 5.3760(2) Å, V = 665.61(5) Å3; β-CsEu(SO4)2: a = 9.5278(5) Å, b = 13.8385(7) Å, c = 5.3783(3) Å, V = 709.13(7) Å3. The asymmetric part of the unit cell contains one-half Rb+/Cs+ ion, one-half Eu3+ ion, both in special sites, and one SO42- ion. Both compounds exhibit nonlinear negative thermal expansion. According to the X-ray structural analysis and theoretical calculations, the polarizing effect of the alkali metal ion has a decisive influence on the demonstration of this phenomenon. Experimental indirect band gaps of β-Rb and β-Cs are 4.05 and 4.11 eV, respectively, while the direct band gaps are 4.48 and 4.54 eV, respectively. The best agreement with theoretical calculations is obtained using the ABINIT package employing PAW pseudopotentials with hybrid PBE0 functional, while norm-conserving pseudopotentials used in the frame of CASTEP code and LCAO approach in the Crystal package gave worse agreement. The properties of alkali ions also significantly affect the luminescent properties of the compounds, which leads to a strong temperature dependence of the intensity of the 5D0 → 7F4 transition in β-CsEu(SO4)2 in contrast to much weaker dependence of this kind in β-RbEu(SO4)2.

First-principles calculations to investigate structural, optical and electronic properties of ZrO2, Zr0.93Si0.07O2 and Zr0.86Si0.14O2 for dye-sensitised solar cells applications

ABSTRACT Initially, density approximations and generalised gradient approximations (GGA) functional of DFT were executed to determine the electronic structure where the non-local functionals of GGA obtained more close value of reference value, so that GGA method was used to calculate the band structures, DOS, PDOS and optical properties of ZrO2, Zr0.93Si0.07O2 and Zr0.86Si0.14O2 using the Perdew–Burke–Ernzerhof (PBE), Revised Perdew–Burke–Ernzerhof (RPBE), Perdew Wang (PW91), and Wu-Cohen. The transferability of electrons is completely related to PseudoPotential (PP), where four methods, such as OTFG ultra-soft, OTFG norm-conserving, ultra-soft, and norm-conserving, were used in this study. Finally, the norm-conserving PP has been selected as the most accurate method for determination of electronic structure. In addition, to assess the nature of orbital, the DOS and PDOS were additionally simulated and even the optical properties of crystals were calculated. Furthermore, the chemical reactivity, such as HOMO and LUMO, quantum properties, is justified where doping can have a vast effect as well as the open circuit voltage (Voc) is increased after doping. Finally, the toxicity data from in silico study refers to us that these are non-toxic or non-carcinogenic materials with not readily biodegradable status.

Density functionl theory study of Cs2InAgCl(6-x)Brx (x = 0–3) halide double-perovskites

Halide double perovskites have recently received much attention due to their application in optoelectronic devices. Cs2InAgCl6 is an exciting compound with a cubic unit cell, Fm-3m space group, and a direct band gap of about 3 eV. This research studied the effect of Br ion doping instead of Cl with the combination of Cs2InAgCl(6-x)Brx (x values between 0 and 3) using density functional theory (DFT) calculations with GGA approximation and norm-conserving pseudopotentials of PBE type. We found that with an increase in Br concentration, the investigated compounds retain their cubic phase. The compounds' calculated lattice constants and band gaps were about 10.44–10.71 Å and 2.12–0.62 eV, respectively. The accurate band gap values were obtained between 3.16 and 1.68 by HSE06 hybrid functional as exchange-correlation. The calculations of the band structure and density of states indicate the direct electron transition from the valence band to the conduction band at the gamma point of the Brillouin zone of the studied compounds. The studied compounds' optical properties show high absorption in the visible and ultraviolet regions, appropriate optical conductivity, and refractive index.

A DFT study of structural, electronic and optical properties of Lead-free and Ge based cubic perovskite RbGeX‌3 (X= I, Br and Cl)

The current research uses density functional theory (DFT) approximations in conjunction with the plane wave-pseudopotential method to investigate structural, electronic, and optical properties of Pb-free cubic perovskite RbGeX3 (X= I, Br and Cl) materials. More specifically, Norm-conserving pseudopotential has been employed to describe the ion and valence electrons interaction, and Perdew-Burke-Ernzerhof (PBE) flavor is used to represent the exchange-correlation part of the energy of the GGA approximation. Our calculated lattice constants are 5.95, 5.55, and 5.29 Å for RbGeX3 (where X=I, Br, and Cl), respectively, and they are of are in good agreement with available empirical and other values. The band structure shows the direct band gap nature of the three compounds under research here and our values of the band gap energy E_g are in good agreement with the other available results. Materials under research show response to the electromagnetic radiation starting from the infrared region to the very high energies (~33 eV). The RbGeI3 has the lowest E_g value at the low region energies and the highest optical response peaks but RbGeCl3 has the highest optical response peaks at energies located near ~20 eV. Our results show that these materials are good candidates for photo electronic applications including solar cells.

Computation of forces and stresses in solids: Towards accurate structural optimization with auxiliary-field quantum Monte Carlo

The accurate computation of forces and other energy derivatives has been a long-standing challenge for quantum Monte Carlo methods. A number of technical obstacles contribute to this challenge. We discuss how these obstacles can be removed with the auxiliary-field quantum Monte Carlo (AFQMC) approach. AFQMC is a general, high-accuracy, many-body total-energy method for molecules and solids. The implementation of back-propagation for pure estimators allows direct calculation of gradients of the energy via the Hellmann-Feynman theorem. A planewave basis with norm-conserving pseudopotentials is used for the study of periodic bulk materials. Completeness of the planewave basis minimizes the effect of so-called Pulay terms. The ionic pseudopotentials, which can be incorporated in AFQMC in exactly the same manner as in standard independent-electron methods, regulate the force and stress estimators and eliminate any potential divergence of the Monte Carlo variances. The resulting approach allows applications of full geometry optimizations in bulk materials. It also paves the way for many-body computations of the phonon spectrum in solids.

Pseudopotential Bethe-Salpeter calculations for shallow-core x-ray absorption near-edge structures: Excitonic effects in α−Al2O3

We present an ab initio description of optical and shallow-core x-ray absorption spectroscopies in a unified formalism based on the pseudopotential plane-wave method at the level of the Bethe-Salpeter equation (BSE) within Green's functions theory. We show that norm-conserving pseudopotentials are reliable and accurate not only for valence, but also for semicore electron excitations. In order to validate our approach, we compare BSE absorption spectra obtained with two different codes: the pseudopotential-based code EXC and the all-electron full-potential code Exciting. We take corundum $\alpha$-Al$_2$O$_3$ as an example, being a prototypical material that presents strong electron-hole interactions for both valence and core electron excitations. We analyze in detail the optical absorption spectrum as well as the Al L$_1$ and L$_{2,3}$ edges in terms of anisotropy, crystal local fields, interference and excitonic effects. We perform a thorough inspection of the origin and localization of the lowest-energy excitons, and conclude highlighting the purely electronic character off the pre-edge of L$_1$ and the dichroic nature of the optical and L$_{23}$ spectra.

DFT investigation of lattice parameters, spontaneous polarization, and refractive indices of tetragonal BaTiO3 and PbTiO3

In the present work, we study the structural, electrical and optical properties of barium titanate (BaTiO3) and lead titanate (PbTiO3) using the plane wave pseudopotential method and the PBESOL approximation in the framework of density functional theory (DFT). We use the open-source program Quantum Espresso to perform the calculation of optimized lattice parameters, spontaneous polarization and refractive indices of tetragonal BaTiO3 and PbTiO3 single crystals. In the calculations, all types of pseudopotentials were considered: ultrasoft pseudopotentials (USPP), norm-conserving pseudopotentials (NCPP) and projector-augmented wave (PAW). The calculated values of the optimized lattice parameters, the spontaneous polarization and the refractive indices of the tetragonal BaTiO3 and PbTiO3 show a good agreement with the experimental values, especially for the USPP and PAW pseudopotentials.

A DFT study of optoelectronic, elastic and thermo-electric properties of the double perovskites Rb2SeX6 (X=Br,Cl)

Thermo-electric (TE) material applications reduce reliance on traditional energy resources by converting heat to electric energy. We have studied, for the first time, the thermo-electric properties of Rb2SeX6 (X=Br,Cl). Using norm-conserving pseudo potentials in a plane wave basis set of Quantum Espresso code, the optoelectronic, elastic and thermo-electric properties of Rb2SeX6 (X=Br,Cl) have been investigated using density functional theory. Generalized Gradient Approximation of Perdew Burke Ernzerhof (GGA-PBE) and Generalized Gradient Approximation of Perdew Burke Ernzerhof adapted for Solid (GGA-PBESol) exchange correlation functionals were employed in all calculations. The band structure plots suggest that the studied double perovskites have indirect band gaps. Rb2SeBr6 band gap values of 1.7574/ 1.569 eV (using GGA-PBE/PBEsol) are remarkably similar to that of two effective inorganic/organic perovskites FAPbI3 and MAPbI3 . Maximum peaks generated from refractive index results indicate possible solar cell uses of the materials because they are in the visible and ultraviolet ranges. The results of other optical properties such as absorption coefficients, electron energy loss, conductivity, and reflectivity concludes that Rb2SeX6 (X=Br,Cl) have good values for electron generation, high potential for applications in the optoelectronic industry and are semiconductor in nature. The calculated shear anisotropy values of Rb2SeBr6/Cl6 are 3.09/1.71, suggesting that they are isotropic materials. With calculated Poisson’s ratio of 0.32 and 0.26, the materials are predicted to be ductile in nature. The two materials are appropriate for thermo-electric applications since their thermal to electrical conductivity ratio are small (the order of 10-5). The calculated minimum values of Seebeck coefficient values of 0.198×103 / 0.166 ×103 (mV/K) at 750 K, for Rb2SeBr6/Cl6 are positive, indicating that they have p-type conduction. Figure of merit values at all temperature range considered are greater than one (ZT > 1) for both Rb2SeBr6 and Rb2SeCl6, suggesting that they are good thermo-electric materials. The results of the calculations provide the basis for the industrial application of Rb2SeBr6/Cl6 as solar cells.

Infrared and Raman vibrational modelling of β-C2S and C3S compounds

Density Functional Theory was applied to obtain a refined IR and Raman theoretical vibrational spectra of the major silicate phases of Portland cement, i.e. dicalcium silicate (β-C2S) and tricalcium silicate (C3S). We investigate the crystal structure and vibrational frequency values in the electronic ground state using a plane wave density functional method with new norm-conserving pseudopotentials. The calculated results show that the optimized geometry can reproduce the crystal structure well, and the theoretical vibrational frequency values show good agreement with the experimental values and allow us to describe some uncertainties in the experimental assignments previously described. The assignment performed for 100–1200 cm−1 region shows that the vibrations of C2S and C3S are presented in the IR spectra, while the vibrations with a minor contribution from the silica groups are observed in the Raman spectra. Satellite bands from C3S are observed in the experimental spectra but there not in the calculated one.

Enhanced photocatalytic performance of a stable type-II PtSe2/GaSe van der Waals heterostructure.

In this investigation, the structural, electronic, and optical properties of two-dimensional van der Waals heterostructure (vdwHS) PtSe2/GaSe with three different configurations have been studied using density functional theory with the generalized gradient approximation. All three optimized vdwHSs PtSe2/GaSe have positive phonon frequencies and hexagonal unit cells. The hybrid exchange-correlation functional has been employed to study the electronic properties of vdwHSs PtSe2/GaSe. The vdwHSs PtSe2/GaSe shows semiconducting behavior with indirect Type-II bandgaps, which have been confirmed by the charge density difference, electrostatic potential, work function, and band edge calculations. Additionally, from the band edge positions, the vdwHSs PtSe2/GaSe are analyzed for photocatalytic activities. The optical properties such as extinction coefficient, refractive index, reflectivity, energy loss spectrum, and absorption coefficient have been studied using norm-conserving pseudo-potentials. The vdwHSs PtSe2/GaSe exhibit consistent absorption from the visible to the ultraviolet region of the electromagnetic spectrum. From the obtained results, we conclude that vdwHSs PtSe2/GaSe could be utilized for H2 production through photocatalytic activity as well as for optoelectronic devices and their application.

Greater transferability and accuracy of norm-conserving pseudopotentials using nonlinear core corrections.

We present an investigation into the transferability of pseudopotentials (PPs) with a nonlinear core correction (NLCC) using the Goedecker, Teter, and Hutter (GTH) protocol across a range of pure GGA, meta-GGA and hybrid functionals, and their impact on thermochemical and non-thermochemical properties. The GTH-NLCC PP for the PBE density functional demonstrates remarkable transferability to the PBE0 and ωB97X-V exchange-correlation functionals, and relative to no NLCC, improves agreement significantly for thermochemical benchmarks compared to all-electron calculations. On the other hand, the B97M-rV meta-GGA functional performs poorly with the PBE-derived GTH-NLCC PP, which is mitigated by reoptimizing the NLCC parameters for this specific functional. The findings reveal that atomization energies exhibit the greatest improvements from use of the NLCC, which thus provides an important correction needed for covalent interactions relevant to applications involving chemical reactivity. Finally we test the NLCC-GTH PPs when combined with medium-size TZV2P molecularly optimized (MOLOPT) basis sets which are typically utilized in condensed phase simulations, and show that they lead to consistently good results when compared to all-electron calculations for atomization energies, ionization potentials, barrier heights, and non-covalent interactions, but lead to somewhat larger errors for electron affinities.

Li/Na atoms' substitution effects on the structural, electronic, and mechanical properties of the CaSnO3 perovskite for battery applications

Calcium stannate perovskite CaSnO3 (CSO) is an interesting material for lithium- and sodium-ion batteries (LIBs and SIBs), due to their high hardness and mechanical resistance. In this work, we systematically investigate the effects of substituting lithium/sodium (Li/Na) atoms on the structural, electronic, and mechanical properties of LixCa1-xSnO3 and NaxCa1-xSnO3 perovskite compounds with concentrations × = 0.00, 0.25 and 0.50 using DFT-based in the generalized gradient approximation (GGA) and norm-conserving pseudopotential. This allows describing which substitution generates optimal conditions for structural, mechanical stability, and electronic properties for the CSO perovskite.Lattice parameters behavior as a function of the concentration shows a gradual decrease (increase) in a and c as Li (Na) amount increases. The Na configurations show a small volumetric expansion, on the other hand, the Li configurations have contractions. According to the cohesive energy analysis, it is demonstrated that Li concentrations are more stable than Na ones. The electronic band structure shows that all the Li and Na configurations have a metallic behavior. Finally, the mechanical properties indicate that Na degrades in a less proportion than the Li-substituted compounds; these results demonstrate that SIBs are more favorable than LIBs.