SciELO - Scientific Electronic Library Online

 
vol.69 número5Exact Green’s functions for localized irreversible potentialsAlternative approach to explore the stability of floating bodies índice de autoresíndice de materiabúsqueda de artículos
Home Pagelista alfabética de revistas  

Servicios Personalizados

Revista

Articulo

Indicadores

Links relacionados

  • No hay artículos similaresSimilares en SciELO

Compartir


Revista mexicana de física

versión impresa ISSN 0035-001X

Rev. mex. fis. vol.69 no.5 México sep./oct. 2023  Epub 28-Oct-2024

https://doi.org/10.31349/revmexfis.69.050501 

Condensed Matter

Structural, mechanical, electronic, thermal, and optical properties of the inverse-Heusler compounds X2RuPb(X = La, Sc): A first-principles investigation

A. Bouazzaa  * 

M. Khiratb 

M. Larbic 

N. Bettahar

D. Rached

a L2GEGI Laboratory, University of Tiaret, 14000 Tiaret, Algeria.

bMagnetic Materials Laboratory Djillali Liabes University of Sidi Bel-Abbes, 22000 Sidi Bel-Abbes, Algeria.

c Physics Department, Djillali Liabes University of Sidi Bel-Abbes, 22000 Sidi Bel-Abbes, Algeria.


Abstract

Topological insulators are novel quantum material states with insulating bulk band gaps and topologically protected metallic surface states that have been extensively studied owing to their intriguing properties for spintronic and quantum-computing applications. The structural, mechanical, electronic, thermal, and optical properties of inverse Heusler compounds La2 RuPb and, Sc2 RuPb in two Hg2 CuTi, Cu2 MnAl-type structures were calculated using the full potential linear muffin-tin orbital simulation methodology as implemented in the computer code,which is based on density functional theory.We employed the local-density approximation for the exchange and correlation potential (XC) terms. Consequently, the optical characteristics of La2 RuPb, Sc2 RuPb and elastic constants C ij and their corresponding elastic moduli were computed for the first time. According to our structural calculations, La2 RuPb is more stable in its Hg2 CuTi-type structure than Sc2 RuPb in its Cu2 MnAl-type structure. However, the mechanical characteristics demonstrate their stability in the final stages of elastic deformation.

Keywords: Inverse-Heusler compounds; density functional theory; topological insulators; mechanical properties; optical properties; elastic constants

1. Introduction

Since Heusler compounds are likely to be fresh prospects for future uses, they have rekindled the attention of researchers. Intermetallic Heusler compounds are ternary ferromagnetic compounds with the general formula X2YZ, which were initially identified in 1903 by German mining engineer Friedrich Heusler while examining the Cu2MnSb alloy [1, 2].

Based on their estimated electronic band structures, a novel class of Heusler alloys was anticipated to be labeled “multifunctional topological insulation (TIs)”. Surface states are topologically safeguarded safeguarded against impurity fusion in this new state of matter [3, 4].

Because of their diverse properties for spintronic applications [5, 6], optoelectronics [7], superconductivity [8], shape-memory [9], giant magnetoresistance spin-valve (GMR) [10], thermoelectric applications [11], and spin injection to semiconductors [12], interest in them is steadily growing as in the thin film deposition technology [13-15].

Several theoretical and experimental attempts have been made to anticipate novel topological insulators (TIs) due to their enormous promise in spintronics and quantum computing [16, 17]. (TIs) are a relatively novel material class. Their narrative begins with Bernevig and König’s [18, 19] initial prediction of (HgTe) as a material (TI). Following this finding, (TIs) materials in binary compounds of (II-VI) group elements have been thoroughly explored.

Full-Heusler alloys are ternary intermetallic compounds with the stoichiometric composition X2YZ [1, 20], where X and Y are transition metal elements and Z is a (III-, IV-, or V-group element). These materials have a cubic crystal shape. Full-Heusler compounds have been found in either Hg2CuTi-form, designated X-form F 4¯3 m symmetry (space group N°216) [21], or the (Cu2MnAl)-form denoted as the L21-type). When the (X) component’s valence electron number is greater than the (Y) ones, the compound adopts the L21-type structure [20, 21], but the X-type structure is preferred in the opposite situation [22]. The unit cell of the L21-type shape belonging to the FmFm symmetry (space group N°225) [12, 23] is made up of four interpenetrating (FCC) (Face Centered Cubic sublattices) with the following Wyckoff coordinates: X 1 = (0, 0, 0), X2=1/2,1/2,1/2,  Y=1/4,1/4,1/4, Z=3/4,3/4,3/4. The occupancy parameters of the X2 and Y components are swapped in the (X-type) structure. Inverse Heusler compounds are complete Heusler compounds with a later prototype architecture.

If the two X particles seem no longer equivalent, ternary and quaternary Full-Heusler alloys are formed. Quaternary Heusler alloys with a stoichiometric composition of (XX’YZ) commonly crystallize as the group N°216 [24-26]. The following Wyckoff coordinates X=1/2,1/2,1/2, X'=(0,0,0),Y=1/4,1/4,1/4, Z=3/4,3/4,3/4 are occupied in the unit cell for the (Y-type) structure. The search for (TIs) has recently expanded to include ternary Heusler and chalcopyrite compounds [3, 27]. Zhang et al. revealed that most inverse Heusler compounds show an inverted band order, a promising characteristic of topological insulators. However, their topological insulating state is reported to be sensitive to uniaxial strain and variations in the lattice constants [28].

Understanding a material’s physical qualities is essential for understanding its features in terms of its operational uses [29-32]. To better understand the behavior of the La2RuPb,Sc2RuPb compounds, we present our calculations on the structural, mechanical, electrical, and optical characteristics of the compounds in their phases Hg2CuTi2,Cu2MnAl-Type, respectively.

These analyses are used to calculate the phase transitions and other essential ground state features of these structures. In addition, electronic band structures are also established to comprehend the properties of the (TIs). Our calculations were carried out in the context of density functional theory (DFT) utilizing the local density approximation (LDA), employing the full-potential linear muffin-tin orbital (FPLMTO) approach as implemented in the computer code (LmtART).

To our knowledge, no theoretical or practical study on the electronic states and characteristics of these Heusler alloys has been conducted.

The following is how the remainder of the paper is organized: Section 2 contains computational details and the calculation technique, Sec. 3 has the findings of the structural, mechanical, electronic, and optical characteristics, as well as discussion, and Sec. 4 concludes with a short conclusion.

2. Computational details

