First principles study on the structural stability, mechanical stability and optoelectronic properties of alkali-based single halide perovskite compounds XMgI3 (X = Li/Na): DFT insight

Metal-halide perovskites are recognized as cutting-edge solar energy technology, boasting remarkable absorption capabilities, minimal environmental impact, and cost-effectiveness. This study delves into the structural stability, mechanical stability, and optoelectronic properties of lead-free halide perovskites, specifically XMgI3 (X = Li/Na), by utilizing the CASTEP and WIEN2k software along with the GGA-PBE and Tran–Blaha modified Becke–Johnson (TB-mBJ) exchange–correlation functions to compare their electronic properties. The structural and mechanical stabilities were confirmed by assessing their tolerance factor and formation energy and by evaluating their elastic constants, respectively. Using the TB-mBJ exchange–correlation potential function, the calculated indirect band gap values for LiMgI3 and NaMgI3 were 2.474 and 2.556 eV, respectively. These band gaps are suitable for solar energy harvesting due to their broad optical absorption ranging from infrared to visible light. The partial density of states and the total density of states were determined to investigate the contribution of individual atoms. Consequently, this study can guide researchers focusing on the experimental synthesis of these materials at the laboratory scale for in-depth exploration, particularly in applications such as photovoltaics and various optoelectronic devices.


