High structural stability, reduced lattice-thermal conductivity, and elevated energy harvesting efficiency in a Lu2CoCrO6
Samia Shahzadi, A. Elfasakhany, S. Nazir
TL;DR
Lu2CoCrO6 is a stable double perovskite oxide with low thermal conductivity and high thermoelectric efficiency, making it promising for energy conversion devices.
Contribution
The study reveals Lu2CoCrO6's structural stability, magnetic properties, and high thermoelectric performance through theoretical analysis.
Findings
Lu2CoCrO6 is thermodynamically, dynamically, and mechanically stable with a negative formation enthalpy and no phonon instabilities.
The material exhibits a semiconducting state with a direct band gap of 1.13 eV and a ferrimagnetic phase due to antiferromagnetic coupling.
Thermoelectric analysis shows a giant figure-of-merit of 1.00 at 700 K due to reduced lattice thermal conductivity.
Abstract
Double perovskite oxides have been emerged as promising candidates for the fast evolving technical frontier, playing a key role in the development of efficient energy conversion devices to address global energy challenges. Therefore, we theoretically examined the structural stabilities, thermoelectric, electronic, and magnetic aspects of the ordered Lu2CoCrO6 structure. The calculated negative formation enthalpy (−4.2 eV per atom), lack of imaginary modes in the phonon curves, and elastic constants that meet the Born conditions, confirms the thermodynamical, dynamical, and mechanical stability of the system, respectively. The material is classified as ductile by Pugh's ratio and Poisson's ratio (ν = 0.308 > 0.26). ratio (0.438), along with a ν value, which affirms ionic bonding. Also, the material exhibits a semiconducting state having direct a band-gap of 1.13 eV. The…
Genes, proteins, chemicals, diseases, species, mutations and cell lines named across the full text — each resolved to its canonical identifier and authoritative record.
Click any figure to enlarge with its caption.
Fig. 1
Fig. 2
Fig. 3
Fig. 4
Fig. 5
Fig. 6
Fig. 7
Fig. 8
Fig. 9
Fig. 10
Fig. 11|
|
|
|
|
|
|
|
|
|
|
|
|
|---|---|---|---|---|---|---|---|---|---|---|---|
| 244.12 | 353.92 | 1.45 | 0.90 | 1.53 | 1.69 | 95.28 | 146.17 | 1.53 | 0.186 | 0.426 | 2.28 |
- —Taif University10.13039/501100006261
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Taxonomy
TopicsMultiferroics and related materials · Heusler alloys: electronic and magnetic properties · Advanced Thermoelectric Materials and Devices
Introduction
1
The global shift towards advanced technology and sustainable energy has led to a growing need for materials that can meet the high-performance requirements of modern applications.^1–3^ Materials that exhibit strong electronic and magnetic aspects, while maintaining long-term stability under varying temperatures (temp.) and environmental conditions, are essential for applications in energy storage, spintronics, and environmental remediation.^4^ The traditional materials frequently fail to meet these requirements, which leads to the investigation of novel chemical groups with improved features. The transition metals (TM) and rare-earth-based double perovskite oxides (DPOs), having the general formula A_2_BB′O_6_, display significant structural diversity, making them hot candidates for a range of applications such as optical, electrical, magnetic, and catalytic uses.^5^ Due to their strong spin polarization, ferromagnetic (FM) phase stability, and excellent Curie temp. (TC), DPO has a number of appealing possibilities.^6^ Furthermore, these materials show remarkable features within their family, including insulating ferrimagnetic (FIM) states,^7^ magneto-dielectric features, magneto-optic aspects,^8^ half-metallic (HM) behavior,^9^ and multi-ferroic^10^ capabilities. Along with this, a number of DPOs have recently been discovered keeping interesting applications in the fields of photovoltaic devices,^11^ photocatalysis,^12^ and photoelectro chemical energy storage systems.^13^ Also due to the intriguing properties, including magnetodielectric effect,^14,15^ colossal magnetoresistance effect,^8^ multiferriocity behavior,^16^ La_2_CoMnO_6_ has emerged as the subject of much matter. Likewise, La_2_CoMnO_6_ and La_2_NiMnO_6_ represents the most extensively examined ordered DPOs, exhibiting FM semiconducting (SC) phase with TC near room temp. (220 K for La_2_CoMnO_6_ and 280 K for La_2_NiMnO_6_).^17,18^ However, Nd_2_CoMnO_6_ has similarities with La_2_CoMnO_6_, yet it has received less attention in research. It exhibits FM and meta-magnetic features,^19^ spin-phonon coupling,^20,21^ MD coupling,^15^ phase separation, and polaronic hopping conduction^22^ along with multiferroic behavior.^23^
Particularly, A_2_FeCrO_6_ (A = Pr, Bi) DPOs with two 3d TM at the B and B′ sites demonstrate the intriguing optoelectronic aspects as they display significant absorbance in the visible range of solar radiation and a favorable band gap energy (Eg).^24–26^ The FM or FIM behavior of the Cr-based ordered A_2_CrMO_6_ (A = Sr,Ca; M = Mo, W, Re) structures have also been reported, whereby Cr^3+^ (3d^3^, S = 3/2) and Mo^5+^ (4d^1^, S = 1/2), W^5+^, Re^5+^ (5d^2^, S = 1) paired antiferromagnetically.^27^ The Sr-based compounds such as Sr_2_CrMO_6_ (M = Mo and W) have high TC 473 and 453 K, respectively, according to their research, but the Ca-based Ca_2_CrMO_6_ (M = Mo/W) have a low TC of around 148/143 K.
Recently, the growing demand for thermoelectric (TE) materials, capable of recovering and converting waste heat into electrical energy at elevated temp., has led to increased interest in DPOs. These materials are favored due to their low thermal conductivity (k), high thermal stability, and environmentally friendly aspects compared to IVA-VIA TE alloys.^28–30^ However, the TE features have been investigated through both experimental and theoretical approaches for a limited number of DPOS, including La_2_CoMnO_6_,^17,31,32^ La_2_NiMnO_6_, Sr_2_BB'O_6_,^33^ and Pr_2_CoFeO_6_,^34^ among others. As the figure of merit (ZT) is a dimensionless parameter used to assess a TE material's overall performance. There are two main ways to increase ZT, improving the Seebeck coefficient (S) by optimizing the material's electronic structure and reducing the overall k, especially by suppressing the lattice k (kl).^35,36^ The chemical composition of the material has a significant impact on the kl, predominantly by means of lattice softening effects and phonon scattering processes [10]. Furthermore, the total number of atoms in the primitive unit cell has an important impact on determining its magnitude.^37^ In this context, Mustafa et al., investigated the halide DPOs and found that the ZT at 800 K for K_2_YAgBr_6_ and K_2_YAgI_6_ were high of 0.68 and 0.74, respectively.^38^
Motivated by the variety of applications of DPOs, a thorough analysis of the Lu_2_CoCrO_6_ (LCCO) has been performed, another possibility for multifunctional applications. To the best of our knowledge, its electronic, magnetic, and TE aspects have not yet been systematically explored, leaving a clear gap in the reported literature. Hence, examining LCCO allows us to understand how the presence of Co^3+^ and Cr^3+^ ions affects the superexchange mechanism and TE behavior, providing valuable comparative insight within the broader class of DPOs. Therefore, we comprehensively investigate its structural stability through formation enthalpy (ΔHf), mechanical, and dynamical behavior. The magnetic ordering is governed by superexchange interactions, together with electronic and magnetic aspects. To determine its potential in energy conversion technologies, we study its TE performance. The objective of this extensive investigation is to reveal the inherent physical features of the LCCO and provide the foundation for its implementation in TE and spintronic applications.
Computational and structural details
2
Spin-polarized (SP) density functional theory calculations were performed using the full potential linearized augmented plane wave (FP-LAPW) approach, which is implemented in the WIEN2k code.^39^ To effectively account for electron correlation effects, the generalized gradient approximation (GGA) is employed alongside an on-site Coulomb interaction (GGA + U), utilizing U value of 5.0/3.5 eV for Co/Cr-3d states.^40^ In the context of wavefunction expansion within atomic spheres, the angular momentum cutoff is established at lmax = 12. Additionally, the plane wave cutoff is specified by Rmt × Kmax = 7 with Gmax set to 24. The integration over the Brillouin zone is conducted by utilizing a 6 × 6 × 4 Monkhorst Pack k-mesh, resulting in 76 k-points within the irreducible wedge, thereby guaranteeing satisfactory convergence. Atomic positions are thoroughly optimized until the Hellmann–Feynman forces acting on all atoms decrease to below 5 mRy/a.u. The convergence criteria for achieving self-consistency are established at 10^−5^ Ry for total energy (Et) and charge density. The TE parameters are evaluated by applying Boltzmann transport theory under the constant relaxation time of 10^−14^ s approximation using the BoltzTraP code.^41^
The LCCO crystallizes in a monoclinic symmetry with the space group P2_1_/n (No. 14). The experimentally determined lattice parameters are a = 5.1416 Å, b = 5.4723 Å, c = 7.4158 Å, and β = 89.0211°.^42^ There are 4Lu, 2Co, 2Cr, and 12O atoms in its primitive unit cell. The atomic positions within the unit cell are defined as follows: Lu occupies the Wyckoff site 4e with coordinates (0.525, 0.571, 0.251); Co/Cr is located at 2d (0.5, 0, 0)/2c (0, 0.5, 0), respectively. The three distinct oxygen atoms O1, O2, and O3 occupy 2e positions at (0.391, 0.970, 0.248), (0.165, 0.154, −0.057), and (0.315, 0.683, −0.059), correspondingly.^42^ Experimentally, the distribution of Co and Cr ions at the B-site is found to be partially disordered. However, in the present study, we assume a perfectly ordered arrangement for computational simplicity. The LCCO crystal structure and spin magnetization density plots are depicted in Fig. 1(a) and (b), respectively.
Schematic illustration of the (a) crystal diagram and (b) spin magnetization density plot in a ferrimegnetic spin ordering for the Lu2CoCrO6 structure.
Results and discussion
3
Structural stabilities
3.1
As the Goldschmidt tolerance factor (τ) is often employed to anticipate the crystallographic phase and evaluate the structural stability of the oxides, which is computed as:^43,44^
in which rLu, rCo, rCr, and rO represent the ionic radii of these cations, having values of 1.03 Å, 0.55 Å, 0.62 Å, and 1.40 Å, respectively. As a result, τ ≈ 0.87, which falls within the standard stability range of 0.8 ≤ τ ≥ 1.0 and is significantly below 0.97, suggests a monoclinic DPO framework as opposed to an ideal cubic structure.^45^
Now, formation enthalpy (ΔHf) is calculated to analyze the thermodynamic stability of the structures as:
where , E^Lu-trigonal^t, E^Co-cubic^t, E^Cr-cubic^t, and represents the Et of the LCCO motif, Lu (R3̄m-166), Co (Fm3̄m-225), Cr (Fm3̄m-225), and O (C2/m-12) molecule in their respective ground states, correspondingly. The computed ΔHf is −4.22 eV per atom, where a negative sign confirms the thermodynamic stability of the structure.
Next, to evaluate the mechanical stability of the structure, we calculated the elastic tensors (Cij) through the application of six finite deformations and an analysis of the stress–strain response.^46,47^ In the monoclinic phase of the structure, we observe 13 independent Cij, which are in accordance with the necessary and sufficient Born criteria for the mechanical stability of the materials as follows:^48^ and
The calculated Cij (see matrix below) for the LCCO DPOs satisfy the specifications given above inequalities, thereby affirming their mechanical stability.
Furthermore, elastic constants play a significant role in examining the stiffness, brittleness, ductility, and anisotropy of the material using Bulk modulus (B)/Young's modulus (Y)/Shear modulus (G). The substantial hardness of the LCCO is indicated by its higher resistance to volume deformation under external pressure, as indicated by its B value of 260.36 GPa. Y = 298.66 GPa and G = 114.1 GPa validate its remarkable resistance to elastic deformation under uniaxial force, further demonstrating its intrinsic stiffness. Also, it is categorically classified as ductile by the Pugh's ratio and Poisson's ratio (ν = 0.308 > 0.26),^49^ both of which surpass the crucial criteria for brittleness. The ν and ratio can be used to clarify the material's bonding nature. The ratio of 0.438 indicates the ionic nature of the system, because ionic materials usually have a value of 0.6, which is closer to ionic bonding. While ν = 0.25 implies ionic bonding, and ν = 0.1 hints at covalent bonding.^49–52^ The computed ν of 0.30881 shows a dominance of ionic bonding, because of its near to 0.25.
Likewise, Fig. 2 depicts the 2D and 3D anisotropic representations of the elastic parameters for the LCCO system along the xy, xz, and yz planes, and their corresponding calculated values are listed in Table 1. The plotted surfaces strong directionality illustrates the compound's mechanical anisotropy. For an ideal isotropic case, these mentioned surfaces are perfectly spherically symmetric; however, the degree of variation from a perfect sphere indicates the degree of elastic anisotropy, is closely linked to structural symmetry and bonding properties. As Y (see Fig. 2(a)) displays moderate anisotropy having value of 1.45, refers that stiffness varies noteably with crystallographic orientation. In contrast, Fig. 2(b)–(d) represents that β/G/ν shows a stronger directional dependence in contrast to Y, demonstrating substantial variation in compressibility, shear resistance, and lateral strain between various planes. These results highlight that LCCO is elastically anisotropic with significant orientation-dependent stiffness and compressibility. Moreover, the bonding nature of the material is explain in terms of charge density map as displayed in Fig. 4. It is clear that charge density is found to be strongly localized around the O atoms, while very little electron density is observed between the metal atoms (Co, Cr, and Lu). The absence of shared electron clouds between the metal–metal sites confirms that there is no direct covalent bonding among them. Instead, the bonding is mainly ionic in character, with oxygen acting as the electron-rich center coordinating the metal cations. This supports and strengthens the conclusion drawn from ν.
Calculated 3D-directional dependent elastic parameters for the Lu2CoCrO6 structure.
Now, the phonon dispersion curve along the high-symmetry route (Γ Z D B Γ A E Z C Y Γ) is computed in Fig. 3 to access the dynamic stability of the LCCO system. Hence, the lack of imaginary phonon frequencies across the Brillouin zone confirms that the LCCO motif is dynamically stable.^53,54^ In principle, the number of atoms in the primitive unit cell determines the number of phonon branches. Each atom contributes three vibrational modes having a total of 3n branches, where n is the number of atoms. These consist of 3n − 3 optical and 3 acoustic modes.^55^ There are 20 atoms in the primitive unit cell, the phonon dispersion shows 60 branches with 57 optical and 3 acoustic. Acoustic modes are represented by the lower frequency portion of the dispersion curves, whereas optical modes are represented by the higher frequencies. Oxygen atoms are the source of the high frequency optical modes because of their low atomic mass, which is inversely correlated with vibrational frequency. On the other hand, heavier atoms like Lu, Co, and Cr provide a larger contribution to the low frequency acoustic modes.
Calculated phonon spectrum of the Lu2CoCrO6 structure.
Finally, the quasi-harmonic Debye model^56,57^ has been used to calculate the thermodynamic aspects in the temp. ranges of 0 to 1000 K at a constant pressure of 0/5/15/25 GPa. Fig. 5 displays the variation of volume (V), Debye temp. (θD), thermal expansion coefficient (α), entropy (S), specific heat capacity at constant volume (CV), and the specific heat capacity at constant pressure (CP). Due to compressive forces and thermal expansion, the unit cell volume decreases with pressure and increases with temp. as presented in Fig. 5(a).
Computed charge density plots of the Lu2CoCrO6 structure.
Computed (a) volume (V), (b) Debye temperature (θD), (c) absorption coefficient (α), (d) entropy (S), and (e) and (f) specific heat capacity at constant pressure/constant (CP/CV) volume of the Lu2CoCrO6 structure.
Along with this, Fig. 5(b) shows that θD decreases as temp. rises, a phenomenon linked to the softening of phonon modes at higher temp. On the other hand, when pressure increases, θD rises, indicating the lattice stiffens, resulting in a less interatomic space. Now, Fig. 5(c) displays how the temp. and pressure affects the α. Up to 300 K, a sharp increase in its value is observed with temp. and a slow or near-saturation point appears beyond 700 K. It shows that pressure has an inverse relation with α: as pressure rises, α falls rapidly. This reduction is attributed to the suppression of atomic displacements under compression, which limits the value of α. This finding demonstrates the material's thermal resilience under various thermodynamic circumstances and provides important information about its volumetric response to temp. Fig. 5(d) displays S, commonly referred to as disorder in systems. It demonstrates the variance of entropy in connection with particular pressures and temp. An increase in lattice vibration brought on by a rise in temp. may excite electrons, increasing the system's entropy further. Fig. 5(e) and (f) illustrates how CP/CV varies with respect to temp. and pressure, which serve as important indicators of phase stability and lattice dynamics. As can be seen, in the low-temp. region (0–400 K), CV/CP grows dramatically with temp., indicating the quick triggering of phonon modes. At temp. over 400 K, the rise in CV slows and approaches the conventional Dulong–Petit limit at 700 K, confirming harmonic behavior at high temp. Both CP and CV exhibit a modest drop in magnitude at any given temp. as pressure increases, which is explained by suppressed lattice vibrations brought on by volume contraction.
Magnetic aspects
3.2
To study the magnetic ground state, the structure was initially optimized by considering three spin-orderings (SO), such as FM, FIM, and antiferromagnetic (AFM) to precisely ascertain the equilibrium lattice parameters, volume, and Et at the ground state using the Murnaghan equation of state. The Et was computed as a function of volume for all three considered configurations as presented in Fig. 6, which confirms that the system is more stable in the FIM SO. The Co and Cr ions are anti-aligned (↑↓) in the FIM SO and remained aligned (↑↑) in the FM ones, and the individual ions of Co/Cr are aligned antiparallel (↑↓)/(↑↓) in AFM SO. The energy differences were calculated as ΔE1 = EFIM – EAFM = −31.54 meV and ΔE2 = EFIM – EFM = −186.56 meV. The negative signs indicate that the FIM SO is more stable than that of the FM and AFM ones, which is in agreement with the recent experimental findings.^42^
Calculated optimized energy-volume relationship for the Lu2CoCrO6 structure in the ferrimagnetic (FIM), ferromagnetic (FM), and antiferromagnetic (AFM) spin ordering.
To further examine the magnetic behavior of the system, the total/partial spin magnetic moment (mt/ms) is computed. The estimated mt is 2.00 µB/f.u. with the ms on the Co/Cr is 3.11/−2.52 µB. The “–” sign on the Cr ms signifies that the magnetic moments of Co and Cr are oriented antiparallel (↑↓) to one another, thereby affirming the existence of AFM coupling, which results in a net FIM SO within the system. Additionally, the computed ms of the Co/Cr ion indicates that they lie in a +3(t^3^2g↑t^1^2g↓e^2^g↑e^0^g↓)/+3(t^3^2g↑t^0^2g↓e^0^g↑e^0^g↓) electronic configuration. Similarly, the computed orbital moment (ml) for the Co/Cr is 0.103/0.042 due to the SOC effect. Further, the 3D spin magnetization density of the system is illustrated in Fig. 1(b), clearly demonstrating the AFM coupling between Co and Cr ions. The distinct colors of the isosurfaces validate their AFM coupling and the slightly more spatial distribution around the Co-site indicates its larger ms value in comparison to that of Cr. The observed asymmetry in the spin density magnitude corresponds with the calculated ms, highlighting the significant influence of the Co on the magnetic behavior of the compound. Further, the distinct separation of spin densities provides additional evidence for the AFM coupling, aligning with the anticipated orbital interactions in this DPO.
For more understanding about the interactions, the various forms of superexchange interactions between magnetic ions mediated by oxygen 2p orbitals are depicted in Fig. 7. The Co(e^2^g↑)–O–Co(e^2^g↑)/Cr(e^0^g)–O–Cr(e^0^g) interactions demonstrate FM coupling, characterized by the alignment of the Co/Cr ms in a ↑↑/↑↑orientation as shown in Fig. 7(a) and (b). Coversely, the charge transfer takes place between the partially filled Co^3+^–t^1^2g↑ and Cr^3+^–t^0^2g↑ orbitals, where electrons hop back and forth across the intermediate O^2−^ 2p states. In this way, a FIM ground state is stabilized by this indirect exchange process, which takes the form of an AFM superexchange via the Co–O–Cr route (see Fig. 7(c)).
Schematic illustration of the ferromagnetic (FM) superexchange interactions of (a and b) Co–O–Co/Cr–O–Cr and (c) antiferromagnetic (AFM) coupling between Co+3 3d6 and Cr+3 3d+3 ions via oxygen in the Lu2CoCrO6 structure.
Electronic structure
3.3
To elucidate the electronic structure of the LCCO motif, the computed SP total density of states (TDOS) and the atom-resolved partial density of states (PDOS) in the FIM SO via the GGA + U method are depicted in Fig. 8(a) and (b), correspondingly. Fig. 8(a) illustrates that the system is a SC with an Eg of 1.13/3.96 eV in the spin-minority/spin-majority channel (N↓)/(N↑). Hence, the material real Eg is 1.13 eV, which further affirms its SC nature.To acquire a deeper understanding of the states adjacent to the Fermi energy (EF), we plotted the Co/Cr-3d and O-2p orbital resolved PDOS in Fig. 8(b). A significant contribution close to (EF) at the edge of the valence band (VB) arises from the Cr-3d orbitals. In a similar manner, the edge of the conduction band (CB) is primarily characterized by significant hybridization between the 3d orbitals of the Cr, along with minor contributions from Co-3d states (see Fig. 8(b)). Additionally, the calculated SP band structures are presented in Fig. 9 along with the high symmetry directions of the monoclinic Brillouin zone, further displaying the SC behavior of the material. A direct Eg of 1.13/3.96 eV in the N↓/N↑ exists and corroborates the computed TDOS in Fig. 8(a). Furthermore, we plotted the GGA + U + SOC calculated TDOS/PDOS in Fig. 1S of the supporting information (SI) in the stable FIM SO. Our results indicate that SOC does not significantly affect the electronic structure of the material; however, a small enhancement of 0.08 eV in the Eg with the inclusion of SOC is revealed. Analysis of the PDOS hints that the primary contributions to the TDOS near the EF are attributed to Cr-3d orbitals (see Fig. 1S(b) of the SI). Moreover, here we would like to mentioned that why we adopted P2_1_/n structure of the LCCO, which corresponds to the ordered Co–Cr arrangement. Although the nominal valences reported there were Co^2+^/Cr^4+^, our DFT results yield Co^3+^ and Cr^3+^. This charge state is consistent with the ms balance in the system: Co^3+^ and Cr^3+^ provide the correct total spin configuration and stabilize the experimentally consistent magnetic ordering. Thus, the P2_1_/n structure is retained, while electron redistribution leads to a Co^3+^/Cr^3+^ state that achieves magnetic and electronic stability.
Calculated spin-polarized (a) total density of states (TDOS) and (b) partial density of states (PDOS), highlighting the role of Co-3d/Cr-3d/O-2p orbitals around the Fermi level for the Lu2CoCrO6 structure.
Calculated spin-polarized electronic band structures of the (a) spin-majority/(a′) spin-minority channel for the Lu2CoCrO6 structure.
Thermoelectric features
3.4
Finally, to study the TE aspects, first κl has been evaluated utilizing Slack's equation as:^58^
In this context, Mav represents the average atomic mass of the crystal, δ denotes the cube root of the average atomic volume, n indicates the number of atoms in the unit cell, T refers to the absolute temp., γ signifies the Grüneisen parameter, derived ν through the following equation:
The factor A and θD are determined as:^59,60^
where ℏ represents the Planck's constant, kB denotes the Boltzmann constant, NA refers to Avogadro's number, ρ indicates the density, M signifies the molecular weight, and vm stands for the mean sound velocity, which is:
where vl and vt denote the longitudinal and transverse sound velocities, respectively, defined as:
The analysis presented in Fig. 10 describes that the calculated κl decreases as temp. increases, which is an advantageous outcome for TE applications. The observed reduction in κl corresponds to increased phonon scattering at higher temp., resulting from increased lattice vibrations that decrease the phonon mean free path. The measured value of κl is 1.66/1.00/0.71 Wm^−1^ K^−1^s at 300/500/700 K, which has an impact on the ZT value at elevated temp.
Computed lattice thermal conductivity (kl) as a function of temperature of the Lu2CoCrO6 structure.
Furthermore, the semi-classical Boltzmann transport theory, which is implemented in the BoltzTraP2 code,^61^ was used to study the transport aspects of the LCCO in the rigid band approximation under the constant scattering time approximation. The evaluation of TE parameters like electrical conductivity per relaxation time , S, electronic κ per relaxation time , , and was conducted as a function of the chemical potential (µ) within the range of −3.5 to 3.5 Ry. In this context, positive values of µ indicate n-type carriers, while negative values correspond to p-type ones. The computed is shown in Fig. 11, demonstrating minimal variations with temp. The peak exhibits at 300 K having a value of 3.9 × 10^19^ (Ω ms)^−1^ in the n-type region and 2.7 × 10^19^ (Ω ms)^−1^ for the p-type coincide with the CB edge, indicating that the majority of the n-type carrier concentration. The S depends on the contributions of both n and p-type charge carriers. As shown in Fig. 11(b), S shows positive values for p-type carriers (E − EF < 0), having the highest value of 1.151 mV/K and negative values for n-type carriers (E − EF > 0) with 1.153 mV K^−1^ at 300 K, where a sign change at the EF signifying the transition between two conduction types. For n-type carriers, S reaches notably high absolute values at energies exceeding the EF. The magnitude of S typically decreases with temp. in both carrier types with in the given temp. range of 300–700 K. As predicted, rises significantly when µ moves into the conduction band (µ > 0), while S shows a peak close to the band edge at µ = 0.01 Ry and subsequently declines as µ increases. The inverse relationship between and S is a defining feature of the TE materials, resulting from the opposing dependence of these features on the carrier concentration.
Computed numerous thermoelectric features as a function of chemical potential (µ) of the Lu2CoCrO6 structure.
Next, the exhibits a similar qualitative trend to , showing maxima in areas of high carrier concentration and demonstrating moderate temp. dependence. The increase in with rising temp. is linked to the simultaneous enhancement in carrier mobility. It has the highest value of 5.19 × 10^14^ W mK^−1^ s^−1^ for n-type and 3.58 × 10^14^ W mK^−1^ s^−1^ for p-type region at 700 K, indicating n-type majority charge carriers. Now, the temp. dependence of the PF has been examined, revealing a linear increase in its value throughout the analyzed temp. range as displayed in Fig. 11(d). The peak is noted at 700 K owing the value of 2.20 × 10^11^ W mK^−2^ s^−1^ in the n-type and 1.06 × 10^11^ W mK^−2^ s^−1^ in the p-type region, indicating that LCCO attains enhanced TE performance at higher temp. Ultimately, the dimensionless ZT (see Fig. 11(e)), aligns with the trend of the PF while also being affected by the overall κe + κl. In both p and n-type regimes, ZT shows notable peaks near the band edges, aligning with the maxima in the PF. The magnitude of ZT increases from 300 to 700 K primarily because the κe rises with temperature and the lattice thermal conductivity (κl) decreases due to enhanced phonon–phonon scattering. Although S decreases at higher temperatures, the combined effect of improved charge carrier transport and suppressed κl results in a net enhancement of ZT.
Conclusion
4
In summary, thermodynamic, mechanical, dynamical, electronic, magnetic, and thermoelectric (TE) aspects of the Lu_2_CoCrO_6_ structure were analyzed using abinitio calculations. The energetic stability is confirmed by the computed negative formation enthalpy, and mechanical robustness is validated by the elastic constants that meet the Born criteria. Likewise, the dynamical stability of the system is confirmed by the phonon dispersion analysis, which reveals that no imaginary frequency modes across the Brillouin zone exist. Furthermore, the computed electronic structure calculations revealed that the system is a semiconductor having a direct band gap of 1.13 eV. Additionally, anti-ferromagnetic superexchange coupling is observed between Co^+3^ 3d^6^ and Cr^+3^ 3d^3^ ions via oxygen, confirms the ferrimagnetic spin ordering. This further affirms the 3D spin-magnitazion density plot and also the spin moments of Co (3.11 µB) and Cr (−2.52 µB) ions are align antiparallel. Also, the TE parameters varies as a function of chemical potential, such as low thermal conductivity, improved electrical conductivity, and a large Seebeck coefficient shows promising performance. Interestingly, a maximum figure of merit of 1.00 at a chemical potential of 0.1 eV under 700 K enhances its potential for TE application along with spintronics.
Author contributions
Samia Shahzadi: writing – original draft, investigations, formal analysis, data curation. A. Elfasakhany: research support and resources. S. Nazir: writing – review and editing, validation, supervision, project administration, conceptualization.
Conflicts of interest
The authors declare no competing interests.
Supplementary Material
RA-015-D5RA07005H-s001
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1Basavarajappa M. G. Chakraborty S. ACS Mater. Au 202226556643685554710.1021/acsmaterialsau.2c 00031 PMC 9928412 · doi ↗ · pubmed ↗
- 2Alam M. Ghosh L. Majumder S. Singh P. Kumar S. V. Dixit S. Kumar D. Anand K. Kumari S. Ghosh A. et al. J. Phys. D: Appl. Phys.202255255003
- 3Cheong S.-W. Mostovoy M. Nat. Mater.2007613201719912110.1038/nmat 1804 · doi ↗ · pubmed ↗
- 4Fiebig M. Lottermoser T. Fröhlich D. Goltsev A. V. Pisarev R. V. Nature 20024198188201239735210.1038/nature 01077 · doi ↗ · pubmed ↗
- 5Chen X. Xu J. Xu Y. Luo F. Du Y. Inorg. Chem. Front.2019622262238
- 6Mandal T. K. Felser C. Greenblatt M. Kübler J. Phys. Rev. B: Condens. Matter Mater. Phys.200878134431
- 7Feng H. L. Arai M. Matsushita Y. Tsujimoto Y. Guo Y. Sathish C. I. Wang X. Yuan Y.-H. Tanaka M. Yamaura K. J. Am. Chem. Soc.2014136332633292455216610.1021/ja 411713 q · doi ↗ · pubmed ↗
- 8Kobayashi K.-I. Kimura T. Sawada H. Terakura K. Tokura Y. Nature 1998395677680