Several approaches may be used to determine the physical characteristics of materials, such as electronic structures. The technique of augmented plane waves (FP-LAPW) [33, 34] and the full-potential linear muffin-tin orbital (FP-LMTO) [35, 36] are two of the most well-known ab-initio approaches that are applicable. Moreover, it may be used to solve the (DFT) equations [37], the most widely used and effective method for many years. It also has the benefit of addressing numerous atoms. The latter is now in the tens of thousands and might rise to the hundreds of thousands in the coming years. The calculations given here were performed using the Perdewe-Wang (P.W) [38] model and the (FP-LMTO) technique built inside the (DFT) as incorporated in the (LmtART MStudio MindLab 5-code) produced by S.Y Savrasov [39].

A portion of the overall energy calculations for the potential of exchange and correlation (XC) (energy/potential) was done using the local density approximation (LDA) [40]. Compared to prior versions of the (LMTO) approach, this version of (FP-LMTO) incorporates significant enhancements. The charge density/crystal potential is not subjected to any shape approximation. Furthermore, in this computational strategy, known as (FP-LMTO), the crystal unit cell is partitioned into the muffin-tin (MT) sphere region and the interstitial region (IR) for mathematical convenience, and both zones are handled on the same footing, increasing the accuracy of the eigenvalue results. Furthermore, a complete foundation is employed in this technique. For the spherical component of the potential, fundamental functions are defined in numerical solutions of the radial Schrödinger equation multiplied by spherical harmonics, while basis functions are given in terms of Fourier series for the (IR) regions. Charge density and the related potential were presented within the (MT) spheres, up to lmax=6, using spherical harmonics to carry out calculations for the (La2RuPb) and (Sc2RuPb) compounds. The k-points integration was carried out using the tetrahedron method [41] over the irreducible edge of the first Brillouin region using a (K-point) grid of (10, 10, 10). At the same time, the criterion of convergence up to (10-5) was met self-consistently for the accurate calculations of total energy and charge density.

Table I RMT (a.u.) of Ru, Pb, and X atoms for X2 RuPb compounds. 

Atom La2 RuPb Sc2 RuPb
X 2.929 2.696
Ru 2.901 2.807
Pb 2.956 2.807

3. Results and discussion

3.1. Phase transition and structural properties

In this part, we estimated the structural characteristics of the two full Heusler alloys La2RuPb,Sc2RuPb and looked at both structural phases Hg2CuTi,Cu2MnAl by utilizing the total energy minimization approach (at T = 0°K) for the determination of the equilibrium lattice constant. Calculating the total energy as a volume function was used to optimize the unit cell volume V. Finally, the change of total energy ETOT with volume V is fitted using Murnaghan’s equation of state [42, 43] to yield the equilibrium lattice constant and the bulk modulus:

E(V)=E0+B0VB'VV0B'1B'-1+1-B0V0B'-1, (1)

where E 0 is the ground-state total energy, V 0 is the optimal unit-cell volume, and B 0 and B’ are the bulk modulus and its pressure derivative.

Table II lists our results and additional information for the optimized lattice constants for both compounds.

Table II Equilibrium lattice constants a 0(A°) and total energy E (Ry/unit cell) for (La RuPb, Sc2 RuPb) compounds. 

Compound Calculation a 0(A°) Energy E (Ry/unit cell)
7.3760 [28] -84911.718
Other (Theory) 7.3417 [29] -84911.720
7.3470 [30] -
La2 RuPb 7.1597 [LDA] X- Type -84854.71837
Our calculation 7.5746 [GGA] L2 - Type -84854.7606
7.3962 [LDA] L2 21 - Type -84911.68004
7,338 [GGA] X - Type -84911.6587
Experimental - -
6.7460 [28] -53978.04810
Other (Theory) 6.7531 [29] -53978.02600
Sc2 RuPb 6.6705 [30] -
6.6342 [LDA] X- Type -53939.81095
Our calculation 6.7263 [LDA] L2 21 - Type -53939.82301
6.7934 [GGA] X - Type -53978.00554
6.7237 [GGA] L2 21 - Type -53939.82461
Experimental - -

Since we are not aware of any experimental data, we compare our results with those obtained in the already-published papers of Y. Lakred and co-workers [28], A. Bahlouli and co-workers [29] for the inverse-Heusler alloys X 2RuPb, and X.M. Zhang and co-authors [30] for the Y 2RuPb investigations. Our computed values are in good correlation and agreement with those obtained by the researchers mentioned above (see Table II), which confirm the precision of our obtained results.

The theoretical values supplied in the table using the (GGA) approximation match with our estimated lattice constants for La2RuPb, Sc2RuPb. Unfortunately, to our knowledge, no experimental data for these compounds have been published.

The variation of total energy as a function of volume for the La2RuPb,Sc2RuPb compounds in the Hg2CuTi,Cu2MnAl-type structures are shown in Fig. 1. Table III shows the obtained values for the a=4V03, E0,V0,B' and B 0 parameters. To the best of our knowledge, there is no data in the scientific literature that can be compared to the findings we’ve achieved.

Figure 1 The variation of the total energy as a function of the volume for the compounds La2RuPb, Sc2RuPb in the Hg2CuTi, Cu2MnAl-type structures respectively. 

Table III The calculated equilibrium lattice parameter a 0(A°), primitive-cell volume V 0(A°)3, bulk modulus B 0 GPa, its pressure derivative B’ and the equilibrium energy E 0 Ryd for the La2RuPb, Sc2RuPb compound in the Hg2CuTi, Cu2MnAl-type structures. 

System structure-type a V 0 B 0 B E 0
La2RuPb Hg2 CuTi 7.1597 91.7538 93.4167 4.72155 -84854.7606
Sc2RuPb Cu2MnAl 6.7263 76.0796 127.6952 4.05872 -53939.82301

As shown in Fig. 1, the Hg2CuTi-type structure is found to be the most stable phase for the La2RuPb compound at zero pressure because the crossing of the E(V) curve of the Hg2CuTi-phase is before the curve of the Cu2MnAl-phase for the lowest energy values, while the Sc2RuPb compound is favored in the reverse case Cu2MnAl-phase.

We believe that the values we discovered and described in the table above are consistent with those that will be found in the future by others when studying the two compounds in this or other ways, even though there was no theory or experimental values to compare them with our knowledge. The calculations of the structural properties of the two structures phases are reliable. These values encourage future experimental work to collaborate with our calculated results.

3.2. Elastic properties

The study of the mechanical properties of materials is significant for practical applications. The elastic properties of a solid are determined by the abovementioned elastic constants, which link the mechanical and dynamic behavior of a crystal and provide crucial information about the nature of the forces acting on the solid. The elastic constants, in particular, are macroscopic quantities that describe the stability, stiffness, and anisotropy of materials in inhomogeneous solids.