Introduction
With the increase in the world's population and economy, the demand for energy is growing. 1 However, non-renewable energy sources such as fossil fuels employed to meet this demand are being depleted, 2,3 and the energy generated from non-renewable energy sources is harmful to the surrounding environment due to CO 2 emissions. 4Therefore, it is necessary to discover safe and reliable new energy technologies.2][13] The most promising solar energy in recent years has been organic-inorganic halide perovskite (OIHP)-based PV technology, which is efficient, inexpensive, and has special optoelectronic effects.However, the low stability of OIHPs hinders their application in the production of solar cells.For example, conventional OIHPs such as MAPbI 3 (ref.14) and MAPbX 3 (X = Cl, Br, and I) 15 have attracted considerable interest due to their high efficiency and suitability for optoelectronic applications, despite concerns regarding their instability and the toxicity associated with the presence of lead (Pb).Most of the metal halide Pb-free materials such as ZYbI 3 (Z = Rb, Cs) have been used as perovskite materials as a result of their high efficiency. 16Mahmood et al. 16 reported the preparation of RbYbI 3 and CsYbI 3 with a narrow direct band of 1.22 eV and 1.12 eV, respectively, which are considered suitable for the fabrication of solar cells.
Interestingly, Ray et al. 17 proposed substituting divalent atoms (such as Ge, Sn, Pb, Mg, Ca, Sr, and Ba) for Pb to address the toxicity concerns associated with perovskite materials in photovoltaic applications.They found that Mg and Ba perovskites are unlikely to form in the cubic, tetragonal, and orthorhombic phases due to their positive formation energies.Alternatively, Ca and Sr perovskites exhibit negative formation energies with respect to the metal-iodide precursors and they possess wide band gaps, making them less favorable candidates for use in photovoltaic devices.They additionally reported that the predictive capability of a local density functional with a nonseparable gradient approximation (NGA) closely resembles that of hybrid functionals in terms of band gap predictions, particularly when M in CsMI 3 represents a p-block element such as Pb, Sn, and Ge.
The stability of perovskite materials was identied to be enhanced by introducing K and Rb atoms, [18][19][20][21] specically KGeX 3 and RbGX 3 (where X = Cl or Br).Due to their exceptional optoelectronic properties, which include a high static refractive index, low reectivity, high absorption coefficient, and excellent stability, they have been thoroughly investigated for their possible uses in solar energy and photoelectric systems.Furthermore, they exhibit relatively low levels of toxicity.Additionally, CsGeX 3 (X = F, Cl, and Br) perovskite has been iden-tied as a promising candidate for optoelectronic applications owing to its exceptional electronic, optical, thermoelectric, and tunable band gap properties. 5,22,23Lead-free metal halide compounds XBeCl 3 (X = Ga/Ag) 24,25 and NaXCl 3 (X = Be/Mg) 26 have been investigated for application in optoelectronic devices due to their unique structural, optical, electronic, and elastic properties.GaBeCl 3 and AgBeCl 3 have band gaps of 3.03 eV and 3.25 eV, while NaBeCl 3 and NaMgCl 3 show band gaps of 4.15 eV and 4.16 eV, respectively.An increase in atomic radii corresponds to a decrease in band gap.These materials also exhibit ductility, enhancing their suitability for use in optoelectronic devices.First-principle calculations were used to comprehensively evaluate these properties.
The electrical and optical properties of different cubic leadfree halide perovskites, CsMgX 3 (X = Cl and Br) compounds, were studied. 27Similar to these compounds, the cubic halide perovskite CsBX 3 (B = Sn, Ge, X = I, Br, and Cl) compounds were reported to be suitable for optoelectronic energy devices. 28urthermore, recently narrow band gap perovskite materials such as BaZrS 3 (ref.6) and CsYbI 3 (ref.16) were reported.Owing to their narrow band gap, the reported materials are highly recommended for perovskite solar cell applications due to their wide absorption spectrum, primarily covering the visible and infrared regions.However, their elements are categorized as transition metals and lanthanides, which are computationally expensive and result in the formation of compounds with a relatively high density.In the calculation of the electronic structure of transition metals, standard density functional theory (DFT) underestimates the results.Therefore, it is essential to consider additional methods such as incorporating the Hubbard-U function or hybrid function.
Caid et al. 29 introduced Cs 2 CdZnCl 6 , a lead-free halide double perovskite, via DFT calculations.This compound exhibits mechanical stability in the cubic nonmagnetic phase and has a direct band gap of 1.43 eV.Its optical conductivity was determined to be up to 13.0 eV, showing a signicant ultraviolet photoresponse.These characteristics suggest that Cs 2 CdZnCl 6 is promising for optoelectronic applications, making it a valuable addition to materials science and technology.In another study, Caid et al. 30 extensively investigated Cs 2 B 0 B 00 Br 6 double perovskites, revealing their potential in various technological elds.These compounds demonstrated robust stability with negative formation energies and favorable elastic constants.In the case of both BeMg and CdBe, Cs 2 B 0 B 00 Br 6 exhibited direct band gaps, whereas indirect band gaps were observed with CdGe, GeMg, GeZn, and MgZn.The optical property assessments highlighted the exceptional dielectric characteristics and strong light absorption in visible and UV spectra of these perovskites, making them promising candidates for solar cells and energy harvesting technologies, indicating signicant potential for future energy and optoelectronic advancements.
Therefore, to address these issues, we were motivated to develop transition metal-free and lightweight compounds within alkali and alkaline-based halide perovskite materials (XMgI 3 , X = Li/Na).Our proposed materials exhibit superior optoelectronic properties compared to previous reports, particularly maintaining a narrow band gap, which enhances the photon absorption region.To the authors' knowledge, there are currently no reports available on alkali metal-based iodide perovskite compounds, namely LiMgI 3 and NaMgI 3 , based on either DFT simulations or experimental studies.Thus, this study aimed to thoroughly investigate the mechanical stability, optoelectronic characteristics, and atomic population analysis of XMgI 3 (X = Li/Na) using DFT and the CASTEP and WIEN2k soware.This article presents XMgI 3 (X = Li and Na) perovskite compounds as novel and inventive materials for solar cell applications that outperforms all previously investigated metalhalide perovskites.

Computational analysis
In this investigation, we applied the plane wave pseudopotential approach utilizing CASTEP soware [31][32][33][34] based on DFT principles.The electron exchange-correlation energy was computed following the Perdew-Burke-Ernzerhof (PBE) scheme within the framework of the generalized gradient approximation (GGA). 35We used an ultraso pseudopotential 36 to describe the atomic potential.The electron congurations of the Li, Na, Mg, and I atoms are 1s 2 2s 1 , [Ne] 3s 1 , [Ne] 3s 2 , and [Kr] 5s 2 4d 10 5p 5 , respectively.For optimizing the crystal structure, we employed a plane wave cutoff energy of 700 eV and a 6 × 6 × 6 k-point grid in the Brillouin zone 37 for both XMgI 3 (X = Li/Na) compounds.Structure optimization and lattice parameter calculations were carried out using the BFGS method.The specic tolerances included a total energy of 5.0 × 10 −6 eV per atom, selfconsistent eld (scf) of 5.0 × 10 −7 eV per atom, maximum force of 0.01 eV Å −1 , maximum stress of 0.02 GP, and maximum displacement of 5.0 × 10 −4 Å.The elastic constants, electron localization function, and Mulliken population analysis were performed using CASTEP soware.The Bader change analysis was performed using Quantum ESPRESSO.
To precisely assess the ground-state characteristics and optoelectronic properties of the materials, the full-potential linearized augmented plane wave (FP-LAPW) 38 method within the WIEN2k code was used. 39The Trans-Blaha (TB-mBJ) exchange-correlation approximations were employed.Implementing the TB-mBJ function yielded a highly accurate band gap value. 40,41For the specied objective, the adjusted fundamental parameters were chosen as the standard conguration in the computational soware, incorporating R MT × K max = 8, with K max denoting the wave vector and R MT representing the muffin-tin sphere radius.A k-mesh of 15 × 15 × 15 was chosen

