First-Principles Studies on the Physical Properties of the Half Heusler RbNbCd and RbNbZn Compounds: A Promising Material for Thermoelectric Applications

: This work focuses on study of the structural, electronic, thermodynamic and thermoelectric properties of RbNbCd and RbNbZn Half Heusler (HH), utilizing a full-potential linearized augmented plane wave (FP-LAPW) approach and the Boltzmann transport equation using a constant relaxation time approximation within the context of density functional theory (DFT) as embedded in the WIEN2k code. The structural analysis employed the generalized gradient approximation (GGA) and considered the Birch Murnaghan equation of state (EOS), which results in the stable phase for RbNbCd and RbNbZn. The positive phonon spectra indicate the dynamical stability of the studied RbNbCd and RbNbZn. The compounds under investigation that have no bandgap are metallic, as evidenced by their electronic properties. Their mechanical and thermal stability as well as their anisotropic and ductile character are conﬁrmed by the various elastic and thermodynamic parameters. The lattice thermal conductivity has been calculated. This thorough analysis demonstrates the applicability of the studied RbNbCd and RbNbZn for thermoelectric applications.


Introduction
Researchers have recently concentrated on the investigation of smart and multifunctional materials. Applications for these materials include sensors, thermoelectricity, optoelectronics, and various forms of renewable energy [1][2][3]. It is quite concerning that the global ecology is in an alarming state due to continuously rising carbon emission rates. The need to find effective, clean and environmentally friendly renewable energy sources has therefore intensified in order to stop the destruction of our environment. Biodiesel, biomass, wind power, solar power and thermoelectric power are considered potential resources in this regard [4][5][6]. Therefore, we require such materials that can function as thermoelectric and photovoltaic generators. In contrast to thermoelectric generators, which utilize heat that has been dissipated to produce electricity, photovoltaics can absorb and transform the tremendous quantity of solar energy that the sun emits continually into useful electrical energy [5,[7][8][9].
Numerous substances have been investigated as prospective components for solar and thermoelectric systems, including chalcogenides, triple point metals, oxide perovskites, skutterudites, and Half Heusler compounds [10][11][12][13][14]. Some of the peculiar properties of these mentioned materials include a simple crystal structure, magnetism, half metallicity, a topological insulator, optoelectronics, and thermoelectricity, etc. [15][16][17]. In 1903, Friedrich Heusler made the discovery of the Cu 2 MnAl prototype Heusler compound [18]. The Heusler compound received significant attention only after its first principles calculations revealed its unique property, the half metal ferromagnet (HMF), which revolutionized the developing field of thermoelectrics and spintronics. The many forms of Heusler compounds have been the subject of several studies, and they have a variety of properties that can be exploited in a number of new technological devices. A half-Heusler compound, which has a total of 18 valence electron counts, is a novel prospective thermoelectric material, with higher electrical conductivity than other members of the Heusler compound family.
With the search of novel compounds for the thermoelectric applications, CrVZ (Z = S, Se, and Te) were investigated by Javed et al., and they reported their mechanical as well as dynamic stability [26]. Furthermore, NbCoSn and NbFeSb HHs that were reported by Zerrouki et al. were found to be mechanically stable, elastically anisotropic, and brittle by nature [27]. The potential use of these materials in solar and thermoelectric applications was also addressed. The half metallicity of RbSrX (X = C, Si, and Ge) was reported by Ahmad et al. [28]. The computed magnetic moment of these compounds was found to be 1 µB. A literature survey indicated that the structural, electronic, magnetic and thermoelectric properties of RbNbCd and RbNbZn materials have not received significant attention. This motivates us to analyze the structural, electronic, elastic, thermodynamic and thermoelectric properties of RbNbCd and RbNbZn. We believe that the findings of this work will be beneficial in future theoretical and experimental studies of additional HHs.
The first section of the manuscript is an introduction, while the second section introduces the computing strategy. The results are then presented together with an analysis and discussion of them in the third section. The conclusion and the key successes are included in the final section.

Computational Method
To analyze the structural as well as the dynamic stabilities of the RbNbCd and RbNbZn including the electronic, magnetic and thermoelectric properties, we employed the fullpotential linearized augmented plane wave approach (FP-LAPW) [29], and relied on the density functional theory (DFT), as embedded in Wien2K [30]. The FP-LAPW technique was used to expand the Kohn Sham orbitals in the muffin-tin (MT) atomic spheres into atomiclike orbitals and plane waves in the interstitial space [29,31]. The analysis of the physical characteristics of the RbNbCd and RbNbZn Perdew-Burke Ernzerhof generalized gradient approximation (PBE-GGA) [32] was taken into consideration for structure optimization.
For computation of the electronic properties, the modified Becke-Johnson potential (mBJ) [33,34] was utilized, as it precisely predicts the ground state electronic properties. Inside the atomic spheres, the valence wave functions were enlarged to l = 10 partial waves. In this computation, the plane wave expansion R MT K max = 7 was used, where K max is the magnitude of the largest K vector and R MT is the size of the smallest muffin tin sphere. A muffin-tin (MT) model was used to predict the crystal potential, and the assumed R MT values for Rb, Nb, Cd and Zn were set at 2.5, 1.80, and 1.63, respectively [35]. The potential and charge density were extended up to G max = 12 a.u −1 . Using the modified tetrahedron approach, a dense mesh of 10 × 10 × 10 uniformly distributed k-points (equivalent to 405 in irreducible BZ) was used inside the Brillouin zone (BZ), where the total energy concentration was set up to 10 −6 Ry. The elastic constants were calculated using the IRELAST method [36], which was employed in the WIEN2k. Furthermore, the Gibbs2 code was utilized to compute the thermodynamic properties of RbNbCd and RbNbZn [37]. The thermoelectric properties of this molecule were examined using the semi-classical Boltzmann theory and the BoltzTrap method [38]. Moreover, by employing the Vienna Ab initio Simulation Package (VASP) [39], we performed the phonon dispersion study in order to study the dynamic stability of the investigated RbNbCd and RbNbZn.

Structural Properties
The analysis of the structural characteristics of RbNbCd and RbNbZn is the focus of this section. This type of research holds great promise because it enables the gathering of data related to the microscopic structure that assists in the forecasting of other properties, such as the electronic, elastic, and thermodynamic. The HHs generally belong to the generic formula XYZ that exists on the non-centrosymmetric cubic structure (#216, F-43m) [40]. However, the HH structure is shown in Figure 1 by the three different cubic face-centered networks, or X, Y, and Z, that are coupled with one another. This kind of structure can exist in any one of the three non-equivalent atomic configurations listed in Table 1. Structural optimization has been carried out for the ferromagnetic (FM) and non-magnetic (NM) phases that were carried out on the three diverse configurations (α, β, γ), as represented in Figure 1. (equivalent to 405 in irreducible BZ) was used inside the Brillouin zone (BZ), where the total energy concentration was set up to 10 −6 Ry. The elastic constants were calculated using the IRELAST method [36], which was employed in the WIEN2k. Furthermore, the Gibbs2 code was utilized to compute the thermodynamic properties of RbNbCd and RbNbZn [37]. The thermoelectric properties of this molecule were examined using the semi-classical Boltzmann theory and the BoltzTrap method [38]. Moreover, by employing the Vienna Ab initio Simulation Package (VASP) [39], we performed the phonon dispersion study in order to study the dynamic stability of the investigated RbNbCd and RbNbZn.

Structural Properties
The analysis of the structural characteristics of RbNbCd and RbNbZn is the focus of this section. This type of research holds great promise because it enables the gathering of data related to the microscopic structure that assists in the forecasting of other properties, such as the electronic, elastic, and thermodynamic. The HHs generally belong to the generic formula XYZ that exists on the non-centrosymmetric cubic structure (#216, F-43m) [40]. However, the HH structure is shown in Figure 1 by the three different cubic facecentered networks, or X, Y, and Z, that are coupled with one another. This kind of structure can exist in any one of the three non-equivalent atomic configurations listed in Table  1. Structural optimization has been carried out for the ferromagnetic (FM) and non-magnetic (NM) phases that were carried out on the three diverse configurations (α, β, γ), as represented in Figure 1.    To determine the equilibrium lattice constant as well as the structural stability of RbNbCd and RbNbZn, Birch Murnaghan EOS [41] was considered, with the variation of total energy as a function of volume, as displayed in Figure 2. Because the energy-volume curve of γ in the FM arrangement is smaller than that of the other two arrangements, this suggests its stability. Table 1 lists the equilibrium energy, bulk modulus, as well as its pressure derivative for RbNbCd and RbNbZn, along with the optimized lattice parameters. The formation energy was used to evaluate the thermodynamic stability that affects the possibility of synthesizing RbNbCd and RbNbZn. The computed formation energies for the RbNbCd and RbNbZn are represented in Table 1. The negative formation energy indicates their stability [42]. To confirm the stability of the RbNbCd and RbNbZn, we estimated the phonon dispersion (PD) against momentum, as shown in Figure 3. The transverse and longitudinal modes in the phonon dispersion diagram typically resemble spaghetti. With a range of 0 to 6 THz, we employed frequency (THz) along the y-axis and momentum along the x-axis. However, practically all of the modes in the RbNbCd and RbNbZn compounds are positive and have genuine phonon branches, reflecting their stability [43,44]. To determine the equilibrium lattice constant as well as the structural stability of RbNbCd and RbNbZn, Birch Murnaghan EOS [41] was considered, with the variation of total energy as a function of volume, as displayed in Figure 2. Because the energy-volume curve of γ in the FM arrangement is smaller than that of the other two arrangements, this suggests its stability. Table 1 lists the equilibrium energy, bulk modulus, as well as its pressure derivative for RbNbCd and RbNbZn, along with the optimized lattice parameters. The formation energy was used to evaluate the thermodynamic stability that affects the possibility of synthesizing RbNbCd and RbNbZn. The computed formation energies for the RbNbCd and RbNbZn are represented in Table 1. The negative formation energy indicates their stability [42]. To confirm the stability of the RbNbCd and RbNbZn, we estimated the phonon dispersion (PD) against momentum, as shown in Figure 3. The transverse and longitudinal modes in the phonon dispersion diagram typically resemble spaghetti. With a range of 0 to 6 THz, we employed frequency (THz) along the y-axis and momentum along the x-axis. However, practically all of the modes in the RbNbCd and RbNbZn compounds are positive and have genuine phonon branches, reflecting their stability [43,44].

Elastic Properties
The elastic properties required for the study of diverse mechanical properties such as strength of material, ductile and brittle nature, influence on stress and strain, as well as interatomic bonding can be determined from the first principle method based on the density functional theory [45,46]. The mechanical behavior of the investigated substance is quantified from the elastic properties. These characteristics have a significant impact on the compounds selected for innovative photovoltaic applications [47]. Three distinct elastic constants exist for cubic compounds (Cij). Table 2 provides an overview of the computed elastic constants for RbNbCd and RbNbZn. The Born-Huang stability condition [48,49] is used to determine the mechanical stability of both of these compounds, which is expressed as C 11 − C 12 > 0, C 11 > 0, C 44 > 0 and (C 11 + 2C 12 ) > 0 (1)

Elastic Properties
The elastic properties required for the study of diverse mechanical properties such as strength of material, ductile and brittle nature, influence on stress and strain, as well as interatomic bonding can be determined from the first principle method based on the density functional theory [45,46]. The mechanical behavior of the investigated substance is quantified from the elastic properties. These characteristics have a significant impact on the compounds selected for innovative photovoltaic applications [47]. Three distinct elastic constants exist for cubic compounds (Cij). Table 2 provides an overview of the computed elastic constants for RbNbCd and RbNbZn. The Born-Huang stability condition [48,49] is used to determine the mechanical stability of both of these compounds, which is expressed as C11 − C12 > 0, C11 > 0, C44 > 0 and (C11 + 2C12) > 0 (1) Satisfying the above condition, the mechanical stability of the RbNbCd and RbNbZnis was established [50,51]. Using Voigt-Reuss-Hill (VRH) approximations [52], these compounds' physical properties, including their bulk modulus (B) and shear modulus (G), were determined. The bulk modulus values derived by applying VRH approximations (Table 2) agree well with the EOS Bulk modulus (Table 1), illustrating the high internal consistency of the research. Other elastic constants such as the Young's modulus (Y), Poisson's ratio (v), and shear anisotropy factor (A) of these compounds are computed as follows:  Satisfying the above condition, the mechanical stability of the RbNbCd and RbNbZnis was established [50,51]. Using Voigt-Reuss-Hill (VRH) approximations [52], these compounds' physical properties, including their bulk modulus (B) and shear modulus (G), were determined. The bulk modulus values derived by applying VRH approximations (Table 2) agree well with the EOS Bulk modulus (Table 1), illustrating the high internal consistency of the research. Other elastic constants such as the Young's modulus (Y), Poisson's ratio (v), and shear anisotropy factor (A) of these compounds are computed as follows: When a material is compressed, its bulk modulus (B) reveals how resistant it is to volume change. The high B value denotes good crystal strength, whereas the shear modulus (G) shows the crystal opposition offered for plastic deformation [53][54][55]. For RbNbCd and RbNbZn compounds, the estimated magnitudes of B are 18.11 GPa and 16.66 GPa, respectively. RbNbZn is more robust to volume change than RbNbCd because it has a higher B value.
Additionally, RbNbZn has a higher G value (10.07) than RbNbCd (9.59). As a result, RbNbZn has higher resistance to the transverse bending. For the determination of brittleness or ductility, it is quantified by Pugh's ratio (B/G) [56,57]. Having more than 1.75 B/G, the substance is ductile or brittle. Table 2 displays the estimated B/G for RbNbCd and RbNbZn compounds, demonstrating the ductility of the compounds in question. Frantsevich et al. [58] used Poisson's ratio (v) to distinguish between brittleness and ductility, with a threshold value 0.26. If v is greater than 0.26 it is considered to be ductile or brittle. The RbNbCd and RbNbZn compounds have estimated values of 0.27 and 0.24, respectively, showing that RbNbCd is ductile by nature, whereas RbNbZn is brittle. This is further supported by Cauchy pressure C p = C 12 − C 44 , as represented in Table 2 [59,60]. A positive (negative) magnitude identifies ionic (covalent) bonding. The positive C p for RbNbCd signifies ionic bonding, whereas a negative C p for RbNbZn signifies covalent bonding. The capacity for creating microcracks in materials has been strongly associated with elastic anisotropy. The anisotropy factor (A) has a value of unity for an entirely isotropic system; the divergence from unity represents the degree of elastic anisotropy [61][62][63]. RbNbZn is anisotropic and more prone to developing microcracks or structural faults during the growth process than RbNbCd according to the computed values of A, which are displayed in Table 2.
Two thermodynamic factors associated with elastic properties: the melting point T m and the Debye temperature (θ D ), were also studied: where n is the number of atoms in a single cell, k is Boltzmann's constant, h is Planck's constant, N A is Avogadro's number, ρ is the material's density, and M is its molecular weight. Table 3 contains the computed Debye temperature. The investigated materials' mean sound speeds were determined using, where the shear modulus and bulk modulus, respectively, are used to determine the longitudinal and transverse components of sound velocity as,  Table 4 presents the computed magnitude of υ t , υ l , V m , θ D , and T m .

Magnetic Properties
When examining the material magnetic properties, the magnetic moment is a crucial consideration. It provides information on the material magnetic field components as well as the magnetic field frequency [23]. The magnetic moments of spins in RbNbCd and RbNbZn compounds were investigated. Table 3 depicts the calculated total and partial magnetic moments with polarized spin in the interstitial sites and the muffin-tin spheres. The total magnetic moment (µ B ) per cell unit for both materials approached an integer Bohr magneton of 4 µ B , confirming their half-metallic character. It is clear that although the magnetic moments of the Rb and Cd/Zn atoms are small and do not significantly contribute to the total moment, they are extremely important for the development of magnetism. Furthermore, the environment of these atoms has an impact on the width of the predominant band gap in the compound RbNbCd and RbNbZn. A significant magnetic moment can be discovered in the interstitial region, as demonstrated in Table 3. This outcome is the result of a large charge transfer from Rb and Nb to Cd/Zn, which can be used to explain the magnetic moment at the interstitial region. This is explained by the fact that the Cd/Zn anion has a higher electronegativity than the Rb and Nb, and is responsible for the magnetic moment in the interstitial region.

Electronic Properties
The dispersion of electrons in the valence and conduction bands and the width of the gap are the primary determinants of the electronic characteristics, band structure, and state density [64,65]. These characteristics are determined for Half Heusler RbNbCd and RbNbZnin equilibrium states using the most stable γ-structure optimized lattice parameter. Figure 4 shows the computed spin-up as well as the spin-down polarized electronic band structures for RbNbCd and RbNbZn. For both RbNbCd and RbNbZn compounds, the results of the band structure investigation demonstrate that both compounds are metallic and that there is no energy difference between the conduction and valence bands. Similar behavior is obtained for the spin down state.
We computed the density of states (Dos, both total and partial) for RbNbCd and RbN-bZn, as, depicted in Figure 5, to better understand the interaction between different orbitals of the atoms. With the exception of the values of their densities of states with changing energies, the computed band structure and DoS for the compounds under consideration behave similarly. In all the materials taken into consideration, the electrons of the d-orbital of the 'Nb' atom make up the majority of the contribution in the lower valence band region from −2.5 eV to 0 eV, while the electrons of the p-orbital of the Y-atom (Y = Zn, Cd) make up the minor portion at the conduction band with small contribution of the Rb atom. The metallic behavior is again confirmed from the DoS plot, with the valance state crossing the Fermi energy [64].

Thermodynamic Properties
Using Gibbs2 [37], we estimated the thermodynamic properties of the investigated RbNbCd and RbNbZn materials. Within the temperature range of 0-800 K, the thermal characteristics of the compounds under consideration were calculated. The effects of pressure in the range of 0-15 GPa were taken into consideration by using the quasi-harmonic model. The essential physical characteristic that indicates whether a material is suitable for usage in high temperature situations is the lattice or phonon thermal conductivity (k ph ) [66]. Figure 6 shows the k ph values determined at room temperature (300 K. Its dependency on temperature between 0 to 1000 K is represented, which demonstrates with increase in temperature the k ph value decreases. At 300 K the k ph is found to be 1.1 W/mk and 0.98 W/mk for RbNbCd and RbNbZn, respectively. These values are significantly lower We computed the density of states (Dos, both total and partial) for RbNbCd and RbNbZn, as, depicted in Figure 5, to better understand the interaction between different orbitals of the atoms. With the exception of the values of their densities of states with changing energies, the computed band structure and DoS for the compounds under consideration behave similarly. In all the materials taken into consideration, the electrons of the d-orbital of the 'Nb' atom make up the majority of the contribution in the lower valence band region from −2.5 eV to 0 eV, while the electrons of the p-orbital of the Y-atom (Y = Zn,Cd) make up the minor portion at the conduction band with small contribution of the Rb atom. The metallic behavior is again confirmed from the DoS plot, with the valance state crossing the Fermi energy [64].   We computed the density of states (Dos, both total and partial) for RbNbZn, as, depicted in Figure 5, to better understand the interaction bet orbitals of the atoms. With the exception of the values of their densities changing energies, the computed band structure and DoS for the compoun sideration behave similarly. In all the materials taken into consideration, t the d-orbital of the 'Nb' atom make up the majority of the contribution in lence band region from −2.5 eV to 0 eV, while the electrons of the p-orbital (Y = Zn,Cd) make up the minor portion at the conduction band with small c the Rb atom. The metallic behavior is again confirmed from the DoS plot, w state crossing the Fermi energy [64].  for usage in high temperature situations is the lattice or phonon thermal conductivity (kph) [66]. Figure 6 shows the kph values determined at room temperature (300 K. Its dependency on temperature between 0 to 1000 K is represented, which demonstrates with increase in temperature the kph value decreases. At 300 K the kph is found to be 1.1 W/mk and 0.98 W/mk for RbNbCd and RbNbZn, respectively. These values are significantly lower than the lattice thermal conductivity of other reported thermoelectric materials; however, experimental work on this will be more beneficial. Variations in the volume of the unit cell versus temperature at a diverse pressure (0, 5, 10, 15 GPa) are shown in Figure 7. It is evident that when the pressure increases and temperature remains constant, the volume decreases. Contrarily, the volume increases with the temperature at a very slow rate at a given pressure for RbNbCd and RbNbZn. Examining the data in Figure 7 reveals that the volume absolute zero pressure and zero temperature is 680 bohr 3 and 775 bohr 3 for RbNbCd and RbNbZn, respectively.
The variability in heat capacity (CV) for various temperatures for the pressure range of 0 to 15 GPa is illustrated in Figure 8. We can infer from the figure that Cv complies with the T 3 law. It is simpler to understand the vibrational properties of the material by taking into account heat capacity, which has multiple uses. At high temperatures, the value of Cv is nearly constant due to the Dulong-Petit limit [67,68], which is the result of the suppression of an anharmonic. The variation of Debye temperature (θD) [69] with diverse temperature and pressure are represented in Figure 9. With the rise in temperature at a given pressure, the pressure reduces almost linearly; however, it increases for increasing pressure at a certain temperature. Therefore, the Debye temperature is significantly affected by pressure as well as temperature in the opposite directions. Variations in the volume of the unit cell versus temperature at a diverse pressure (0, 5, 10, 15 GPa) are shown in Figure 7. It is evident that when the pressure increases and temperature remains constant, the volume decreases. Contrarily, the volume increases with the temperature at a very slow rate at a given pressure for RbNbCd and RbNbZn. Examining the data in Figure 7 reveals that the volume absolute zero pressure and zero temperature is 680 bohr 3 and 775 bohr 3 for RbNbCd and RbNbZn, respectively.  The variability in heat capacity (C V ) for various temperatures for the pressure range of 0 to 15 GPa is illustrated in Figure 8. We can infer from the figure that C v complies with the T 3 law. It is simpler to understand the vibrational properties of the material by taking into account heat capacity, which has multiple uses. At high temperatures, the value of C v is nearly constant due to the Dulong-Petit limit [67,68], which is the result of the suppression of an anharmonic. The variation of Debye temperature (θ D ) [69] with diverse temperature and pressure are represented in Figure 9. With the rise in temperature at a given pressure, the pressure reduces almost linearly; however, it increases for increasing pressure at a certain temperature. Therefore, the Debye temperature is significantly affected by pressure as well as temperature in the opposite directions.   These findings suggest that temperature and pressure have an impact on particle vibration frequency across all computed compounds. The fluctuation of the coefficient of thermal expansion (α) with temperature and pressure is shown in Figure 10. It can be seen in the figure that α rises sharply between 0 and 200 K, then rises gradually, and above 300 K, particularly in high-pressure locations, it rises almost linearly, with a low slope. The thorough examination of the thermodynamic characteristics will shed light on the experimental efforts.

Thermoelectric Properties
The capacity of thermoelectric materials to transform thermal energy into usable energy and vice versa has increased interest in them significantly around the globe [70]. Because of their capacity for energy production, dielectric properties, and refrigeration, these the efficiency of these materials finds a significant demand in the thermoelectric realm. The BoltzTraP code [38] was employed to compute thermoelectric parameters using an approximation of charge carrier constant relaxation times. This section focuses on analyzing several transport coefficients of RbNbCd and RbNbZn in both up-and down-spin states, including electrical conductivity (σ/τ), Seebeck coefficient (S), and the merit factor (ZT). For RbNbCd and RbNbZn, the variation in Seebeck coefficient (S) is represented in Figure 11. Over the whole temperature range, it can be seen that the Seebeck coefficient for RbNbCd is positive for the up-spin state, whereas S for RbNbZn is negative up to 250 K, which ends up being positive after a further increment in the temperature. In the spin-down condition, we also see that the Seebeck coefficient for the RbNbCd and RbNbZn increases linearly with temperature in the entire temperature range. These findings suggest that temperature and pressure have an impact on particle vibration frequency across all computed compounds. The fluctuation of the coefficient of thermal expansion (α) with temperature and pressure is shown in Figure 10. It can be seen in the figure that α rises sharply between 0 and 200 K, then rises gradually, and above 300 K, particularly in high-pressure locations, it rises almost linearly, with a low slope. The thorough examination of the thermodynamic characteristics will shed light on the experimental efforts. These findings suggest that temperature and pressure have an impact on particle vibration frequency across all computed compounds. The fluctuation of the coefficient of thermal expansion (α) with temperature and pressure is shown in Figure 10. It can be seen in the figure that α rises sharply between 0 and 200 K, then rises gradually, and above 300 K, particularly in high-pressure locations, it rises almost linearly, with a low slope. The thorough examination of the thermodynamic characteristics will shed light on the experimental efforts.  ing several transport coefficients of RbNbCd and RbNbZn in both up-and down-spin states, including electrical conductivity (σ/τ), Seebeck coefficient (S), and the merit factor (ZT). For RbNbCd and RbNbZn, the variation in Seebeck coefficient (S) is represented in Figure 11. Over the whole temperature range, it can be seen that the Seebeck coefficient for RbNbCd is positive for the up-spin state, whereas S for RbNbZn is negative up to 250 K, which ends up being positive after a further increment in the temperature. In the spindown condition, we also see that the Seebeck coefficient for the RbNbCd and RbNbZn increases linearly with temperature in the entire temperature range. Furthermore, Figure 12 depicts the fluctuation of electrical conductivity (σ/τ) with temperature for RbNbCd and RbNbZn. For both Half Heuslers, the electrical conductivity computed for the spin-up state and that discovered for the spin-down state are different. When the temperature is between 100 K-800 K, the σ/τ is observed to increase with temperature for the RbNbY for spin-up. However, for spin-down, the σ/τ unexpectedly dropped between 200-800 K. Finally, to measure the thermoelectric efficiency of these compounds, we examined the fluctuation of the figure of merit throughout a temperature range from 0 to 800 K, as represented in Figure 13. For RbNbCd, the magnitude of ZT rises with temperature for the spin-up state, whereas an opposite trend is observed for the RbNbZn. For RbNbCd and RbNbZn, the ZT values at 300 K for the spin-up state are around 0.007 and 0.035, respectively. The magnitude of ZT rises as the temperature rises for the spin-down state. The computed results indicate that RbNbCd and RbNbZn are suitable for thermoelectric applications at low temperatures. Furthermore, Figure 12 depicts the fluctuation of electrical conductivity (σ/τ) with temperature for RbNbCd and RbNbZn. For both Half Heuslers, the electrical conductivity computed for the spin-up state and that discovered for the spin-down state are different. When the temperature is between 100 K-800 K, the σ/τ is observed to increase with temperature for the RbNbY for spin-up. However, for spin-down, the σ/τ unexpectedly dropped between 200-800 K. Finally, to measure the thermoelectric efficiency of these compounds, we examined the fluctuation of the figure of merit throughout a temperature range from 0 to 800 K, as represented in Figure 13. For RbNbCd, the magnitude of ZT rises with temperature for the spin-up state, whereas an opposite trend is observed for the RbNbZn. For RbNbCd and RbNbZn, the ZT values at 300 K for the spin-up state are around 0.007 and 0.035, respectively. The magnitude of ZT rises as the temperature rises for the spin-down state. The computed results indicate that RbNbCd and RbNbZn are suitable for thermoelectric applications at low temperatures.

Conclusions
We analyzed the structural, electronic, elastic, thermodynamic and thermoelectric features of RbNbCd and RbNbZn employing first-principle computations and Boltzmann transport theory. The elastic properties demonstrate that RbNbCd and RbNbZn are mechanically stable. Furthermore, the computed lattice thermal conductivity values for RbNbCd and RbNbZn at 300 K are 1.1 (0.98) Wm −1 K −1 . The best thermoelectric properties are believed to be produced by the lattice thermal conductivity with the least magnitude. We also discovered that the Seebeck coefficient and electrical conductivity yield good results to match the performance requirements of thermoelectric devices with low thermal conductivity. According to the calculated results, RbNbCd and RbNbZn are suitable for thermoelectric applications at low temperatures. The aforementioned study will motivate

Conclusions
We analyzed the structural, electronic, elastic, thermodynamic and thermoelectric features of RbNbCd and RbNbZn employing first-principle computations and Boltzmann transport theory. The elastic properties demonstrate that RbNbCd and RbNbZn are mechanically stable. Furthermore, the computed lattice thermal conductivity values for RbNbCd and RbNbZn at 300 K are 1.1 (0.98) Wm −1 K −1 . The best thermoelectric properties are believed to be produced by the lattice thermal conductivity with the least magnitude. We also discovered that the Seebeck coefficient and electrical conductivity yield good results to match the performance requirements of thermoelectric devices with low thermal conductivity. According to the calculated results, RbNbCd and RbNbZn are suitable for thermoelectric applications at low temperatures. The aforementioned study will motivate researchers to pursue this avenue; however, rubidium is an element that is extremely sensitive to air and moisture, hence any experimental operation has to be performed with special precautions in place.  Data Availability Statement: The raw/processed data can be made available with reasonable request to corresponding author. Acknowledgments: One of the authors (Debidatta Behera) is thankful to the Birla Institute of Technology, Ranchi for the award of an institute Fellowship.

Conflicts of Interest:
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.