Understanding interatomic interactions are possible because of the influence of hydrostatic pressure. Mechanical stability and phase transition mechanisms stress to deformations allow for the evaluation of the wave’s elastic straps’ speed of propagation, as well as the assessment of mechanical properties such as ductility/brittleness, specific heat, hardness, melting point, Debye temperature, and thermal expansion coefficient [44]. Furthermore, elastic constants may be used to calculate the bulk modulus, shear modulus, Young’s modulus, and Poisson’s ratio, which are all important in calculating the strength of materials [45].

By studying the impact of hydrostatic pressure, we may learn about interatomic interactions, phase transition, and mechanical stability processes. Stresses to deformations may be used to measure mechanical qualities such as ductility/brittleness, specific heat, hardness, melting point, Debye temperature, and thermal expansion coefficient, as well as the speed of propagation of the wave’s elastic straps [44]. Furthermore, elastic constants may be used to calculate the bulk modulus, shear modulus, Young’s modulus, and Poisson’s ratio, which are all critical in calculating the strength of materials [45].

Using stress-strain equations, the second-order elastic constants (C ij) for the La2RuPb, Sc2RuPb compounds in the Hg2CuTi,Cu2MnAl-type respectively, structure were calculated in this work. Total energy calculations are performed with regard to volume conserving stresses using the Mehl approach [46] as described in Ref. [47] to obtain the (C ij) values. Only three elastic constants C11,C12 and (C 44) are required to demonstrate the mechanical behavior of a cubic crystal and to satisfy the Born Huang criteria for positive assuredness of the stiffness matrix [48]. C44>0,C11>C12,C11+2C12>0 are wellknown requirements for cubic crystal mechanical stability [49]. We utilized the averaged values of the Voigt approximation (the Voigt shear modulus G v ) [50] and the Reuss approximation (the Reuss shear modulus G r ) [51] to derive the bulk structural parameters.

For investigation of the material hardness, the elastic properties, e.g., bulk modulus B, Voigt-Reuss-Hill shear modulus GvRH [52] and the elastic Zener anisotropy factor A [53] are given by the following formulas for cubic crystal:

Bv=Br=C11+2C123, (2)

B=Br+Bv2, (3)

Gv=C11-C12+3C445, (4)

Gr=5C11-C12C443C11-C12+4C44, (5)

GVRH=Gr+Gr2, (6)

E=9BG3B+G, (7)

v=3B-2G2(3B+G), (8)

A=2C44C11-C12. (9)

Table IV shows the set of values that were obtained using the local density approximation (LDA) for the elastic constants and elastic moduli C11,C12,C44 Young modulus E, shear modulus G, Poisson’s ratio v and Anisotropy parameter A, the B/G ratio.

It is worth noting that each elastic constant meets the conditions listed above in the case of the compounds of concern La2RuPb, Sc2RuPb. Bulk modulus also meets the requirement C12<B<C11 [54], indicating that both compounds are mechanically stable at ambient conditions with elastic deformation in the two Hg2CuTi,Cu2MnAl-type structures.

The Poisson’s ratio plays a significant role in distinguishing between covalent and Mod ionic materials, as a small value of Poisson’s ratio of around v = 0.1 characterizes covalent materials and around v = 0.25 characterizes ionic materials [55, 56]. In this study, note that the determined Poisson’s ratio values v in both types of structures of the two compounds are more significant than the critical value 0.25 of the ionic materials, indicating that these materials are ionic.

Furthermore, we refer to Frantsevich and co-authors [55], who differentiate the ductility and brittleness of materials by the B/G ratio to comprehend the ductile and brittle nature of any compound. Frantsevich [57] further links the (ductility/brittleness) behavior to the Poisson’s ratio value v, which is smaller than 0.33 for brittle materials and more than 0.33 for ductile materials.

The calculated value of the Poisson’s ratio v, for La2RuPb is less than (0.33), classifying this compound as brittle in the Hg2CuTi-type structure, whereas v for Sc2RuPb is more than 0.33, qualifying this compound as ductile in the Cu2MnAl-type structure.

The critical value for distinguishing ductile and brittle materials is 2.76 for the B/G ratio. If the B/G ratio is greater than 2.76, the material behaves in a ductile manner. Otherwise, the material behaves in a brittle manner.

In the Hg2CuTi,Cu2MnAl-type structure, the calculated values of B/G ratio are (2.4048, 2.6879) for both compounds La2RuPb, Sc2RuPb, respectively, where these values are less than the critical value of 2.76, demonstrating the brittle nature of both these materials in this type structures.

However, these compounds may be categorized as ductile materials using Pugh’s criteria (B/G) ratio [58], which distinguishes ductility from brittleness for materials with a value of 1.75.

These findings are still theoretical research, which we intend to contradict shortly with experimental experiments.

We must know the value of the shear anisotropic factor A to comprehend the elastic anisotropic degree possessed by that composite. The elastic anisotropy A has a significant meaning in industrial science as it identifies micro-fractures in materials.