Nanoscale Advances Paper
to ensure precise energy calculations.The self-consistent eld method with an energy convergence threshold of 10 −4 Ry was employed in the iteration process.The comparison of the GGA-PBE function employed in CASTEP and WIEN2k was conducted comprehensively in this study.Moreover, to validate the signicance of FP-LPAW in conjunction with the mBJ hybrid function, an analysis employing the Heyd-Scuseria-Ernzerhof 2006 functional (HSE06) within the framework of Quantum ESPRESSO (QE) was conducted. 42In the process of the scf calculation, a convergence threshold of 60 Ry was established for the E cutoff wave function, and a denser k-mesh of 14 × 14 × 14 in the Brillouin zone was utilized.The phonon dispersion analysis was conducted through the utilization of the phonopy module integrated in the WIEN2k soware.Moreover, ab initio molecular dynamics (AIMD) simulations were executed using the Forcite tools, wherein the NVE thermodynamic ensemble was designated to assess the thermal stability of the compounds under investigation at room temperature.

Structural analysis
The investigated alkali-based single halide perovskite compounds, denoted as XMgI 3 (X = Li/Na), exhibit a cubic ABX 3 structure, conforming to the Pm 3m crystallographic space group (no.221).The cubic primitive unit cell of XMgI 3 (where X = Li/ Na) comprises ve atoms.Specically, the X atoms are positioned at the corners in the 1a Wyckoff position (0, 0, 0), the Mg atom occupies the body center in the 1b Wyckoff position (0.5, 0.5, 0.5), and the I atoms are positioned at the face center in the Wyckoff position (0.5, 0.5, 0.0).The optimized bulk crystal structure of XMgI 3 (X = Li/Na) is given in Fig. 1.Table 1 displays the computed the lattice parameters (a = b = c), band gap (E g ), cell volume (V), tolerance factor (s), and formation energy (E f ) values for the investigated compounds.To further conrm the structural stability, specically the lattice parameters, of the cubic perovskite materials, we performed the optimization utilizing the CASTEP and WIEN2k code tools.The optimized lattice constants for LiMgI 3 and NaMgI 3 were calculated to be 5.8022/5.8297Å and 5.8207/5.8408Å utilizing CASTEP and WIEN2k with the corresponding cell volumes of 194.2461/ 198.0969Å 3 and 195.3316/199.2664Å 3 , respectively, as indicated in Table 1 and Fig. 2A.The lattice parameters increased by approximately 0.47% and 0.34% for LiMgI 3 and NaMgI 3 , respectively, while the volume increased by nearly 2% for both materials.These results indicate that their cubic structures remain stable even when optimized with different tools.Due to ionic differences between the Li and Na atoms, variations in both their lattice constants and cell volume were observed when utilizing the same tool (either CASTEP or WIEN2k).
The Goldschmidt tolerance factor, denoted as "s", is widely used to investigate the structural stability of ABX 3 perovskites, giving insight to determine the compatibility of the cations within their structure, which is dened as follows: where r A , r B , and r X are the A-site, B-site, and X-site ionic radii, respectively.In general, perovskite structures are stable in cubic structures when their tolerance factor is in the range of 0.8 < s < 1. [45][46][47][48] Alternatively, when s < 0.8, the perovskite adopts an orthorhombic structure, whereas for s > 1, it is categorized as having a hexagonal crystal structure. 49The ionic radii used in this calculation were 1.13, 1.39, 0.72, and 2.20 Å for Li, Na, Mg, and I, respectively. 47Thus, the calculated tolerance factors of 0.81 for LiMgI 3 and 0.87 for NaMgI 3 fall within the acceptable range, as listed in Table 1, resulting in the conrmation of the ideal cubic structural stability of the investigated compounds.Furthermore, the formation energy (DE f ) of the compounds was calculated to conrm their thermodynamic stability at room temperature using eqn (2), as follows: 50 where E tot (XMgI 3 ) represents the total energy of the XMgI 3 compound, N is the number of atoms, and E i (X), E i (Mg), and E i (I) are the total energy of the individual atoms, including X (where X = Li and Na), Mg, and I, respectively.It is evident that the obtained negative values of the formation energy in Table 1 indicate the stability of the proposed perovskite compounds under ambient conditions.This information is valuable for experimentalists considering the synthesis of these materials at the laboratory scale.To the best of our knowledge, these materials have not been experimentally reported.Nevertheless, the comparable works reported using DFT, focusing on the lattice parameters and band gap, are compiled in Table 1 to validate the credibility of our study.Furthermore, phonon dispersion calculations offers valuable insights into the dynamical stability of the materials being studied. 51The phonon dispersion curves were obtained by computing the phonon frequencies along the high symmetry directions (X / G / R / M / G) of the Brillouin zone, as illustrated in Fig. 2B.Typically, stable phonon structures exhibit real frequencies, as denoted by positive signs, whereas unstable structures are characterized by imaginary frequencies, which are indicated by negative signs.The computed phonon spectra revealed that both the LiMgI 3 and NaMgI 3 perovskite materials exhibit real frequencies at the gamma (G) point, conrming their dynamic stability.However, at the X-point in the rst Brillouin zone, a less negative phonon frequency (∼−0.54THz) was observed, suggesting localized instabilities or specic vibrational modes affecting certain directions or regions within the crystal lattice. 52,53Positive phonon frequencies at the G-point indicate stable vibrational modes at the center of the Brillouin zone, typically corresponding to the collective motions of the entire crystal lattice or large-scale atomic displacements.These stable modes contribute to the overall stability of the material by ensuring well-dened lattice vibrations without leading to global instabilities.
To evaluate the thermal stability of LiMgI 3 and NaMgI 3 , ab initio molecular dynamics simulations (AIMD) were carried out at room temperature over a period of 500 ps, utilizing a time step of 1 fs.The results, as depicted in Fig. 2C, reveal that minimal energy uctuations occurred during the simulation.Furthermore, both perovskite frameworks preserved their original congurations without exhibiting any discernible structural deformities.Thus, these ndings undeniably conrm the stability of LiMgI 3 and NaMgI 3 at room temperature. 26

