Abstract

We present a comprehensive investigation of the electromagnetic properties of CdAB2 compounds, where A represents Cr, Mn, or Fe, and B denotes P or As. To investigate the spin-polarized behavior of these compounds the A atoms were substituted at the Group IV (Ge) position in CdGeB2 in the chalcopyrite crystal structure. Our results reveal that all the CdAB2 compounds exhibit compelling spin-splitting of energy states near the Fermi level (EF). Notably, CdAB2 materials with A = Cr and Mn exhibit intriguing half-metallic ferromagnetic (HMF) characteristics, with the calculated total magnetic moments of 2.00 and 3.00 µB/f.u., respectively. The HMF properties originated in CdAB2 (A = Cr and Mn; B = P, As) these compounds owing to the hybridization of partially filled -3d(t2g) states of A atoms with the p-states of B (P, As) atoms, with minor contributions from Cd’s-like states. In contrast, CdFeB2 displays distinct behavior, demonstrating spin-splitting of energy levels around the EF indicative of a stable ferromagnetic (FM) state and the absence of HMF at their equilibrium volume. The calculated total magnetic moments for CdFeP2 and CdFeAs2 are about 1.83 (1.64 µB/f.u.) and 1.94 µB/f.u. (1.84 µB/f.u.) under generalized gradient approximation (GGA) (local spin density approximation (LSDA)) approximations, respectively. Perhaps these CdAB2 compounds (A = Cr and Mn; B = P, As) with HMF characteristic within both LSDA and GGA formalisms makes them highly promising candidates for spin injectors in the spintronic device applications. Furthermore, their semiconducting nature renders CdCrB2 and CdMnB2 materials compatible with silicon and other semiconducting lattices, enhancing their potential practical applications in the spintronic technologies. In conclusion, this study presents a thorough exploration of the robust electronic and magnetic properties of CdAB2 chalcopyrites, offering exciting prospects for their utilization in the future spintronic applications.

1. Introduction

The quest for half-metallic ferromagnetic (HMF) materials with Curie temperature (TC) surpassing room temperature has garnered significant attention for the potential spin-electronic device applications [1]. Many materials have been explored for their HMF properties, such as Ti-doped ZnS [2], Mn-doped InAs and GaAs [3], V-doped BeB (B = S, Se, Te) [4], CrB (B = S, Se, Te) [5, 6], Mg1−xVxY (Y = Se, Te) [7], and Sc, Ti, V, Mn, and Cr-doped CdS [8]. These systems exhibit HMF behavior due to the substitution of transition metal (TM) ions, inducing small local magnetic moments in the nonmagnetic host materials with TC up to 200 K [3, 9]. Consequently, they have attracted great interest in spintronic applications [18]. However, their TC does not exceed 200 K [3], making them less favorable for applications in spin-based solid-state devices.

For practical spintronic applications, materials displaying room-temperature ferromagnetism (FM) are highly preferred. Ternary II–IV–V2 semiconductors, crystallizing in the body-centered tetragonal (BCT) chalcopyrite structure with the space Group I-42d, have emerged as promising candidates for various technological applications [1013]. By replacing a small fraction of constituent atoms with TM ions such as Cr or Mn, II–IV–V2 materials can be engineered for spintronic purposes. The doping of TM ions induces a strong spin splitting of energy levels around the Fermi level (EF) with relatively high TC [1014]. The appearance of HMF property relies on the concentration of magnetic ions and temperature. Experimental results have revealed room temperature FM (TC = 320 K) in Mn-doped CdGeP2 due to the high-carrier solubility of Mn atoms in CdGeP2, where Mn2+ easily occupies the divalent Cd site without altering electrical neutrality.

The substitutional site of TM ions in ternary compounds plays a crucial role in determining FM or antiferromagnetic (AFM) properties. When M (A = Mn, Cr, or V) atoms are substituted at the Group II (divalent) site of Ternary II–IV–V2 chalcopyrite, AFM spin arrangement is favored. Conversely, when the same M atoms are substituted at the tetravalent (Group IV) site, FM configuration with hole formation is observed [15]. Notably, stable FM states have been achieved with TC = 355 K by substituting Mn at the Ge site in CdGeP2, making the experimental fabrication of such materials feasible. Moreover, TM-doped or replaced chalcopyrite (II–IV–V2) semiconductors have been shown to exhibit FM above room temperature, drawing attention in the field of spintronics [13, 14, 16, 17].