Where:

  1. The composite is termed isotropic if the value of A = 1 [59, 60].

  2. The composite is deemed anisotropic if it deviates from the unit value for A (High / Low A < 1 or A > 1.

So, for La2RuPb, Sc2RuPb composites in the two phases’ Hg2CuTi,Cu2MnAl respectively, the values in Table IV above for the factor shear anisotropic A and computed by the (LDA) approximation, show anisotropic behavior and have a low probability of developing micro-cracks or structural defects during their growing process.

Table IV The calculated single crystal elastic constants C ij, and the polycrystalline elastic modulus: [shear modulus G, Voigt shear modulus G v , Reuss shear modulus G r , Young’s modulus E, the Lame’s first parameter λ in GPA, Poisson’s ratio v, B/G ratio and the shear anisotropic factor A for the La2RuPb, Sc2RuPb compounds in the Hg2CuTi, Cu2MnAl-type structures. 

Parameter La2RuPb in the Hg2CuTi-type structure Sc2RuPb in the Cu2MnAl-type structure
C 11 117.4 159.3
C 12 88.8 121.5
C 12 79.8 93.3
G 40.888 49.889
E 107.732 133.155
v 0.3174 0.3345
A 5.58 4.9365
B r 98.333 134.1
B v 98.333 134.1
G r 28.17 36.238
G v 53.6 63.54
λ 71.074 100.84
B/G 2.4048 2.6879

3.3. Thermal properties (Debye temperature, Sound velocity)

The thermal characteristics dictate the best circumstances for crystal formation to begin and sustain quality. We calculate the Debye temperature θ D, an essential basic parameter for a solid’s thermodynamic characteristics, and it is linked to numerous physical properties such as specific heat, melting temperature, and elastic stiffness constants [61].

It may also distinguish between high and low-temperature areas in solids, with a greater θ D indicating higher thermal conductivity and melting temperature. It is found in equations describing properties arising from atomic vibrations and phonon theories; the determination of θ D using elastic constants at low temperatures will be the same as that of specific heat measurements, and the vibrational excitations at low temperatures are due to acoustic modes, i.e., when θ D is connected with lattice vibrations.

Because θ D is proportional to the average sound velocity vm, one of the traditional techniques of determining the Debye temperature is by using constant elastic data [60, 61]. As a result, the following classical equation [62, 63] may be used to estimate it:

θD=hkB3n4πNAρM13vm, (10)

where h is Planck’s constant, k B is Boltzmann’s constant, N A is Avogadro’s number, ρ is the density, M is the molar mass and n is the number of atoms in the unit cell.

From specific heat measurements the θ D can be calculated by the following relation [64]:

θD=200V1/3TmM, (11)

where V the molar volume in cm3/mol, M is the molar mass in g/mol, and T m the melting point in Kelvin °K.

For the polycrystalline material, the average sound velocity vm can be calculated by the following relation [62, 65]:

vm=132vt3+1vl3-13, (12)

where vt and vl represent the transverse and longitudinal sound velocities, which are obtained by using the shear modulus G, the bulk modulus B and the density ρ by employing the Navier’s equation as found in Refs. [67, 68] and is given below:

vl=3B+4G3ρ, (13)

vt=Gρ. (14)

Table V shows the calculated results of the vt,vl,vm, and the θ D for the present compounds La2RuPb, Sc2RuPb in the Hg2CuTi, Cu2MnAl phases.

Table V The calculated longitudinal, transverse, and average sound velocity v l , v t , and v m , in ms-1 and the Debye temperature θ D, in °K for the La2RuPb, Sc2RuPb compounds in the Hg2CuTi, Cu2MnAl-type structures. 

Compound Type-structure v t v l v m (θ D )
La2RuPb Hg2CuTi 1503.126 2906.234 1682.709 209.226
Sc2RuPb Cu2MnAl 1780.955 3571.375 1998.211 269.598

There has never been a complete investigation of the characteristics of these essential chemicals in the literature before. Consequently, our current findings may anticipate future trials and other theoretical investigations, revealing its true potential for appropriate applications under various situations.

One can clearly observe from this table that the θD of the studied compounds are found to be in the following order La2RuPb < Sc2RuPb in the Hg2CuTi,Cu2MnAl-type structure respectively. This behavior is related to the atomic bonding of these materials, and this is due to the difference observed in the atomic number Z of the (Sc, La) atoms, where the (La) has a large number Z = 57 of electrons, whereas the Sc has a small number Z = 21, leading to a radius increase and thus a decrease in temperature.

3.4. Electronic and optical properties

3.4.1. Electronic properties

The study of electronic properties is crucial because it allows us to investigate and comprehend the nature of the bonds between the material’s atoms.

In addition, the state density (DOS) is a necessary physical quantity for identifying a material’s physical characteristics. It may offer access to a material’s electronic and transport characteristics and an estimate of the number of binding and anti-binder states in the population with a particular energy. It also gives us access to data on electronic conductivity. Furthermore, for each atom, the electronic density is projected onto spherical harmonics of types p, d, or f, inside a sphere of a specified radius. As a result, partial state densities are produced, allowing the chemical bonding between atoms in a crystal to be determined. The total state density projections depend on the sphere radius on which the partial state densities are projected and hence only provide qualitative information.

The electronic properties of the La2RuPb, Sc2RuPb compounds in their stable phases Hg2CuTiCu2MnAl are investigated in the following section by calculating the electronic energy band structure at zero pressure along with the principal symmetry points in the Brillouin Zone, as well as the total density of states (TDOS) and the atomic site-projected I-decomposed partial densities of states (PDOS). For this, we used the (LDA) approximation to inject the theoretical mesh parameter at equilibrium previously determined in the (LmtART) code.

Figures 2a), 2b) and 3a), 3b) show the band structure and plots of the (TDOS) and (PDOS) functions calculated using the LDA of La2RuPb, Sc2RuPb compounds in both crystallographic phases Hg2CuTi, Cu2MnAl. The Fermi level enables the determination of a material’s conduction characteristics. This research aims to figure out how the electronic population varies at Fermi’s level E F .

Figure 2 The bands structures along the symmetry lines of the Brillouin zone at the equilibrium lattice constant for the a) La2RuPb and b) Sc2RuPb compounds in the Hg2CuTi, Cu2MnAl-type structures. 

Figure 3 The total density of state (TDOS) and the partial density of state (PDOS) for the a) La2RuPb, b) Sc2RuPb compounds in the Hg2CuTi, Cu2MnAl-type structure respectively obtained by (LDA) approximation. 

Although it is a recognized fact that the self-consistent (DFT) energy-gap calculations at the level of the local density approximation (LDA) exchange-correlation functional are severely underestimated, it is clear from the band structure nature of this compound, as shown in Fig. 2a), that this compound La2RuPb is a semiconductor with the direct energy bandgap within the following value 0.6044 eV at the Γ point. Moreover, we note that our compound exhibits the Dirac cone-shaped surface energy levels responsible for the topological properties of the (TIs) materials.

It is evident that in Fig. 2b), the non-existence of a gap at the Fermi level for the Sc2RuPb compound (majority) band’s structure has metallic intersections at the Fermi level confirms the assertive metallic behavior and indicates the presence of conducting features. Whereas in Fig. 2a), the La2RuPb compound (minority) band’s structure exhibits semiconductor behavior, we note an indirect energy bandgap between KC-ΓV points, and it takes the following value 0.398 eV and a direct energy bandgap along ΓC-ΓV point 0.6044 eV. To the best of our knowledge, there are no theoretical calculations exploring the electronic band structures using the LDA approach for La2RuPb, Sc2RuPb compounds.