Atomic populations analysis
Mulliken atomic population (MAP) is suitable for gaining insight into the atomic bonding (such as chemical bonds, bond overlap, Table 1 The computed optimized lattice constants (a = b = c), band gap (E g ) using CASTEP, WIEN2k, and QE, cell volume (V), formation of energy (E f ), and tolerance factor (s) for the cubic structure of XMgI 3 (X = Li/Na) halide single perovskite compounds and some reported data

Nanoscale Advances Paper
and movement of charge between atoms) in compounds. 54,55The calculated bond populations, charges, and bond lengths are summarized in Table 2 for the XMgI 3 (X = Li/Na) compounds.The determined positive (Li, Mg, and Na) and negative (I) charges in the XMgI 3 (X = Li/Na) compounds, as shown in Table 2, suggest that charge is transferred from the alkali metals (Li and Na) and alkaline earth metallic (Mg) atoms to the halogen (I).In both the LiMgI 3 and NaMgI 3 compounds, charge transfer of 0.46 e and 0.28 e to iodine (I), respectively, was observed.Moreover, we conducted theoretical investigations into the charge-transfer mechanism utilizing Bader analysis, 56 as depicted in Table 2.The Bader charge analysis utilizes electron density to determine the electronic charges associated with individual atoms. 57It helps in understanding the nature of chemical bonding and intermolecular interactions.The ability of the I atom to accept electrons from Li/ Na and Mg suggests its ionization nature.Substituting Li with Na resulted in slight alterations in the Bader charges, as shown in Table 2.In the LiMgI 3 structure, the average Bader charge on the I atoms was −0.816 e, while that in the NaMgI 3 structure was −0.784 e.The variation observed in the charge distribution computed using Mulliken and Bader analyses likely stems from their differing methodologies.Bader analysis emphasizes the topological aspects of electron density, whereas Mulliken analysis relies on population analysis.Additionally, Table 2 provides data on the bond populations, population ionicity, and bond lengths for the Li-I, Mg-I, Na-I, and I-I bonds within the XMgI 3 (X = Li/Na) perovskite structure.Positive and negative values signify bonding and antibonding states, respectively, with a higher positive bond population denoting increased covalency in the bond. 58For example, the lower positive bond population in Li/Na-I (associated with longer bond lengths) and Mg-I (associated with shorter bond lengths) in Table 2 suggests a high level of ionic character and insignicant covalent bonding in both compounds.Additionally, the negative bond population of I-I indicates the formation of antibonding interactions.