Theoretical and experimental results have demonstrated room-temperature FM (TC = 320 and 340 K, respectively) in MnGeP2 and MnGeAs2 chalcopyrite materials [13, 1719], motivating us to conduct first-principles calculations on the structural, electronic, and magnetic properties of CdAB2 (A = Cr, Mn, Fe; B = P, As) chalcopyrite’s.

Building on the previous research, [19, 20] which explored spin-polarized calculations on ZnMX2 (A = Sc, Ti, V, Cr, Mn, Fe; X = P, As) and CdAB2 (A = Sc, Ti, V; X = P, As) compounds, revealing the possibility of FM upon M atom substitution at the Ge site of II(Zn/Cd)IV(Ge)–VI(P/As)2 compounds, our present work extends this study by substituting M atoms (A = Cr, Mn, Fe) at the Ge site of CdGeX2. Utilizing the full-potential linearized augmented plane wave (FP-LAPW) method with two correlation functionals, generalized gradient approximation (GGA) [21] and local spin density approximation (LSDA) [22, 23], we investigate the introduction of deep levels by TM impurities in the band gap of CdGeB2 semiconductors, resulting in attractive electronic and magnetic properties in these chalcopyrite’s. To understand the magnetic states of CdAB2 compounds, we perform total energy vs. volume calculations for NM, FM, and AFM phases. Additionally, we present spin-dependent band structures, charge density plots, and density of states (DOS) at equilibrium conditions for all the compounds.

This study contributes to the growing understanding of robust room temperature ferromagnetism in CdAB2 chalcopyrite through transition metal substitution, offering valuable insights for the potential spintronic applications.

2. Method of Calculation and Crystal Structure

In this study, we conducted first-principal calculations using the Full-Potential Linearized Augmented Plane Wave (FP-LAPW) method [22], known for its accuracy in electronic structure calculations of the periodic solids. The WIEN2k code [24] was employed for the calculations of the structural, electronic, and magnetic properties of CdAB2 compounds. Both GGA [24] and LSDA [22] were used for the exchange-correlation functions in these calculations.

The FP-LAPW technique employs a muffin-tin potential to represent the crystal potential. This assumes that the crystal potential is spherically symmetric within the muffin-tin sphere and remains constant in the interstitial space. The core electrons within the atomic sphere are treated with full relativistic considerations and self-consistency using a spherical approximation. Meanwhile, the valence electrons in the interstitial region are treated as self-consistently as possible, but with a semi-relativistic approach. Inside the atomic sphere, the crystal potential, electronic wave functions, and charge density are expanded using spherical harmonics up to a maximum angular momentum quantum number of lmax = 10. In the interstitial region, a Fourier series of plane waves is employed, with a cutoff defined as the product of the muffin-tin radius (RMT) and the maximum wave vector (Kmax), set to 7.0. For specific elements like Cd, Cr, Mn, Fe, P, and As, atomic sphere radii are designated as 2.08, 2.23, 2.24, 2.17, 2.20, and 2.09, respectively. The charge density is expanded in a Fourier series up to a maximum reciprocal lattice vector denoted as Gmax = 12. To integrate over k-points in the irreducible wedge of the first Brillouin zone, an 8 × 8 × 4 k-mesh is used. Through an iterative process, the crystal’s self-consistent total energy is brought to convergence, achieving a value lower than 0.01 m Ry.

CdGeB2 semiconductor adopts the chalcopyrite structure with BCT lattice type and space group I-42d (122). The atomic positions of the atoms are as follows: Cd1 (0, 0, 0), Cd2 (0, 1/2, 1/2); Ge1 (1/2, 1/2, 0); Ge2 (1/2, 0, 1/4), B1 (u, 1/4, 1/8); B2(u, 3/4, 1/8); B3 (3/4, u, 7/8); B4 (1/4, u, 7/8), where u represents the anion displacement (internal parameter). Ternary CdAB2 (A = Cr, Mn, Fe; B = P, As) compounds are derived from CdGeB2 by replacing Ge with 3d TM (Cr, Mn, Fe) ions. The ionic, atomic radii and electronegativities of M [25] atoms are compared to Ge, ensuring the retention of the chalcopyrite crystal structure for the all substituents (A = Cr, Mn, Fe) in CdGeB2 compounds. In the chalcopyrite structure, each B (P/As) anion is coordinated with two Cd and two Ge/M cations, while each cation is tetrahedrally coordinated by four B (P/As) anions. Slight displacement of the anions from their ideal tetrahedral site results in two unequal cation–anion bond lengths. Maintaining the tetragonal distortion ɳ = c/2a ∼ 1 as in CdG2, the chalcopyrite crystal symmetry is preserved in all CdAB2 compounds.