We note that the band profile of the calculated band structure for the La2RuPb compound is similar to that obtained by Zhang and co-authors [28], with a slight difference due to the absence of spin-orbit in our calculation, which plays a vital role in the band inversion mechanism [34, 35, 68], and in particular, the atoms in (Y-Z) zinc-blende are much heavier than those in (X - X). We should also highlight that, as Zhang et al. [28] point out, the band topology in such compounds is very sensitive to variations in lattice constant and uniaxial strain. The La2RuPb compound combination has a positive energy difference between the Γ6 and Γ8 bands ΔE=EΓ6-EΓ8, making it a trivial semiconductor that can be transformed into a nontrivial semiconductor by including spin-orbit coupling (SOC) and controlling the degree of hybridization by firm strain control (constant lattice expansion). We also highlight that our findings are quite similar to Zhang et al. [28]. The electronic states’ nature (atom and orbital) that create the energy bands may be determined from the spectra (PDOS). The (TDOS) spectrum of the La2RuPb is separated into two main areas, as shown in Fig. 3a): The first represents the valance band and indicates the mean contribution of the 4d states of the Ru atoms, as well as the low contribution of the 6p states of the Pb atoms included in the combination. The conduction band is represented by the second area, which exhibits a significant contribution from the La atoms’ 4f states, with a bit of contribution from the Ru atoms’ 4d states.

Concerning the spectrum (TDOS) of the compound Sc2RuPb shown in Fig. 3b), which is the opposite of the first compound, it can also be divided into two main sections: The first region represents the valance band, which shows the vital contribution of the 4d states of the Ru atoms in addition to the low contribution of the 3d, 6P states of the Sc and Pb atoms, respectively. The second region represents the conduction band; it shows the mean contribution of the 3d states of the (Sc) atoms with a small contribution of the 6P, 4d states of the (Pb, Ru) atoms, respectively.

3.4.2. Optical properties

The various methods of light interacting with matter in solid-state physics, including absorption, transmission, reflection, and broadcast, are of significant interest. In addition, the study of solid optical characteristics is a valuable tool in our knowledge of material electronic properties. Understanding the optoelectronic behavior of a material, in particular, requires an understanding of its optical characteristics. The frequency-dependent complex permittivity of light (dielectric function) is the most significant characteristic for describing a solid’s optical properties, and the following relationship provides it [69]:

ε(ω)=ε1(ω)+iε2(ω), (15)

where ε1(ω) and iε2(ω) are used to represent the real and imaginary parts of the dielectric function, respectively. The imaginary part is given by [70-72]:

ε2(ω)=e2Πω2m2v,c×BZMcv(k)2δωcv(k)-ωd3k. (16)

The integral is over the first Brillouin Zone. The momentum dipole elements: Mcv(k)=uck|δ.|uvk, where δ is the potential vector defining the electric field, are matrix elements for direct transitions between valence uvk(r) and conduction-band uck(r) states, and the energy ωcv(k)=Eck-Evk is the corresponding transition energy.

To deduce the real part of the complex dielectric function ε1(ω) from the imaginary part ε2(ω), we use the transformation of Kramers-Kronig [70-72]:

ε1(ω)=1+2πP0ω'ε2ω'ω'2-ω2dω', (17)

where P implies the principal value of the Cauchy integral. From the real and imaginary parts of the frequency dependent dielectric function, we can calculate the refractive index n(ω) following this relation [73-75]:

n(ω)=ε1(ω)2+ε12(ω)+ε22(ω)2 . (18)

In this study, we used the optimized structure under equilibrium condition, which is at zero pressure, to compute the refractive index n(ω) and the frequency-dependent imaginary and real parts of the dielectric function of the La2RuPb, Sc2RuPb compounds in the (Hg2CuTi2, Cu2MnAl)-type structures, respectively. Transferring unoccupied bands from occupied bands is a typical process in the electronic energy band structure, especially near the Brillouin Zone’s high symmetry points, to explain the optical spectra. The real and imaginary parts of the dielectric function that are frequency dependent are presented in Figs. 4a) and 4b). Photons with energies between [0 and 13.6] eV are used to study the optical response of these materials.

There are multiple significant peaks in the imaginary part of the dielectric function for the two compounds La2RuPb, Sc2RuPb between 1.36-3.538 eV and 0.272-2.449 eV, respectively. These peaks illustrate how the optical transitions from the valence band to the conduction band coincide.

The spectrum also indicates the development of the first optical critical point for the dielectric function at 0.6044 eV for the material La2RuPb, as shown in Fig. 4a). The fundamental absorption edge is placed at the splitting level ΓC-ΓV, representing the direct optical transition threshold between the valence band’s maximum and the conduction band’s minimum. The spectrum curve overgrows away from this point until it reaches its peak value of 51.763 for 1.36 eV, owing to an unexpected rise in the number of points contributing to the imaginary part. The spectrum of the dielectric function’s real part takes the value of 18.1676 at zero 0 eV, reaches a maximum of 38.057 at 1.08 eV, and terminates at zero at 2.1 eV.

For the metallic alloy Sc2RuPb, we have taken in consideration the strong contribution of the intraband transitions between VBM and CBM [76]. The obtained results show, for low energies, that there is a large difference between the curves of the real part and the imaginary part (with/without intraband), which confirms the strong contribution of the intraband transitions and confirms the strong metallic behavior of this alloy. For the highest energies, we see that the curves are almost similar, and it is obvious, because these transitions are mainly due to the transitions in the bands that overlap (VBM and CBM) close to the fermi level.

As a result, we can see that both compounds La2RuPb, Sc2RuPb have a particular anisotropic property.

We show the calculated dispersion of the refractive index n(ω)in the same spectral range in Fig. 5a) and b), which is one of the most critical factors for defining light propagation through the optical medium. Our computed static refractive index n(0)for the La2RuPb compound is 4.262 , whereas the maximum value of the predicted refractive index was 6.602 at 1.088 eV of photon energy. The static refractive index n(0)of the second compound Sc2RuPb was found to be 54.7671 for the curve with intraband transitions and 9.16597 for that without intraband. These latter high values are due to the metallic behavior of this alloy.

Figure 5 The real and imaginary parts of the dielectric function for the a) La2RuPb, b) Sc2RuPb compounds in the Hg2CuTi2, Cu2MnAl-type structure obtained by (LDA) approximation. 

Figure 5 The refractive index n as a function of energy for the a-La R2 R Pb, b - Sc2RuPb compounds in the Hg2CuTi2, Cu2MnAl-type structure respectively obtained by (LDA) approximation. 

4. Conclusion

The main goal of this research was to investigate the structural (phase stability), mechanical, elastic, thermal, electronic, and optical properties of Heusler compounds X2RuPb, where X = (La, Sc). The first principle (FP-LMTO) approach was used in the context of the (DFT) implemented in the (LmtART) code to carry out this research. First, the exchange-correlation potential (XC) is calculated using the LDA approximation parameterized by Perdew-Wang.