Electron localization function
The electron localization function (ELF) serves as a valuable tool for analyzing the electron distribution and localization during the formation of a crystalline structure. 59ELF quanties the probability of nding electrons with identical spins and xed positions within a reference electron eld, offering insights into the electron localization within the lattice space and revealing the bonding characteristics in materials.Fig. 3 shows the 3D and 110 crystal dimension ELF maps for the XMgI 3 compounds (X = Li/Na).The right side of the ELF maps illustrates the ELF scaling, including the coordination scaling, where regions with lower electron density are depicted in blue and that with higher electron density in red.The examination of the ELF maps for the XMgI 3 compounds revealed a heightened degree of electron localization around the I ions, which is attributed to their signicant electronegativity.This observation is consistent with the ndings from the Mulliken population analysis and the Bader charge analysis.

Elastic constants and mechanical properties
The mechanical and structural stability of materials can be determined by measuring their elastic constant (C ij ).These measurements indicate that materials deform under stress but return to their original shape once the stress is removed. 60Thus, to ensure the materials' mechanical stability, we calculated C ij using the nite strain theory, providing insights into their strength and response to external forces.In this investigation, we rst performed elastic constant calculations, followed by determining the mechanical and structural stability   3. The mechanical properties of these materials include hardness, soness, compressibility, ductility, brittleness, and strength, which are vital for industrial applications.In cubical crystals, there are three independent elastic constants, C 11 , C 12 , and C 44 , 61 which were calculated using the nite strain theory. 62According to the symmetry results, C 11 = C 22 = C 33 , C 12 = C 23 = C 13 , and C 44 = C 55 = C 66 . 635][66] In addition, the elastic constant is used to determine the mechanical stability of the compounds according to the Born stability criterion, 67 where the crystal deformations in a compound depends on a specic criterion. 68In the case of the cubic phase, three independent elastic constants are required to fulll the Born stability criterion, as given in eqn (3), to determine the stability of a compound: 69 The elasticity of length and shape is related to the C 11 and C 12 constants, respectively.For a cubic structure, C 11 must be considerably larger than C 12 . 61Interestingly, the estimated elastic constants for both the LiMgI 3 and NaMgI 3 compounds listed in Table 3 conrm they satisfy the Born stability criteria, indicating their mechanical stability.The Pugh ratio (B/G) and Cauchy pressure are useful parameters to assess a material's brittleness and ductility. 60In this case, a material is determined to be ductile when B/G is equal to or greater than 1.75, whereas when it is less than 1.75, the material is brittle. 70According to this study, the (B/G) ratios for LiMgI 3 (2.1) and NaMgI 3 (2.03) are higher than 1.75, suggesting that these materials exhibit ductile behavior, as shown in Table 3.
The Cauchy pressure for the cubic LiMgI 3 and NaMgI 3 determines their ductility and brittleness based on the negative and positive values of C 12 − C 44 , as previously reported. 71A positive value indicates ductility, whereas a negative value suggests brittleness.In this study, the difference between C 12 and C 44 (Table 3) for both compounds are positive values, resulting in indicating ductility.Alternatively, the Voigt-Reuss-Hill (VRH) method is an approach used to calculate the B, G, E, and s of the materials.The Voigt and Reuss methods are used to determine the bulk 72 and shear 73 moduli of materials, respectively, while Hill's 64 method is used to calculate the average values of these properties.The Reuss bulk modulus (B R ) and Voigt bulk modulus (B V ) are given by eqn (4): The shear modulus (G V ) of a material according to Voigt 72 is presented in eqn ( 5): moreover, the shear modulus (G R ) was calculated by Reuss using eqn ( 6): and their average value was computed by Hill, 64 as indicated in eqn ( 7) and ( 8), respectively: The values for B, G, and E are presented in Table 3.The B (25.08, 25.08), G (11.93, 12.36), and E (30.88, 31.85) of LiMgI 3 and NaMgI 3 , respectively, indicate that both materials exhibit comparable resistance to changes in volume and shape.Moreover, the observed E and s listed in Table 3 can be calculated using the following expression, respectively: 69 The resistance of a material, which is indicated by its stiffness, is quantied by E, where materials with high E values exhibit increased resistance.Therefore, as shown in Table 3, LiMgI 3 and NaMgI 3 exhibit identical and high E values, suggesting that they possess a signicant level of stiffness.Alternatively, s, which is limited to values between 0.25 and 0.50, determines the behavior of bond forces.A s value greater than 0.26 indicates that the material is ductile, whereas a value less than 0.26 indicates that it is brittle. 74Hence, as indicated in Table 3 The calculated elastic constant C ij (GPa), bulk modulus B (GPa), shear modulus G (GPa), Young's modulus E (GPa), Poisson's ratio (s), Pugh ratio (B/G), Cauchy pressure (C P = C 12 − C 44 ) (GPa), anisotropy factor (A), sound velocity (V), Debye temperature q D (K), and melting temperature T m (K) of XMgI 3 (X = Li/Na) single perovskite compounds Table 3, the bond force of the LiMgI 3 and NaMgI 3 compounds is close to the center of the above-mentioned range (0.29).As a result, both compounds exhibit a characteristic of ductility that agrees well with the results of the B/G ratio and Cauchy's pressure equation.Moreover, the elastic anisotropy factor (A), which is used to determine the direction-dependent properties of a system, was computed by using eqn (11): 69