3. Results and Discussion

Before delving into electronic and magnetic properties, we performed full structural optimization for the CdAB2 (A = Cr, Mn, Fe; B = P, As) compounds. The chalcopyrite structure (BCT phase) with space group I-42d (# No. 122) was described using unit cell lattice parameters (ao, co) and one internal parameter (u). Structural optimization was carried out for all compounds within both the GGA and LSDA for the nonmagnetic (NM), ferromagnetic (FM), and antiferromagnetic (AFM) phases. Since the results obtained under GGA and LSDA schemes were found to be identical, the calculated total energy as a function of the lattice parameter for GGA is represented in Figure 1.

To determine the equilibrium values of u, ao, co, and bulk modulus (Bo) of CdAB2, the total energies were fitted to Birch–Murnaghan’s equation of state [26] as a function of relative volume within GGA and LSDA for the NM, FM, and AFM states of each compound. The calculated values of u, ao, co, and Bo for CdGeB2 and CdAB2 compounds are presented in Tables 1 and 2. Remarkably, the calculated ground state properties of CdGeB2 are in excellent agreement with the experimental and previously reported results [27]. However, it is noted that the calculated lattice constants are overestimated by GGA and underestimated by LSDA, as expected. The internal parameter u was found to vary from 0.268 to 0.306, as observed from the tables.

Figure 1 showcases the total energy–volume curves of CdAP2 and CdAAs2 (A = Cr, Mn, Fe) compounds under GGA approximation. The energy–volume curves for CdAB2 (A = Cr, Mn, Fe; B = P, As) compounds indicate that the FM phase exhibits the minimum energy, signifying their ferromagnetic stability at equilibrium. To further assess the stability of the magnetic state, we calculated the spin-polarized energy differences ΔE1 (ΔE1 = EFM −ENM) and ΔE2 (ΔE2 = EFM −EAFM), i.e., the energy gains due to the formation of magnetic moments in the systems. The negative values of total energy differences ΔE1 and ΔE2 for CdAB2 (M = Cr, Mn, Fe) compounds indicate the desirability of the FM state in these materials. The calculated total energy differences ΔE1 and ΔE2 under GGA and LSDA can be found in Tables 1 and 2.

The equilibrium bulk modulus (Bo) values of CdAB2 (A = Cr, Mn, Fe; B = P, As) compounds were calculated for their NM, FM, and AFM phases using Equation (1) [28]. The calculated Bo values are presented in Tables 1 and 2 under the GGA and LSDA schemes.

Bulk modulus is a crucial mechanical property that characterizes a material’s resistance to volume changes during expansion or compression. With the substitution of transition metal (TM) atoms in the host CdGeB2 compounds, the lattice parameters (ao, co) may increase or decrease, leading to corresponding changes in the Bo values due to the displacement of the electric charge of core electrons via valence electrons. Notably, GGA tends to yield larger lattice parameters and smaller bulk modulus values compared to the LSDA.

The chalcopyrite crystal structure is a superlattice of the zinc blende (ZnS) structure. Doubling the unit cell of ZnS results in a change in the anion As/P position along the z-direction. This leads to tetragonal distortion in CdAB2 compounds, forming different bond lengths. The relaxed anion–cation bond lengths (RCd-X and RM-X) for the host CdGeB2 and CdAB2 compounds were calculated for the NM, FM, and AFM states using Equations (2) and (3) [29] under GGA and LSDA and are presented in Tables 3 and 4.

In CdAB2 compounds with M = Cr and Mn, the calculated bond length RM-B is shorter than RCd-B due to the larger electronegativity of B [22] atoms. Consequently, the M atoms move closer to the B atom, resulting in decreased RM-B in these compounds with u values > 0.25. Conversely, for u < 0.25, RCd-B becomes reduced compared to RFe-B in CdFeB2, where RFe-B increases due to the decrease in the electro-positivity of Fe atoms. In this case, the Cd atom moves toward the B atom. It is important to highlight that bond lengths exhibit a reduction in the GGA method as opposed to the LSDA approach. This reduction in bond lengths intensifies the rigidity among Cd, M, and B atoms and amplifies the electron concentration in the valence area. In order to enhance our comprehension of electron behavior within the valence region, additional calculations were performed to determine the compounds’ band structure, density of states (DOS), and electron charge density.

Additionally, the heat of formation (ΔH) for all the CdAB2 compounds in the NM, FM, and AFM phases were calculated using Equation (4) [29], and the results are provided in Tables 3 and 4.

The computed ΔH (enthalpy) values were discovered to be negative using both the LSDA and GGA methodologies. This signifies that all the compounds are expected to be thermodynamically stable in the chalcopyrite structure. These findings provide valuable insights into the mechanical and thermodynamic stability of CdAB2 chalcopyrite, contributing to a deeper understanding of their potential applications in various fields of materials science.

The spin-resolved electronic band structures of CdAB2 compounds (where A = Cr, Mn, Fe, and B = P, As) were computed along high-symmetry directions within the first Brillouin zone at their equilibrium volume using both the GGA and LSDA. Electronic band structures for CdAB2 compounds (A = Cr, Fe; B = P, As) are depicted in Figures 2 and 3 for the GGA approach. These figures illustrate the overall profile of the band structure, which remains quite similar for all CdAB2 compounds under both approximations. The introduction of -d states near the Fermi level (EF) due to the substitution of elements A = Sc, Ti, V, Cr, Mn, and Fe in the ZnGeB2 structure leads to a crystal field resonance. This resonance causes a splitting of the -d states into suborbitals labeled as “eg” and “t2g”. As a result of the hybridization between -d states and -p states from the nearest neighboring B atom, an indirect band gap emerges in the spin-down channel. Meanwhile, the, eg and t2g states within the majority spin channel intersect the Fermi level, indicating a HMF characteristic within these compounds.

Figure 2 illustrates the spin-polarized band structure of CdCrP2, highlighting the noticeable spin polarization of energy states in the vicinity of the Fermi level (EF). Specifically, the bonding -3d, eg states and the anti-bonding t2g states associated with the Cr atom is notably concentrated near the EF. This arrangement results in the formation of an indirect band gap between the valence band maximum (VBM) and the conduction band minimum (CBM) within the spin-down channel. In contrast, in the spin-up channel, the eg, t2g states of the -3d Cr atom and -3p states of P atoms extend and cross the EF, leading to HMF nature in this compound. Similar phenomena were observed in CdCrAs2 and CdMnB2 (B = P, As) compounds. The calculated HMF gap (EHM) and spin-splitting energy gaps (Eg↓) of CdAB2 (A = Cr, Mn; B = P, As) are listed in Table 5 for both GGA and LSDA.

On the other hand, CdFeP2 (Figure 3) shows spin-splitting of energy levels around the EF, indicating the occurrence of stable ferromagnetism. However, the VBM crosses the EF in both the spin-up and spin-down channels, signifying the absence of HMF property in these compounds. The decrease in exchange spin-splitting energy is attributed to the coupling of magnetic t2g and eg suborbital of -3d states of Fe atoms and the repulsion between the -3d (t2g and eg) states of Fe and -p states of P atoms in the VBM. This leads to an increase in the VB width, reducing the localization of the electron in the VB region and increasing the spin-polarization around the EF.

The charge density contour plots (Figures 4 and 5) offer insight into the chemical bond nature in CdAP2 (A = Cr, Fe) compounds. The plots provide evidence of charge transfer involving valence electrons shifting from Cd to P and from M to P ions. This transfer signifies a blend of both ionic and covalent attributes in the atomic bonds between Cd-P and M-P, holding for both spin channels across all compounds. An examination of the electronic charge distribution surrounding the P atom reveals the establishment of Cd-P and M-P bonds. Particularly strong p-d exchange interactions between Cr (d) and P (p) states are evident, driving a substantial covalent character between P and Cr within CdCrP2. Likewise, the compounds CdCrAs2 and CdMnB2 (where B = P, As) showcase pronounced charge transfer between M and P ions. This results in a robust p-d(t2g) interaction between the band carriers and localized spins, leading to a HMF nature in these compounds.

Total and partial density of states (PDOS) calculations (Figures 6 and 7) for CdCrP2 and CdFeP2 within GGA reveal significant changes in the DOS near the EF, leading to spin-polarization of the valence bands (VBs). In CdCrP2, the eg orbitals are well-localized and form an indirect band gap between the VB and CB in the spin-down channel. In the spin-up channel, the t2g states cross the EF, leading to HMF behavior. A similar phenomenon was observed in CdCrAs2 and CdMnB2 (B = P, As), where the p-d (t2g) states of B and A (Cr, Mn) atoms play a dominant role in the formation of the gap between valence and conduction bands in the spin-down channel.

Conversely, CdFeP2 and CdFeAs2 do not exhibit HMF behavior, as the strong overlapping and repulsive reaction between -3d states of Fe (eg and t2g) and -p states of P/As reduce the magnetic moment of Fe, resulting in spin-splitting of energy states around the EF without an observed band gap in the minority spin channel.

Using the spin-polarized calculations, the total and local magnetic moments of CdAB2 compounds (where A = Cr, Mn, Fe and B = P, As) have been computed within both the LSDA and GGA methods. These results are summarized in Table 5. It is evident from the findings that the primary source of the total magnetic moments stems from the partially filled -3d states of the M atoms. For instance, in CdCrB2, the Cr atom contributes some of its valence electrons to the majority spin band during bond formation, resulting in a total magnetic moment of 2.00 µB/f.u., as shown in Table 5.

Similarly, the Mn atom in CdMnB2 contributes three unpaired d electrons for a total magnetic moment of 3.00 µB per formula unit. However, in CdFeB2, the Fe atom’s magnetic moment gets reduced due to hybridization and repulsion, resulting in stable ferromagnetism without HMF behavior. The magnetic moments calculated by GGA and LSDA significantly differ, with GGA yielding higher values. Our obtained spin-resolved electronic band structures and charge density analysis provide valuable insights into the electronic and magnetic properties of CdAB2 chalcopyrites, revealing their potential for applications in various fields of materials science.

4. Conclusion

In pursuit of new HMF materials, we conducted comprehensive first-principles calculations on CdAB2 compounds in the chalcopyrite structure within the GGA and LSDA schemes. Ground state properties were evaluated from the energy–volume relation. The CdAB2 compounds were found to be ferromagnetically stable from the total energy calculations. The structural parameters of CdGeB2 host compounds agreed well with existing experimental and theoretical data, no comparative data were available for CdAB2 compounds. The negative formation energies (ΔH) indicated that these compounds are energetically favorable when A occupies the Group IV (Ge) position in CdGeB2.

The spin-polarized calculations unveiled intriguing spin-splitting of energy states near the Fermi level (EF) upon the substitution of A atoms at the Group IV (Ge) position in CdGeB2. For CdAB2 compounds with A = Cr and Mn, a distinct energy gap was observed around the EF in the spin-down channel, while the VB crossed the EF in the spin-up channel, confirming the presence of HMF nature with 100% spin polarization around the EF. Therefore, the HMF property in CdAB2 (A = Cr, Mn; B = P, As) arose due to the hybridization of A -3d (t2g) and B -3/4p states, with a minor contribution from Cd -s states and their corresponding calculated total magnetic moment were found to be 2.00 µB and 3.00 µB per formula unit, respectively.

Conversely, in CdFeB2 compounds, the overlapping of wave functions of -3d (t2g and eg) states of Fe and repulsion with B -np states were more pronounced, leading to a reduction in the space charge density of the Fe atom. This effect diminished the possibility of HMF due to the broadening of VBs. The calculated total magnetic moments for CdFeP2 and CdFeAs2 were determined to be 1.83 (1.64 µB/f.u.) and 1.94 µB/f.u. (1.84 µB/f.u.), respectively, under GGA (LSDA) approximations. Consequently, the ground state of CdFeB2 compounds did not exhibit significant HMF characteristics when using LSDA and GGA functionals. However, considerable magnetic properties, such as stable ferromagnetism, were observed in these compounds. In conclusion, our first-principles study sheds light on the potential HMF nature in CdAB2 (A = Cr, Mn; B = P, As) chalcopyrites. These findings hold significance in the ongoing search for novel HMF materials, offering promising avenues for the future spintronics applications.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.