We started by looking at the structural characteristics of our X 2RuPb compounds, which showed that they are stable in the Hg2CuTi,Cu2MnAl phases. Then, at zero pressure and ambient temperature, we investigated the ground state properties (lattice parameters, bulk modulus B and its pressure derivative B’, total energy, and mechanical properties), where all elastic constants were calculated while satisfying the mechanical stability criteria for our compounds. The B/G ratios for the X2RuPb compounds indicate ductile materials. In addition, the obtained Debye temperature of La2RuPb is lower than that of Sc2RuPb

We employed the (LDA) approach to study electronic properties, including band structure and state density. With the direct energy bandgap, our compound displays topological insulator characteristics and demonstrates that the gap at E F is formed by hybridization between the d-states of the transition atoms. Furthermore, the structure of the band’s findings shows metallic behavior for Sc2RuPb and semiconducting behavior for La2RuPb, with a direct gap in the ΓC-ΓV point and an indirect gap in the KC-ΓV point. As a result, our materials seem to be semi-metallic.

In addition to their electrical properties, the optical properties of the X2RuPb compounds, including the imaginary parts of the dielectric function and the refractive index n(ω), are also examined in this study.

We discovered satisfactory results for all of the calculated physical attributes in general. The (FP-LMTO) method’s approximation (LDA) seems to be a good fit for the Heusler investigation. Finally, to the best of the authors’ knowledge, no data on these compounds can be found in the literature.

References

1. F. Heusler, Über magnetische manganlegierungen. Verhandlungen der Deutschen Physikalischen Gesellschaft 5 (1903) 219. [ Links ]

2. M. Khirat, et al., First principles calculations of the band ordering conversion in the La x Sc 2-x RuPb compounds. Solid State Communications 276 (2018) 14, https://doi.org/10.1016/j.ssc.2018.03.018 [ Links ]

3. S. Chadov, et al., Tunable multifunctional topological insulators in ternary Heusler compounds. Nature Mater 9 (2010) 541, https://doi.org/10.1038/nmat2770 [ Links ]

4. H. Lin, et al., Half-Heusler ternary compounds as new multifunctional experimental platforms for topological quantum phenomena. Nature Mater 9 (2010) 546, https://doi.org/10.1038/nmat2771 [ Links ]

5. R. A. De Groot, et al., New class of materials: half-metallic ferromagnets. Physical review letters 50 (1983) 2024, https://doi.org/10.1103/PhysRevLett.50.2024 [ Links ]

6. C. Felser, et al., Investigation of a novel material for magnetoelectronics: Co2Cr0. 6Fe0. 4Al. Journal of Physics: Condensed Matter 15 (2003) 7019, https://doi.org/10.1088/0953-8984/15/41/010 [ Links ]

7. D. Kieven, et al., I-II-V half-Heusler compounds for optoelectronics: Ab initio calculations. Phys. Rev. B 81 (2010): 075208. https://doi.org/10.1103/PhysRevB.81.075208 [ Links ]

8. L. Kautzsch, et al., Are AuPd TM (T= Sc, Y and M= Al, Ga, In), Heusler Compounds Superconductors without Inversion Symmetry?. Materials 12 (2019) 2580, https://doi.org/10.3390/ma12162580 [ Links ]

9. Ch. G. F. Blum, et al., Exploring the details of the martensite-austenite phase transition of the shape memory Heusler compound Mn 2 NiGa by hard x-ray photoelectron spectroscopy, magnetic and transport measurements. Applied Physics Letters 98 (2011) 252501, https://doi.org/10.1063/1.3600663. [ Links ]

10. B. Dieny, et al. Giant magnetoresistive in soft ferromagnetic multilayers. Phys. Rev. B. 43 (1991) 1297, https://doi.org/10.1103/PhysRevB.43.1297. [ Links ]

11. Heusler, Fr. Magnetismus und Kristallstruktur bei Manganaluminiumkupfer. Zeitschrift für anorganische und allgemeine Chemie 161 (1927) 159, https://doi.org/10.1002/zaac.19271610113. [ Links ]

12. Y. M. Cho, et al. Effects of rapid thermal annealing on the ferromagnetic properties of sputtered Zn1-x(Co 0.5 Fe 0.5) x O thin films. Applied Phys. Let. 80 (2002) 3358, https://doi.org/10.1063/1.1478146 [ Links ]

13. A. Bouazza, Sputtering of semiconductors, conductors, and dielectrics for the realization of electronics components thin-films. Int. J. Thin Film Science and Technology 11 (2022) 225, https://doi.org/10.18576/ijtfst/110210 [ Links ]

14. A. Bouazza, Deposition of Thin Films Materials used in Modern Photovoltaic Cells. Int. J. Thin. Film. Sci. Tec. 11 (2022) 313, https://doi.org/10.18576/ijtfst/110308 [ Links ]

15. R. Salah Eddine Chouaib, A. Bouazza, and Y. Belhadji, 3D sputtering simulations of the CZTS, Si and CIGS thin films using Monte-Carlo method. Monte Carlo Methods and Applications 27 (2021) 373, https://doi.org/10.1515/mcma-2021-2094 [ Links ]

16. L. Fu, L. K. Charles, and J. M. Eugene Topological insulators in three dimensions. Phys. Rev. Let. 98 (2007) 106803, https://doi.org/10.1103/PhysRevLett.98.106803 [ Links ]

17. J. E. Moore, The birth of topological insulators. Nature 464 (2010) 194, https://doi.org/10.1038/nature08916. [ Links ]

18. Bernevig, B. Andrei, Taylor L. Hughes, and Shou-Cheng Zhang. Quantum spin Hall effect and topological phase transition in HgTe quantum wells. Science 314 (2006) 1757, https://doi.org/10.1126/science.1133734. [ Links ]

19. M. Konig, et al. Quantum spin Hall insulator state in HgTe quantum wells. Science 318 (2007) 766, https://doi.org/10.1126/science.1148047. [ Links ]

20. O. Heusler, Kristallstruktur und Ferromagnetismus der Mangan-Aluminium-Kupferlegierungen. Annalen der Physik 411 (1934) 155, https://doi.org/10.1002/andp.19344110205. [ Links ]

21. T. Graf, F. Claudia and S.P.P. Stuart, Simple rules for the understanding of Heusler compounds. Progress in solid state chemistry 39 (2011) 1, https://doi.org/10.1016/j.progsolidstchem.2011.02.001. [ Links ]

22. F. Ahmadian, New Heusler Compounds Cr 2 YSb (Y= Fe, Co and Ni): Magnetic and Electronic Properties. J. superconductivity and novel magnetism 26 (2013) 1737, https://doi.org/10.1007/s10948-012-1893-6. [ Links ]