Material property
For isotropic compounds, A is ideally 1, indicating uniformity, whereas any deviation from 1 indicates the presence of anisotropy. 75,76As shown in Table 3, the computed anisotropy factors for both compounds deviated from unity (0.81 for LiMgI 3 and 0.84 for NaMgI 3 ), indicating the presence of anisotropy in these materials.
The calculation of G and B enables the determination of the transverse sound velocity, v t (eqn ( 12)), and longitudinal sound velocity, v l (eqn ( 13)), of XMgI 3 (X = Li/Na) crystals, where r refers to the density of the investigated compounds as follows: [77][78][79][80] and Moreover, utilizing eqn (10) and (11), eqn (12) allows the calculation of the mean sound velocity (v m ), as follows: The Debye temperature (q D ) is used to determine the thermal characteristics of a material including its melting temperature and thermal expansion.Employing eqn (14), q D of LiMgI 3 and NaMgI 3 can be obtained, as expressed in eqn (15): where h represents the Planck constant, k represents the Boltzmann constant, N A denotes the Avogadro constant, M signies the molecular mass, and n indicates the number of atoms in the compound.This calculation yielded valuable insights into the thermal properties of the investigated crystals.
The calculated values (v t , v l , and v m ) and q D for the LiMgI 3 and NaMgI 3 compounds are shown in Table 3.The calculated parameters for LiMgI 3 and NaMgI 3 compounds approximately show similar q D , making both compounds have similar thermal conductivity and melting temperature.The melting temperature (T m ) is a crucial parameter for determining the bond characteristics and thermal conductivity of XMgI 3 (X = Li/Na).The melting temperatures of these cubic materials, as listed in Table 3, were calculated using eqn (16): [81][82][83] T m = 553 + 5.91C 11 (16)   As illustrated in Table 3, both LiMgI 3 and NaMgI 3 exhibit high and comparable melting temperatures (T m ), indicating their suitability for high-temperature optoelectronic applications.Thus, the high computed melting temperatures for these compounds imply their strong bonding characteristics and excellent thermal conductivity.