23. J. Drews, E. Ursula and S. Hans-Uwe, Optische Untersuchungen an farbigen Intermetallischen phasen. Journal of the Less Common Metals 116 (1986) 271, https://doi.org/10.1016/0022-5088(86)90235-3. [ Links ]

24. N. Tezuka, et al. Improved tunnel magnetoresistance of magnetic tunnel junctions with Heusler Co 2 FeAl 0.5 Si 0.5 electrodes fabricated by molecular beam epitaxy. Applied Phys. Lett. 94 (2009) 162504, https://doi.org/10.1063/1.3116717. [ Links ]

25. X. Dai, et al. New quarternary half metallic material CoFeMnSi. J. Applied Phy. 105 (2009) 07E901. https://doi.org/10.1063/1.3062812. [ Links ]

26. E. Ursula, W. Seelentag, and S. Hans-Uwe, Zur Kenntnis farbiger ternärer und quaternärer Zintl-Phasen/Coloured Ternary and Quaternary Zintl-Phases. Zeitschrift für Naturforschung B 35 (1980) 1341, https://doi.org/10.1515/znb-1980-1103. [ Links ]

27. W. Feng, et al. Three-Dimensional Topological Insulators in I- III- VI 2 and II- IV- V 2 Chalcopyrite Semiconductors. Phys. Rev. Lett. 106 (2011) 016402, https://doi.org/10.1103/PhysRevLett.106.016402. [ Links ]

28. X. M. Zhang, et al. Prediction of topological insulating behavior in inverse Heusler compounds from first principles. Computational materials science 70 (2013) 145, https://doi.org/10.1016/j.commatsci.2012.12.013. [ Links ]

29. A. Bouazza, Investigation using Monte-Carlo codes simulations for the impact of temperatures and high pressures on thin films quality. Rev. Mex. Fis. 69 (2023) 021501, https://doi.org/10.31349/RevMexFis.69.021501 [ Links ]

30. A. Bouazza, Simulation of the Deposition of Thin-Film Materials Used in the Manufacturing of Devices with Miniaturized Circuits. J. Sur.Investigation: X-ray, Synchrotron and Neutron Techniques 16 (2022) 1221, https://doi.org/10.1134/S1027451022060283. [ Links ]

31. A. Bouazza and S. Abderrahmane Understanding the contribution of energy and angular distribution in the morphology of thin films using Monte Carlo simulation. Monte Carlo Methods and Applications 24 (2018) 215, https://doi.org/10.1515/mcma-2018-0019. [ Links ]

32. A. Bouazza, 3D Visualization of the Effect of Plasma Temperature on Thin-Film Morphology. Bull. Lebedev Phys. Inst., 50 (2023) 7, https://doi.org/10.1134/S1027451022060283. [ Links ]

33. Y. Lakred, S. Bahlouli, and M. Elchikh. Electronic and topological properties in heavy-element based Heusler compounds X 2 RuPb (X= Sc, Y, La): first-principle method. The European Physical Journal B 94 (2021) 135, https://doi.org/10.1140/epjb/s10051-021-00141-8. [ Links ]

34. K. Mun Wong, et al. First-principles investigation of the size-dependent structural stability and electronic properties of O-vacancies at the ZnO polar and non-polar surfaces. J. Appli. Phys. 113 (2013) 014304, https://doi.org/10.1063/1.4772647 [ Links ]

35. S. Yu. Savrasov, Linear-response theory and lattice dynamics: A muffin-tin-orbital approach. Phys. Rev. B 54 (1996) 16470, https://doi.org/10.1103/PhysRevB.54.16470. [ Links ]

36. M. Larbi, et al. First Principle Investigation of Structural, Electronic and Optical Properties of Quaternary B x In y Ga 1-x-y N Compounds. Spin. 10 World Scientific Publishing Company, 2020. https://doi.org/10.1142/S2010324720500241. [ Links ]

37. W. SKohn, Density functional and density matrix method scaling linearly with the number of atoms. Phys. Rev. Lett. 76 (1996) 3168, https://doi.org/10.1103/PhysRevLett.76.3168. [ Links ]

38. J. P. Perdew, K. Burke, and M. Ernzerhof, Generalized gradient approximation made simple. Phys. Rev. Lett. 77 (1996 3865, https://doi.org/10.1103/PhysRevLett.77.3865. [ Links ]

39. S. Y. Savrasov, Program LMTART for electronic structure calculations. Zeitschrift für Kristallographie-Crystalline Materials 220 (2005) 555, https://doi.org/10.1524/zkri.220.5.555.65067. [ Links ]

40. J. P. Perdew, and W. Yue Pair-distribution function and its coupling-constant average for the spin-polarized electron gas. Phys. Rev. B 46 (1992) 12947, https://doi.org/10.1103/PhysRevB.46.12947. [ Links ]

41. P. E. Blöchl, Projector augmented-wave method. Phys. Rev. B 50 (1994) 17953. https://doi.org/10.1103/PhysRevB.50.17953. [ Links ]

42. F. D. Murnaghan, The compressibility of media under extreme pressures. Proceedings of the National Academy of Sciences 30 (1944) 244, https://doi.org/10.1073/pnas.30.9.24. [ Links ]

43. V. G. Tyuterev, and Nathalie Vast, Murnaghan’s equation of state for the electronic ground state energy. Computational materials science 38 (2006) 350, https://doi.org/10.1016/j.commatsci.2005.08.012. [ Links ]

44. K. Benyelloul, and A. Hafid Elastic constants of austenitic stainless steel: Investigation by the first-principles calculations and the artificial neural network approach. Comp. Mat. Sci. 67 (2013) 353, https://doi.org/10.1016/j.commatsci.2012.09.005. [ Links ]

45. Y. Ouyang, et al. First-principles calculations of elastic and thermo-physical properties of Al, Mg and rare earth lanthanide elements. Physica B: Condensed Matter 404 (2009) 2299, https://doi.org/10.1016/j.physb.2009.04.032. [ Links ]

46. M. J. Mehl, et al., Structural properties of ordered high-melting-temperature intermetallic alloys from first-principles total-energy calculations. Phys. Rev. B 41 (1990) 10311, https://doi.org/10.1103/PhysRevB.41.10311 [ Links ]

47. H. Rached, et al. First-principles calculations of structural, elastic and electronic properties of Ni2MnZ (Z= Al, Ga and In) Heusler alloys. Physica Status Solidi (b) 246 (2009) 1580, https://doi.org/10.1002/pssb.200844400 [ Links ]

48. F. I Fedoras, Theory of Elastic Waves in Crystals. (1985) 213. [ Links ]

49. M. Born, and K. Huang, Dynamical theory and experiment I. Publishers, Berlin (1982). [ Links ]

50. W. J. T. L. Voight, Lehrbuch der kristallphysik. Teubner, Leipzig (1928). [ Links ]

51. A. Reuß, Berechnung der fließ grenze von mischkristallen auf grund der plastizitätsbedingung für einkristalle. ZAMM-Journal of Applied Mathematics and Mechanics/Zeitschrift für Angewandte Mathematik und Mechanik 9 (1929) 49, https://doi.org/10.1002/zamm.19290090104 [ Links ]

52. R. Hill, The elastic behaviour of a crystalline aggregate. Proceedings of the Physical Society. Section A 65 (1952) 349, https://doi.org/10.1088/0370-1298/65/5/307 [ Links ]

53. B. Mayer, et al., Ab-initio calculation of the elastic constants and thermal expansion coefficients of Laves phases. Intermetallics 11 (2003) 23, https://doi.org/10.1016/S0966-9795(02)00127-9 [ Links ]

54. M. Rajagopalan, Full potential linear augmented plane wave study of the elastic properties of XPt3 (X= V, Cr, Mn, Fe, Co, Ni). Physica B: Condensed Matter 405 (2010) 2516, https://doi.org/10.1016/j.physb.2010.02.060. [ Links ]

55. J. M. Haines, and G. Bocquillon. Annu. Rev. Mater. Res 31 (2001) 1. [ Links ]

56. H. Fu, et al. Ab initio calculations of elastic constants and thermodynamic properties of NiAl under high pressures. Comp. Mat. Sci. 44 (2008) 774, https://doi.org/10.1016/j.commatsci.2008.05.026. [ Links ]

57. I. N. Frantsevich, F. F. Voronov, and S. A. Bokuta. Elastic Constants and Elastic Moduli of Metals and Insulators, ed. IN Frantsevich. (1983). [ Links ]

58. S. F. Pugh, XCII. Relations between the elastic moduli and the plastic properties of polycrystalline pure metals. The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science 45 (1954) 823. https://doi.org/10.1080/14786440808520496. [ Links ]

59. Z. E. Biskri, et al. Computational study of structural, elastic and electronic properties of lithium disilicate (Li2Si2O5) glass-ceramic. Journal of the mechanical behavior of biomedical materials 32 (2014) 345, https://doi.org/10.1016/j.jmbbm.2013.10.029. [ Links ]

60. S. Ouardi, et al. Electronic and crystallographic structure, hard x-ray photoemission, and mechanical and transport properties of the half-metallic Heusler compound Co 2 MnGe. Phys. Rev. B. 84 (2011) 155122. https://doi.org/10.1103/PhysRevB.84.155122. [ Links ]

61. Y.-Z. Zhang, et al. Flexible supercapacitors based on paper substrates: a new paradigm for low-cost energy storage. Chemical Society Reviews 44 (2015) 5181, https://doi.org/10.1039/C5CS00174A. [ Links ]

62. O. L. Anderson, A simplified method for calculating the Debye temperature from elastic constants. Journal of Physics and Chemistry of Solids 24 (1963) 909, https://doi.org/10.1016/0022-3697(63)90067-2. [ Links ]

63. I. Johnston, Solid state physics simulations. Wiley, 1996. [ Links ]

64. R. E. Newnham, Properties of materials: anisotropy, symmetry, structure. Oxford University Press on Demand, 2005. [ Links ]

65. J. Feng, et al. Stability, thermodynamic and mechanical properties of the compounds in the Ag-Sn-O system. Physica B: Condensed Matter 404 (2009) 2461, https://doi.org/10.1016/j.physb.2009.05.004. [ Links ]

66. E. Schreiber, O. L. Anderson and N. Soga, Elastic Constants, and Their Measurements, McGraw-Hill Companies, Inc., New York, 1973. [ Links ]

67. S. Daoud, et al. Effect of hydrostatic pressure on the structural, elastic and electronic properties of (B3) boron phosphide. Pramana 79 (2012) 95. https://doi.org/10.1007/s12043-012-0283-8 [ Links ]

68. H. Lin, et al. Nature Mater. 9 (2010) 546. [ Links ]

69. H. Bouafia, et al. Insight into elastic anisotropy, mechanical and dynamical stability, electronic properties, bonding and weak interactions analysis of LuAuSn Half-Heusler. Solid State Sciences 118 (2021) 106677, https://doi.org/10.1016/j.solidstatesciences.2021.106677. [ Links ]

70. C. Ambrosch-Draxl, and J. O. Sofo, Linear optical properties of solids within the full-potential linearized augmented planewave method. Computer physics communications 175 (2006) 1, https://doi.org/10.1016/j.cpc.2006.03.005. [ Links ]

71. M. Bendjemai, et al. Insight into the role of weak interactions on optoelectronic properties of LiGaTe2-chalcopyrite under pressure effect: DFT-D3, NCI and QTAIM investigations. Physica B: Condensed Matter 599 (2020) 412463, https://doi.org/10.1016/j.physb.2020.412463. [ Links ]

72. B. Sana, et al. Study of magnetic and optoelectronic properties of BaCmO3-cubic perovskite after the estimation of Hubbard interaction and Hund’s exchange parameters: GW and DFT+ U investigations. Optik 168 (2018) 196, https://doi.org/10.1016/j.ijleo.2018.04.064 [ Links ]

73. X-G., Wang, et al. The hematite (α - Fe203)(0001) surface: Evidence for domains of distinct chemistry. Phys. Rev. Let. 81 (1998) 1038. https://doi.org/10.1103/PhysRevLett.81.1038 [ Links ]

74. A. Delin, et al. Optical properties of the group-IVB refractory metal compounds. Phys. Rev. B. 54 (1996) 1673. https://doi.org/10.1103/PhysRevB.54.1673. [ Links ]

75. Y. U. Peter, and M. Cardona, Fundamentals of semiconductors: physics and materials properties. Springer Science and Business Media, 2010. [ Links ]

76. B. S. Mendoza, and W. Luis Mochán. Ab initio theory of the Drude plasma frequency. JOSA B 38 (2021) 1918, https://doi.org/10.1364/JOSAB.416741 [ Links ]

Data availability statement. All relevant data are included in the paper All the data used in the manuscript are produced from the numerical simulations by the authors.

Received: September 08, 2022; Accepted: March 25, 2023

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License