Electronic properties
In this section, we present the analysis of the energy dispersion and electronic states of the compounds, which determine their electronic properties.The electronic bandgap (E g ) of the materials was obtained through band structure analysis, while their partial density of states (PDOS) and total density of states (TDOS) were determined, which describe the distribution of electronics at different energy levels.The electronic properties of the XMgI 3 (X = Li/Na) single perovskite materials were investigated using the CASTEP and WIEN2k codes.Fig. 4 and Table 1 demonstrate the narrow indirect band gaps of 1.174/ 1.176 eV, 1.181/1.185eV, 2.474/2.556eV, and 2.219/2.239eV for LiMgI 3 /NaMgI 3 , which were calculated employing PBE-CASTEP, PBE-WIEN2k, mBJ-WIEN2k, and HSE06-QE, respectively.It is clear that the standard DFT PBE function consistently underestimated the electronic structure, particularly the band gap, in both CASTEP and WIEN2k, as shown in Fig. 4a and b, respectively.However, despite this underestimation, the results suggest that the proposed compounds exhibit sufficient stability, given that they yielded similar outcomes across different computational codes.Furthermore, the mBJ correlation function (Fig. 4c) accurately predicted the band gap, which correlates well with the experimental values.This suggests that LiMgI 3 and NaMgI 3 have a favorable band gap, which shows their potential for use as semiconductor materials in PV devices and other optoelectronic applications.Further examination was conducted to validate the prediction by the mBJ correlation function employing the alternative robust HSE06 hybrid functional (Fig. 4d), leading to the sustained accurate correlation with the mBJ prediction.The stability and reliability of both LiMgI 3 and NaMgI 3 are evident from the predictions by the mBJ and HSE06 correlation functions, suggesting their suitability for future empirical investigations.
Fig. 4a-d show that both the conduction band minimum (CBM) and valence band maximum (VBM) of the LiMgI 3 and NaMgI 3 compounds are located at G and R, respectively, in the rst Brillion zone, which agrees with the previous nding. 84The s-block VB metals are at between the R and M points, suggesting the presence of heavy holes, whereas the p-block metals have a greater dispersed CB at the G-point. 85ig. 5 and 6 show the analysis of the partial density of states (PDOS) and total partial density of states (TDOS) of the LiMgI 3 and NaMgI 3 compounds, respectively.In these gures, the TDOS and PDOS were analyzed to identify their effects on the VB

Paper
Nanoscale Advances to CB electronic characteristics.Fig. 5a illustrates the dominance of the Li p-state in the VB due to its high intensity near to the Fermi level.The proximity of the peak intensity to the Fermi level is crucial for electron transitions from the VB to the CB.Therefore, the Li p-state is the dominant peak in the electron transfer from VB to CB.In the case of the Mg atom, the p-state contributes to the VB, while the s-state contributes to the CB, as depicted in Fig. 5b.Moreover, Fig. 5c clearly shows that the prominent peak near the Fermi level within the VB primarily originates from the I p-state.In summary, the TDOS in Fig. 5d shows that the VB is mainly contributed by the p-states of Li, Mg, and I, while the s-states of Li and Mg contribute to the CB with a smaller contribution from the I p-states in the LiMgI 3 compound.Notably, when replacing the Li atom with Na, similar trends were observed, as depicted in Fig. 6.These similarities arise from the shared properties of Li and Na as elements within the same family.Our investigation is supported by the previous reports on similar families. 20,866 Optical properties Details on the optical properties of perovskite compounds can be obtained by calculating various parameters, including the dielectric constant, absorption coefficient, extinction coefficient, refractive index, reectance, transmission, and loss function.The expression for the complex dielectric function, 3(u), with respect to angular frequency is provided as follows: [77][78][79][80] where 3 1 (u) and 3 2 (u) represent the real and imaginary components of the complex dielectric function, respectively, and u is the angular frequency of the incident photon.The absorption coefficient, a(u), extinction coefficient, k(u), refractive index, n(u), reectivity, R(u), conductivity, s(u), and loss function, L(u), were calculated based on eqn (17), as follows: and Fig. 7 shows the optical properties of the investigated LiMgI 3 and NaMgI 3 compounds.Based on the photon radiation on the surface of the compounds, the dielectric function determines their optical properties.Fig. 7a and b show the dielectric functions of the real and imaginary parts of XMgI 3 (X = Li/Na), respectively.As shown in Fig. 7a, according to the static frequency of 3 1 (0), the real dielectric constant is 3.26 for LiMgI 3 and 3.19 for NaMgI 3 .The 3 1 (u) of LiMgI 3 and NaMgI 3 reached the maximum peak intensities of 3.96 and 3.86 at the 3.10 eV edge of the visible light region, respectively.However, as shown in Fig. 7b, LiMgI 3 and NaMgI 3 revealed 0 eV at 3 2 (0), indicating the absence of dissipated energy.The imaginary part of a dielectric material, 3 2 (u), increases until energy dissipation occurs in the photon energy absorption region.Thus, the imaginary parts of a dielectric material reect the behavior of its absorption coefficient.
Fig. 7c shows the absorption coefficient of the cubic structure of LiMgI 3 and NaMgI 3 .In Fig. 7c, it is evident that both the LiMgI 3 and NaMgI 3 compounds initiate photon absorption at 1.3 eV and display a wide absorption peak spanning the infrared (IR) to ultraviolet (UV) region.Therefore, the proposed materials are highly suitable as novel lead-free halide perovskite materials for solar energy harvesting applications.The extinction coefficient demonstrates how the light intensity attenuated on the investigated compounds, LiMgI 3 and NaMgI 3 , as displayed in Fig. 7d, where both compounds exhibit similar behavior.Thus, the light-absorbing properties of both compounds are similar and consistent with the discussion in Fig. 7c.
Fig. 7e shows the evaluation of the optical characteristics of the refractive index, which explains the bending or refracting behavior of light in a medium.The static refractive indices of LiMgI 3 and NaMgI 3 were determined to be 1.81 and 1.78, respectively.This suggests that the square of the refractive index (n(u)) 2 corresponds to 3 1 (u) at a static frequency, as represented in Fig. 7a.Additionally, both compounds reached their maximum at the edge of higher energy of visible light (3.1 eV).The detection of the refractive index peaks within the low-energy segment of the spectrum is a crucial aspect, especially in the visible region, providing valuable insights for applications in optoelectronic devices.Thus, examining the reectivity of perovskites offers insights into their surface attributes.Fig. 7f illustrates the reectance spectra of the LiMgI 3 and NaMgI 3 perovskite compounds as a function of photon energy.Both compounds consistently demonstrate low reectivity (less than 0.10) across the infrared to visible region spectrum, which is negligible.Remarkably, their reectivity remains consistently low in the visible to infrared range, making these materials well-suited for solar cell applications.
Fig. 7g shows how electromagnetic (EM) radiation induces optical excitation in the compounds, displaying the conductivity s(u) of XMgI 3 (X = Li and Na).In the case of LiMgI 3 and NaMgI 3 , their s(u) remains zero in the photon energy range of 0 eV to 0.98 eV due to the absence of optical excitation, indicating the absence of interaction between incident EM radiation and the electrons in these materials in this energy range.As LiMgI 3 and NaMgI 3 absorb energy, their electrons shi from the VB to the CB.The consistency of their spectra with the absorption coefficient analysis (Fig. 7c) suggests that absorbing light within the visible part of the electromagnetic spectrum is particularly benecial for solar cell applications.Another signicant consideration involves measuring the energy loss function, denoted as L(u), which accounts for inter-band, intra-band, and plasmonic interactions.Additionally, the optical loss serves as an indicator of energy dissipation through scattering, dispersion, and thermal effects, as depicted in Fig. 7h.Fascinatingly, for both LiMgI 3 and NaMgI 3 , their optical loss in the visible region (1.55-3.10eV) is nearly insignicant compared to that in the ultraviolet region (3.10-4.0eV), conrming the suitability of these materials to realize an excellent optoelectronic performance within the visible range of the electromagnetic spectrum.

Conclusion
This investigation explored the structural, electronic, optical, and mechanical attributes of cubic XMgI 3 (X = Li/Na) through rst-principle calculations by employing CASTEP and the GGA-PBE exchange-correlation function.Our examination of the elastic constants and mechanical properties indicated the notable stability and ductile nature of each perovskite material.Additionally, we investigated the electron charge distribution using Mulliken atomic population analysis and the electron localization function.The developed compounds exhibited advantageous indirect narrow band gaps (2.474 eV for LiMgI 3 and 2.556 eV for NaMgI 3 ) and a broad-spectrum absorption coefficient using the mBJ correlation function in WIEN2k.Based on our ndings, especially the band gap and optical properties, these alkali-based halide perovskite materials are suitable for several applications, such as photovoltaics, photocatalysts, and optoelectronic devices.DFT analysis suggested their highly favorable band gap for photovoltaic applications.Consequently, further efforts are required for the synthesis and exploration of these alkali-based halide perovskite materials for realistic applications and future advancements.

Fig. 6
Fig. 6 Calculated PDOS of individual atoms in the NaMgI 3 compound including (a) Na, (b) Mg, and (c) I and (d) TDOS of the NaMgI 3 perovskite structure.