- 1Department of Physics, Faculty of Sciences, Dr. Tahar Moulay University of Saida, Saida, Algeria
- 2Instrumentation and Advanced Materials Laboratory, University Center of Nour Bachir El-Bayadh, El-Bayadh, Algeria
- 3Technology Department, University Center of Nour Bachir El-Bayadh, El-Bayadh, Algeria
- 4Modelling and Simulation in Materials Science Laboratory, Physics Department, Djillali Liabes University of Sidi Bel-Abbes, Sidi Bel-Abbes, Algeria
- 5University of Artois, CNRS, Centrale Lille, Univ. Lille, UMR 8181 – UCCS – Unité de Catalyse et Chimie du Solide, Lens, France
The first-principles computations of density functional theory are employed to characterize the structural properties, electronic structures, and ferromagnetism induced by Cr impurities in Ca1-xCrxO compounds at concentrations x = 0. 25, 0.5, and 0.75. The dynamic stability is performed by the phonon spectra calculations. The structural parameters are computed by using Wu-Cohen generalized gradient approximation, while the electronic and magnetic properties are determined by the accurate Tran–Blaha-modified Becke–Johnson exchange potential. The crystal field, direct and indirect exchange splittings were investigated to determine the origin and stability of ferromagnetic state configuration. The Ca1-xCrxO systems have right half-metallicities, which are verified by the spin polarization of 100% and the integer values of total magnetic moments. The Ca0.75Cr0.25O, Ca0.5Cr0.5O, and Ca0.25Cr0.75O are half-metallic ferromagnetic with flip-gaps of 1.495, 0.888, and 0.218 eV, respectively. Therefore, the Ca1-xCrxO materials are suitable candidates for possible applications of spin-injection in future semiconductors spintronics.
Introduction
In recent years, the exploitation of new materials for technological development has become an attractive area of research for the scientific community. Spintronics is a modern field emerging from quantum electronics relating to exploration and control of the electron spin; it is currently an important research topic because the spin-based devices have interesting multifunctional characteristics such as non-volatility, data-processing speed, higher integration densities (Elilarassi and Chandrasekaran, 2017; Rai et al., 2018). In the last years, the half-metal ferromagnetic (HMF) materials have attracted considerable interest due to their wide application in spintronic devices as magnetic sensors, tunnel junctions, and spin injection sources of electrons (Özdemir and Merdan, 2019). The diluted magnetic semiconductors (DMSs) are obtained from the III–V, II–VI, or IV–VI semiconductors doped with magnetic elements, inducing changes in their structural, electronic, and magnetic properties, which are related to the concentration of magnetic impurities. These materials have received a lot of attention for scientific research because they are promising candidates for innovative applications in spintronics, as an example for non-volatile memories, optical switches, and magneto-optical devices (Ohno, 1998; Dietl et al., 2000; Wolf et al., 2001; Zuti et al., 2004; Sato et al., 2010; Kumar et al., 2014).
Alkaline earth chalcogenides are important materials with a wide technological application ranging from catalysis to microelectronics as well as in the area of electro- and photo-luminescent devices (Hakamata et al., 2005; Speziale et al., 2006; Murtaza et al., 2014). Among these compounds, the CaX (X = O, S, Se) has a closed-shell ionic environment, crystallizing in the rock-salt NaCl-type (B1) structure at ambient conditions (Doll et al., 1996; Mishra et al., 2012). Under high pressure, they undergo a first-order structural phase transition from the NaCl-type (B1) to the CsCl-type (B2) crystal structure (Luo et al., 1994; Doll et al., 1996; Dadsetani and Doosti, 2009; Mishra et al., 2012). The CaO calcium oxide can be a potential candidate as a DMS material according to several theoretical studies (Kenmochi et al., 2004; Kemmochi et al., 2005; An Dinh et al., 2006). Kenmochi et al. (2004) have studied a new class of DMSs based on CaO without transition metals, where the substituted CaO by C and N at the oxygen sites shows room-temperature ferromagnetism with half-metallic property. The ferromagnetism in C-doped CaO is explained by the exchange coupling constants, where it is induced by the strong hybridization between 2p states of carbon and 2p states of oxygen (An Dinh et al., 2006). Further, Kemmochi et al. (2005) have found that the stability of ferromagnetic state in the C-doped CaO is originated from the Zener's double-exchange mechanism.
In this study, we have investigated the structural parameters, electronic performance, magnetism, and exchange splittings of Ca1-xCrxO materials at various concentrations x = 0, 0.25, 0.5, and 0.75. The calculations are performed by the use of first-principle approaches of density functional theory (DFT) (Hohenberg and Kohn, 1964; Kohn and Sham, 1965) as implemented in the WIEN2k package (Blaha et al., 2001). The phonon spectra are used to verify the dynamic stability of Ca1-xCrxO materials. The exchange and correlation potential is evaluated by the generalized gradient approximation of Wu and Cohen (GGA-WC) (Wu and Cohen, 2006), and the Tran–Blaha-modified Becke–Johnson potential (TB-mBJ) (Tran and Blaha, 2009).
Computational Methods
The calculations of structural, electronic, and magnetic properties of Ca1-xCrxO are carried out using the WIEN2k code (Blaha et al., 2001) based on the first-principle computations of DFT and the full-potential linearized augmented plane-wave (FP-LAPW) method (Singh and Nordstrom, 2006). We have computed the structural parameters using GGA-WC approximation (Wu and Cohen, 2006), while the electronic and magnetic properties are determined by employing the TB-mBJ exchange potential (Becke and Johnson, 2006; Tran and Blaha, 2009). The phonon frequencies are determined by the phonopy code (Togo and Tanaka, 2015) based on the density functional perturbation theory (Giannozzi et al., 1991; Gonze and Lee, 1997). The phonon spectra are calculated using the WIEN2k and the phonopy packages to verify the dynamic stability of Ca1-xCrxO doping compounds.
The averages of Muffin-tin spheres radii are selected for Ca, Cr, and O atoms in condition that Muffin-tin spheres do not overlap. The cutoff of −6 Ry is utilized to separate between the valence and core states. The basic functions and the potentials are extended in combination of spherical harmonics around the atomic sites, that is to say, the atomic spheres with a cutoff lmax = 10, and in the interstitial region are extended into plane waves with a cutoff RMTkmax = 8 (where RMT is the average radius of the Muffin-tin spheres). We have expanded the charge density in Fourier up to Gmax = 14 (a.u.)−1, where Gmax represents the largest vector in the Fourier expansion. The calculations were performed for structures of eight (8) atoms of Ca1-xCrxO at concentrations x = 0, 0.25, 0.5, and 0.75 such as CaO, Ca0.75Cr0.25O, Ca0.5Cr0.5O, and Ca0.25Cr0.75O supercells. The integration of the Brillouin zone was carried out by the Monkhorst–Pack mesh (Monkhorst and Pack, 1976) using the special k-points of (10 × 10 × 10) for CaO and (5 × 5 × 5) for Ca0.75Cr0.25O, Ca0.5Cr0.5O, and Ca0.25Cr0.75O compounds. The self-consistent is achieved when the total energy convergence condition was set to be 10−1 mRy.
Results and Discussion
Structural Properties
Crystal Structures
The calcium oxide CaO is known as a wide-band gap II–VI semiconductor, crystallizing in rock-salt NaCl (B1) structure at ambient conditions (Doll et al., 1996; Mishra et al., 2012). The primitive unit cell of CaO contains two atoms such as the calcium (Ca) and the oxygen (O), which belong, respectively, to the main groups II and VI of the periodic table. Our study is based on the conventional NaCl structure of eight atoms of CaO such as Ca4O4, where the Ca is located at (0, 0, 0) position and the O atom at (0.5, 0.5, 0.5) with space group of No. 225. The Ca3CrO4, Ca2Cr2O4, and CaCr3O4 supercells are created by substituting one, two, and three Cr impurities at Ca sites, respectively. Consequently, we have obtained three supercells such as Ca0.75Cr0.25O at concentration x = 0.25, Ca0.5Cr0.5O at x = 0.5 and Ca0.25Cr0.75O at x = 0.75. The Ca0.75Cr0.25O and Ca0.25Cr0.75O have cubic lattices with space group of No 221, while the Ca0.5Cr0.5O supercell has a tetragonal lattice with space group of P4/mmm No. 123. We have noticed that all symmetry operations have been considered in the initialization of structures of Ca1-xCrxO supercells for computations of different properties.
Structural Stability
The structural stability of the Ca1-xCrxO compounds in the ferromagnetic rock-salt NaCl (B1) structure is verified by computing the formation energies. The formation energies (Eform) of the Ca4–yCryO4 materials are determined using the following expression (Bai et al., 2011; Doumi et al., 2015a):
where Etotal(Ca4–yCryO4) corresponds to the minimum total energy of Ca4–yCryO4 per atom, and E(Ca), E(Cr), and E(O) refer to the minimum total energies per atom of bulks Ca, Cr, and O, respectively. The y = 1, 2, and 3 are the number of substituted Cr impurities in Ca4–yCryO4 supercells. We have found that the formation energies are −5.99, −6.37, and −6.76 eV, respectively, for Ca0.75Cr0.25O, Ca0.5Cr0.5O, and Ca0.25Cr0.75O. Therefore, the negative values of formation energies suggest that these compounds are thermodynamically stable in the ferromagnetic rock-salt structure.
Furthermore, we have used the phonon spectra to demonstrate the dynamic stability of the Ca1-xCrxO doping compounds. The diagrams of phonon frequencies at high symmetry points in the Brillouin zone of Ca0.75Cr0.25O, Ca0.5Cr0.5O, and Ca0.25Cr0.75O are shown by Figures 1–3, respectively. These figures depicted that all branches of phonons have positive frequencies, and there are no negative frequencies (Deng et al., 2019), confirming the dynamic stability of Ca0.75Cr0.25O, Ca0.5Cr0.5O, and Ca0.25Cr0.75O compounds.
Equilibrium Structural Parameters
The variations of total energies as a function of volumes of Ca1-xCrxO compounds are fitted by Murnaghan (1944) equation of state to determine the equilibrium structural parameters. Table 1 summarizes the lattice constants, the bulk modules, and their first derivatives with other experimental data (Mammone et al., 1981; Kaneko et al., 1982; Richet et al., 1988) and theoretical values (Tran et al., 2007; Wu et al., 2014; Cinthia et al., 2015; Fan et al., 2015; Nejatipour and Dadsetani, 2015; Santana et al., 2016; Yang et al., 2016; Salam, 2018) computed by the generalized gradient approximation of Perdew–Burke–Ernzerhof (GGA-PBE) (Perdew et al., 1996), revised GGA of Perdew–Burke–Ernzerhof (GGA-PBEsol) (Perdew et al., 2008) and the local density approximation (LDA) (Perdew and Zunger, 1981; Perdew and Wang, 1992).
Table 1. Computed structural parameters such as lattice constants (a), bulk modules (B), and their pressure derivatives (B′) for Ca1-xCrxO materials at concentrations x = 0, 0.25, 0.5, and 0.75 with other experimental and theoretical data.
The results of lattice constants and bulk modulus of CaO are closed to theoretical values of Tran et al. (2007) and Fan et al. (2015) calculated with GGA-WC potential (Wu and Cohen, 2006). We see that there is a good agreement between our calculations and experimental results (Mammone et al., 1981; Kaneko et al., 1982; Richet et al., 1988) and theoretical value found by GGA-PBE sol (Perdew et al., 2008). On the other hand, the obtained values of these parameters with GGA-WC (Wu and Cohen, 2006) are better than the calculations of Tran et al. (2007), Wu et al. (2014) Cinthia et al. (2015), Nejatipour and Dadsetani (2015), Yang et al. (2016), Salam (2018), Tran et al. (2007), Fan et al. (2015), and Santana et al. (2016) found, respectively, by GGA-PBE (Perdew et al., 1996) and LDA approximation (Perdew and Zunger, 1981; Perdew and Wang, 1992). The performance of GGA-WC approximation for predicting structural parameters is due to the fourth-order gradient expansion of exchange-correlation functional (Doumi et al., 2015b; Sajjad et al., 2015). For the Ca1-xCrxO doping compounds, the ionic radius of Cr is lower than that of Ca atom, leading to decrease in the lattice constant with increasing Cr concentration. Consequently, the bulk modulus of Ca1-xCrxO increases with increasing concentration of Cr. The Ca1-xCrxO doping compound becomes harder as the concentration of Cr increases.
Electronic Properties
Band Structures and Half-Metallicity
We have used the predicted structural parameters to characterize the electronic structures and the half-metallicity in Ca1-xCrxO compounds. The TB-mBJ potential is used to determine the electronic structures and accurate band gaps of Ca1-xCrxO compounds at concentration x = 0, 0.25, 0.5, and 0.75. The band structures of CaO, Ca0.75Cr0.25O, Ca0.5Cr0.5O, and Ca0.25Cr0.75O are shown, respectively, by the Figures 4–7.
Figure 4. Band structures of CaO. (A) Majority spins (up) and (B) minority spins (dn). The Fermi level is set to zero (horizontal dotted line).
According to Figure 4, the CaO is a semiconductor because both majority and minority spins revealed the analogous band structures with indirect gaps. The CaO has an indirect gap(EΓX) between the valence bands maximum and conduction bands minimum situated, respectively, at Γ and X high symmetry points. Figures 5–7 demonstrate that the Ca1-xCrxO for all concentrations are metallic and semiconductors, respectively, for majority- and minority-spin bands. Consequently, the Ca1-xCrxO doping materials are half-metallic ferromagnets. The metallic character is due to localized partially filled 3d-Cr states of majority spins around Fermi level (EF), which interact with carriers of valence states of host semiconductor. For minority spins, the Ca1-xCrxO doping systems keep a semiconductor behavior as CaO compound because the 3d-Cr empty minority-spin states arise in the bottom of conduction band minimum at Γ high symmetry point far than EF. Thus, the gap nature shifted from indirect band gap for CaO to direct gap for Ca1-xCrxO doping systems.
Figure 5. Band structures of Ca0.75Cr0.25O. (A) Majority spins (up) and (B) minority spins (dn). The Fermi level is set to zero (horizontal dotted line).
Figure 6. Band structures of Ca0.5Cr0.5O. (A) Majority spins (up) and (B) minority spins (dn). The Fermi level is set to zero (horizontal dotted line).
Figure 7. Band structures of Ca0.25Cr0.75O. (A) Majority spins (up) and (B) minority spins (dn). The Fermi level is set to zero (horizontal dotted line).
The minority spins of Ca1-xCrxO systems are characterized by direct wide-band gaps called half-metallic ferromagnetic (HMF) gaps located at the Γ high symmetry point. The minimal energies that separate the Fermi level (zero eV), and the conduction band minimums determine the half-metallic (HM) gaps. Table 2 shows the computed band gaps such as the direct HMF gaps (EHMF) and HM gaps (EHM) of Ca0.75Cr0.25O, Ca0.5Cr0.5O, and Ca0.25Cr0.75O, and the indirect gap (EΓX) of CaO with other theoretical (Fan et al., 2015; Nejatipour and Dadsetani, 2015; Yang et al., 2016; Tran and Blaha, 2017; Salam, 2018) and experimental data (Whited et al., 1973). The result of wide-indirect gap of 5.392 eV of CaO is in good agreement with recent calculated value of 5.35 eV of Tran and Blaha (2017) study using the same TB-mBJ potential (Becke and Johnson, 2006; Tran and Blaha, 2009), and it is better than the values ranging from 3.437 to 3.67 eV of Fan et al. (2015), Nejatipour and Dadsetani (2015), Yang et al. (2016), Tran and Blaha (2017), Salam (2018) found by the LDA (Perdew and Zunger, 1981; Perdew and Wang, 1992), GGA-WC (Wu and Cohen, 2006), GGA-PBE (Perdew et al., 1996), and GGA-PBEsol exchange potentials (Perdew et al., 2008). The TB-mBJ provides accurate gap for CaO because he is known for its performance in the calculation of electronic structures of insulators and semiconductors with respect to LDA and all forms of GGA approximations (Bhattacharjee and Chattopadhyaya, 2017a,b; Chattopadhyaya and Bhattacharjee, 2017; Berriah et al., 2018).
Table 2. Computed band gaps such as indirect band gap (EΓx) for CaO, half-metallic ferromagnetic (HMF) gaps (EHMF), and half-metallic (HM) gaps (EHM) of minority-spin bands for Ca1-xCrxO at concentrations x = 0.25, 0.5, and 0.75 with other theoretical and experimental data.
Moreover, the half-metallic gap or flip-gap is an interesting performance for exploring the spin-injection for a half-metallic material in spintronics. The Ca0.75Cr0.25O, Ca0.5Cr0.5O, and Ca0.25Cr0.75O compounds have the HM gaps of 1.495, 0.888, and 0.218 eV, respectively, situated between the conduction bands minimum (CBM) and EF. The CBM moves toward EF due to widening of 3d-Cr empty states at the bottom of conduction band, and hence, the HM gap decreases from Ca0.75Cr0.25O, Ca0.5Cr0.5O, to Ca0.25Cr0.75O. The HM gap behavior is maintained for different concentrations, meaning that the Ca1-xCrxO doping systems are right half-metallic materials.
Densities of States and Spin Polarization
In order to understand the origin of half-metallic feature induced by the substituting effect of Cr, we have plotted the total and partial densities of states (DOS) of Cr, Ca, and O atoms of Ca1-xCrxO compounds. Figures 8–10 show, respectively, the DOS of Ca0.75Cr0.25O, Ca0.5Cr0.5O, and Ca0.25Cr0.75O. We clearly see that all compounds are metallic for majority spin and semiconductors for minority spin, where their bonding, non-bonding, and anti-bonding states are distinguished from each other by the contributions of the DOS formed by different orbitals of each Ca, O, and Cr atoms.
Figure 8. Spin-polarized total and partial densities of states (DOS) of Ca0.75Cr0.25O. The Fermi level is set to zero (vertical dotted line). (A) Total DOS of Ca, Cr, and O and (B) Partial DOS of p-Ca, p-O, d-Cr, t2g_d-Cr, and eg_d-Cr.
Figure 9. Spin-polarized total and partial densities of states (DOS) of Ca0.5Cr0.5O. The Fermi level is set to zero (vertical dotted line). (A) Total DOS of Ca, Cr, and O and (B) DOS of p-Ca, p-O, d-Cr, t2g_d-Cr, and eg_d-Cr.
Figure 10. Spin-polarized total and partial densities of states (DOS) of Ca0.25Cr0.75O. The Fermi level is set to zero (vertical dotted line). (A) Total DOS of Ca, Cr, and O and (B) Partial DOS of p-Ca, p-O, d-Cr, t2g_d-Cr, and eg_d-Cr.
The partial DOS of Figures 8–10B revealed that the 3d-Cr states are splitted into the low-lying t2g and high-lying eg states owing to the effect of octahedral crystal field created by neighboring oxygen (O) ions. For both minority and majority spins, the bonding states are mainly contributed by the p-O states and by the small p-Ca and 3d-Cr states in the ranges of −6.5 to −4, −7.5 to −4.5, and −8.25 to −5.25 eV, respectively, for Ca0.75Cr0.25O, Ca0.5Cr0.5O, and Ca0.25Cr0.75O. For the majority spins of Ca1-xCrxO, the non-bonding states are principally dominated by the t2g-(3d-Cr) fully filled levels that hybridize largely with p-O orbitals, while the anti-bonding states are mostly contributed by the hybridization of the eg-(3d-Cr) levels and p-O states around EF. However, the upper part of valence bands of majority spins is originated from both t2g-(3d-Cr) non-bonding states and eg-(3d-Cr) anti-bonding states, which extend strongly in the gap owing to broadening of 3d-Cr levels around EF.
Furthermore, the upper part of valence bands of minority spins are formed by the bonding states originated from the p-O orbitals, which are centered at −5.25, −6, and −6.75 eV for Ca0.75Cr0.25O, Ca0.5Cr0.5O, and Ca0.25Cr0.75O, respectively. The bottom of conduction band is generated predominantly by the 3d-Cr empty anti-bonding states, which move to lower energies toward Fermi level from Ca0.75Cr0.25O, Ca0.5Cr0.5O, to Ca0.25Cr0.75O with increasing Cr concentration. Therefore, the distinguish half-metallic ferromagnetic and half-metallic gaps are created around EF in minority spins for Ca0.75Cr0.25O, Ca0.5Cr0.5O, and Ca0.25Cr0.75O materials.
The spin-polarized material is known by its spin polarization of electronic structures, which is an important factor for characterizing a half-metal. The polarization (P) of half-metal compound results from the different contributions of numbers of DOS around EF. It is determined from the following expression (Soulen et al., 1998; Wang et al., 2017):
where N ↑ (EF) and N ↓ (EF) indicate, respectively, the DOS of majority and minority spins at EF. From Figures 8–10A, the density of states N ↓ (EF) equals zero because the DOS do not cross EF, leading to a spin polarization of 100%. Therefore, the Ca1-xCrxO compounds appear to be suitable candidates for possible spin-injection in spintronics applications.
Magnetic Properties
Magnetic Moments
The spin-polarized electronic structures show an asymmetric character, indicating that the Ca1-xCrxO compounds are magnetic in nature. The magnetism is induced by the magnetic spins arising from the localized 3d-Cr majority-spin states around EF. The Cr impurity contributes two electrons to host carriers of valence bands, resulting in Cr+2 ion. Thus, the 3d-Cr states become partially filled with four electrons. Therefore, the magnetism is originated from the four itinerant electrons of Cr+2 ions in Ca1-xCrxO doping materials.
From the DOS of majority spins, one can see that the non-bonding states enclose three low-lying t2g levels, which are completely occupied by three electrons. In contrast, the two high-lying eg states have an anti-bonding nature, which are partially occupied by one electron because they are crossed by Fermi level in the middle. According to the Hund's rule, the 3d-Cr states generate a total magnetic moment of 4 μB per Cr ion. Consequently, the Ca0.75Cr0.25O, Ca0.5Cr0.5O, and Ca0.25Cr0.75O compounds have total magnetic moments of 4, 8, and 12 μB, respectively, for one, two, and three Cr doping impurities.
The calculated total and partial magnetic moments of Ca, Cr, and O atoms and in the interstitial sites for Ca0.75Cr0.25O, Ca0.5Cr0.5O, and Ca0.25Cr0.75O are given in Table 3. The total magnetic moment per Cr atom is an integer number of 4 μB. It is a typical feature of half-metallic materials (Wang et al., 2017; Zhang et al., 2017). All compounds have negative partial magnetic moments for O atoms, revealing that carriers of host valence bands enclosing p-O states interact anti-ferromagnetically with 3d-Cr states. Besides, the positive partial magnetic moments of Cr and Ca atoms describe the ferromagnetic interaction between Cr and Ca magnetic spins.
Table 3. Computed total and partial magnetic moments of the relevant Cr, Ca, and O atoms and in the interstitial sites (in Bohr magneton μB) for Ca1-xCrxO at concentrations x = 0.25, 0.5, and 0.75.
Exchange Splitting Mechanisms
The ferromagnetism is created from the localized electrons of t2g and e2g of 3d-Cr states mediated by holes, where the empty holes of eg states occur at the top of valence band and generate an acceptor hole type carrier hybridizing with carriers of host valence band (Doumi et al., 2013). The degenerate 3d-Cr into t2g and e2g is due to crystal field generated from the electrostatic environment of surrounding oxygen (O) ions, where the 3d-Cr majority-spin levels contain four unpaired electrons (Noor et al., 2018). According to the Hund's rule, this kind of process mediates magnetism in Ca1-xCrxO materials (Yaqoob et al., 2019). It is understood that both exchange splitting of 3d-Cr states and crystal field take part in the contribution of ferromagnetism in Ca1-xCrxO doping materials.
The magnitude of crystal field energy (ΔECF) is defined as a difference between Eeg and Et2g energies of eg and t2g states (ΔECF = Eeg − Et2g). The direct exchange splitting energy Δx(d) = d↓ − d↑ is determined from the separation between the empty d-Cr (d↓) minority-spin and occupied d-Cr (d↑) majority-spin peaks. Another factor is used to determine the magnitude of ferromagnetism such as the indirect exchange splitting resulted from the energy difference between and valence band maximums of minority- and majority-spin channels, respectively.
The values of ΔECF, Δx(d), and Δx(pd) exchange splitting are shown in Table 4. The obtained values of (ΔECF and Δx(d)) are (1.71 and 4.98 eV), (2.97 and 4.90 eV), and (1.74 and 5.17 eV), revealing that the direct exchange splittings are more dominated compared to the crystal field splittings. Therefore, the ferromagnetism state is mainly favored by the contribution of direct exchange splitting mechanism than that of crystal field (Amin et al., 2018). Besides, the partial magnetic moment of Cr increases from 3.652 μB for Ca0.75Cr0.25O, 7.203 μB for Ca0.5Cr0.5O, to 10.541 μB for Ca0.25Cr0.75O with decreasing Δx(pd), meaning that magnetic moment of Cr is directly related to the indirect exchange splitting. On the other hand, the attraction nature of spin-polarized electrons is measured from the indirect exchange splitting (Verma et al., 2011; Mahmood et al., 2016, 2018). The indirect exchange splittings Δx(pd) have negative values because the valence band maximums of minority-spins are situated at lower energies with respect to the valence band maximums of majority-spins. This means that the minority-spins have more attractive (negative) potential than that of majority-spins, leading to perfectly localized magnetic states. Consequently, the ferromagnetism is preferred in the Ca1-xCrxO compounds. In the DMS based on III–V and II–VI semiconductors doped with transition metals, when the d states of magnetic dopants are partially occupied (Sato and Katayama-Yoshida, 2002; Sato et al., 2003a,b), the stability of ferromagnetism is described by the double-exchange mechanism (Akai, 1998). According to this rule, the anti-bonding states of partially filled 3d-Cr atoms associated with double-exchange mechanism stabilize the ferromagnetic state configuration in Ca1-xCrxO.
Table 4. Computed energies of crystal field ΔECF, direct Δx(d), and indirect Δx(pd)exchange splittings for Ca1-xCrxO at concentrations x = 0.25, 0.5, and 0.75.
The magnitudes of valence and conduction band edge splitting are measured by the two sp-d exchange constants N0α and N0β. The N0α parameter determines the s-d exchange interaction between the s carriers of conduction bands and d-Cr levels, while the N0β describes the p-d exchange coupling between the p states of carriers of valence bands and d-Cr levels. We have used the mean field theory (Sanvito et al., 2001; Raebiger et al., 2004) to calculate the exchange constants of Ca1-xCrxO compounds given by the following relations.
The and difference energies correspond, respectively, to the valence band-edge and conduction band-edge spin-splittings at Γ high symmetry point. The 〈s〉 and x are the half total magnetic moment per Cr and the concentration of Cr, respectively (Sanvito et al., 2001). The calculated ΔEv, ΔEc, N0α, and N0β are summarized in Table 5. Yaqoob et al. (2019) have predicted that the negative values of N0α and N0β are due to the quantum confinement effect or structural symmetry changes that modify the coupling of different states. For Ca1-xCrxO materials, the negative values of N0α and N0β describe the anti-ferromagnetic coupling between the valence and conduction bands and the 3d states of magnetic Cr atoms.
Table 5. Computed conduction band-edge ΔEc(eV) and valence band-edge ΔEv(eV) spin-splittings, and sp-d exchange constants N0α and N0β for Ca1-xCrxO at concentrations x = 0.25, 0.5, and 0.75.
Conclusions
We have characterized the structural properties, the spin-polarized electronic structures, and ferromagnetic performance induced by the Cr impurities in Ca1-xCrxO compounds at concentrations x = 0, 0.25, 0.5, and 0.75. The Ca1-xCrxO doping materials are dynamically stable due to their positive phonon frequencies. The structural parameters and indirect gap of CaO are very significant with respect to the recent calculations found by GGA-WC and TB-mBJ. For Ca1-xCrxO, the lattice constant decreases with increasing Cr concentration, leading to the increase in bulk modulus, where the sizes of the ionic radii of Ca and Cr atoms are the factors determining these changes. The Ca1-xCrxO systems have a half-metallic ferromagnetic feature, where the ferromagnetism is favored by the large contribution of direct exchange splitting than that of the crystal field. The analysis of magnetic properties shows that both double-exchange and p-d exchange mechanisms participate to stabilize ferromagnetic state configuration. The half-metallic ferromagnets Ca1-xCrxO compounds maintain a half-metal gap feature for all concentrations with spin polarization of 100%, making them suitable spin-injection candidates for possible exploration in semiconductors spintronics.
Data Availability Statement
All datasets presented in this study are included in the article/supplementary material.
Author Contributions
With the submission of this manuscript, all authors declare that: The manuscript is original. All authors of this research paper have contributed to this manuscript. All authors of this paper have read and approved the final version submitted. The contents of this manuscript have not been copyrighted or published previously. The contents of this manuscript are not now under consideration for publication elsewhere.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
Akai, H. (1998). Ferromagnetism and its stability in the diluted magnetic semiconductor (In, Mn) As. Phys. Rev. Lett. 81, 3002–3005. doi: 10.1103/PhysRevLett.81.3002
Amin, B., Majid, F., Saddique, M. B., Haq, B. U., Laref, A., Alrebdi, T. A., et al. (2018). Physical properties of half-metallic AMnO3 (A= Mg, Ca) oxides via ab initio calculations. Comp. Mater. Sci. 146, 248–254. doi: 10.1016/j.commatsci.2018.01.033
An Dinh, V., Toyoda, M., Sato, K., and Katayama-Yoshida, H. (2006). Exchange interaction and Tc in alkaline-earth-metal-oxide-basedDMS without magnetic impurities: first principle pseudo-SIC and monte carlo calculation. J. Phys. Soc. Jpn. 75:093705. doi: 10.1143/JPSJ.75.093705
Bai, J., Raulot, J.-M., Zhang, Y., Esling, C., Zhao, X., and Zuo, L. (2011). The effects of alloying element Co on Ni–Mn–Ga ferromagnetic shape memory alloys from first-principles calculations. Appl. Phys. Lett. 98:164103. doi: 10.1063/1.3582239
Becke, A. D., and Johnson, E. R. (2006). A simple effective potential for exchange. J. Chem. Phys. 124:221101. doi: 10.1063/1.2213970
Berriah, K., Doumi, B., Mokaddem, A., Elkeurti, M., Sayede, A., Tadjer, A., et al. (2018). Theoretical investigation of electronic performance, half-metallicity, and magnetic properties of Cr-substituted BaTe. J. Comput. Electron. 17, 909–919. doi: 10.1007/s10825-018-1192-y
Bhattacharjee, R., and Chattopadhyaya, S. (2017a). Effects of barium (Ba) doping on structural, electronic and optical properties of binary strontium chalcogenide semiconductor compounds-A theoretical investigation using DFT based FP-LAPW approach. Mater. Chem. Phys. 199, 295–312. doi: 10.1016/j.matchemphys.2017.06.057
Bhattacharjee, R., and Chattopadhyaya, S. (2017b). Effects of doping of calcium atom (s) on structural, electronic and optical properties of binary strontium chalcogenides-A theoretical investigation using DFT based FP-LAPW methodology. Solid State Sci. 71, 92–110. doi: 10.1016/j.solidstatesciences.2017.06.010
Blaha, P., Schwarz, K., Madsen, G. K., Kvasnicka, D., and Luitz, J. (2001). In WIEN2k, An Augmented Plane Wave + Local Orbitals Program for Calculating Crystal Properties. ed K. Schwarz (Wien: Techn. Universitaet Wien).
Chattopadhyaya, S., and Bhattacharjee, R. (2017). Theoretical study of structural, electronic and optical properties of BaxPb1−xS, BaxPb1−xSe and BaxPb1−xTe ternary alloys using FP-LAPW approach. J. Alloy. Compd. 694, 1348–1364. doi: 10.1016/j.jallcom.2016.10.096
Cinthia, A. J., Priyanga, G. S., Rajeswarapalanichamy, R., and Iyakutti, K. (2015). Structural, electronic and mechanical properties of alkaline earth metal oxides MO (M= Be, Mg, Ca, Sr, Ba). J. Phys. Chem. Solids 79, 23–42. doi: 10.1016/j.jpcs.2014.10.021
Dadsetani, M., and Doosti, H. (2009). The linear optical properties for NaCl phase of calcium mono chalcogenides by density functional theory. Comp. Mater. Sci. 45, 315–320. doi: 10.1016/j.commatsci.2008.10.003
Deng, J., Liu, N., Guo, J., and Chen, X. (2019). Large spin gaps in the half-metals MN4 (M= Mn, Fe, Co) with N2 dimers. Phys. Rev. B 99:184409. doi: 10.1103/PhysRevB.99.184409
Dietl, T., Ohno, H., Matsukura, F., Cibert, J., and Ferrand, E. D. (2000). Zener model description of ferromagnetism in zinc-blende magnetic semiconductors. Science 287, 1019–1022. doi: 10.1126/science.287.5455.1019
Doll, K., Dolg, M., and Stoll, H. (1996). Correlation effects in MgO and CaO: cohesive energies and lattice constants. Phys. Rev. B 54:13529. doi: 10.1103/PhysRevB.54.13529
Doumi, B., Mokaddem, A., Ishak-Boushaki, M., and Bensaid, D. (2015a). First-principle investigation of magnetic and electronic properties of vanadium-and chromium-doped cubic aluminum phosphide. Mat. Sci. Semicon. Proc. 32, 166–171. doi: 10.1016/j.mssp.2015.01.014
Doumi, B., Mokaddem, A., Sayede, A., Dahmane, F., Mogulkoc, Y., and Tadjer, A. (2015b). First-principles investigations on ferromagnetic behaviour of Be1−xVxZ (Z= S, Se and Te)(x= 0.25). Superlattice Microst. 88, 139–149. doi: 10.1016/j.spmi.2015.09.008
Doumi, B., Tadjer, A., Dahmane, F., Mesri, D., and Aourag, H. (2013). Investigations of structural, electronic, and half-metallic ferromagnetic properties in (Al, Ga, In)1−xMxN (M= Fe, Mn) diluted magnetic semiconductors. J. Supercond. Novel. Magn. 26, 515–525 doi: 10.1007/s10948-012-1808-6
Elilarassi, R., and Chandrasekaran, G. (2017). Optical, electrical and ferromagnetic studies of ZnO: Fe diluted magnetic semiconductor nanoparticles for spintronic applications. Spectrochim. Acta A 186, 120–131. doi: 10.1016/j.saa.2017.05.065
Fan, Q., Chai, C., Wei, Q., Yang, Y., Qiao, L., Zhao, Y., et al. (2015). Mechanical and electronic properties of Ca1−xMgxO alloys. Mat. Sci. Semicon. Proc. 40, 676–684. doi: 10.1016/j.mssp.2015.07.035
Giannozzi, P., De Gironcoli, S., Pavone, P., and Baroni, S. (1991). Ab initio calculation of phonon dispersions in semiconductors. Phys. Rev. B 43:7231. doi: 10.1103/PhysRevB.43.7231
Gonze, X., and Lee, C. (1997). Dynamical matrices, Born effective charges, dielectric permittivity tensors, and interatomic force constants from density-functional perturbation theory. Phys. Rev. B 55:10355. doi: 10.1103/PhysRevB.55.10355
Hakamata, S., Ehara, M., Kominami, H., Nakanishi, Y., and Hatanaka, Y. (2005). Preparation of CaS: Cu, F thin-film electroluminescent devices with an emission including purple region. Appl. Surf. Sci. 244, 469–472. doi: 10.1016/j.apsusc.2004.10.104
Hohenberg, P., and Kohn, W. (1964). Inhomogeneous electron gas. Phys. Rev. 136:B864. doi: 10.1103/PhysRev.136.B864
Kaneko, Y., Morimoto, K., and Koda, T. (1982). Optical properties of alkaline-earth chalcogenides. I. Single crystal growth and infrared reflection spectra due to optical phonons. J. Phys. Soc. Jpn. 51, 2247–2254. doi: 10.1143/JPSJ.51.2247
Kemmochi, K., Seike, M., Sato, K., Yanase, A., and Katayama-Yoshida, H. (2005). New class of high-TC diluted ferromagnetic semiconductors based on CaO without transition metal elements. J. Supercond. 18, 37–40. doi: 10.1007/s10948-005-2147-7
Kenmochi, K., Seike, M., Sato, K., Yanase, A., and Katayama-Yoshida, H. (2004). New class of diluted ferromagnetic semiconductors based on CaO without transition metal elements. Jpn. J. Appl. Phys. 43, L934–L936. doi: 10.1143/JJAP.43.L934
Kohn, W., and Sham, L. J. (1965). Self-consistent equations including exchange and correlation effects. Phys. Rev. 140, A1133–A1138. doi: 10.1103/PhysRev.140.A1133
Kumar, S., Negi, N., Katyal, S., Sharma, P., and Sharma, V. (2014). Structural, morphological and magnetic analysis of Cd–Co–S dilute magnetic semiconductor nanofilms. J. Magn. Magn. Mater. 367, 1–8. doi: 10.1016/j.jmmm.2014.04.065
Luo, H., Greene, R. G., Ghandehari, K., Li, T., and Ruoff, A. L. (1994). Structural phase transformations and the equations of state of calcium chalcogenides at high pressure. Phys. Rev. B 50, 16232–16237. doi: 10.1103/PhysRevB.50.16232
Mahmood, Q., Alay-E-Abbas, S., Hassan, M., and Noor, N. (2016). First-principles evaluation of Co-doped ZnS and ZnSe ferromagnetic semiconductors. J. Alloy. Compd. 688, 899–907. doi: 10.1016/j.jallcom.2016.07.302
Mahmood, Q., Ali, S., Hassan, M., and Laref, A. (2018). First principles study of ferromagnetism, optical and thermoelectric behaviours of AVO3 (A= Ca, Sr, Ba) perovskites. Mater. Chem. Phys. 211, 428–437. doi: 10.1016/j.matchemphys.2018.03.019
Mammone, J., Mao, H., and Bell, P. (1981). Equations of state of CaO under static pressure conditions. Geophys. Res. Lett. 8, 140–142. doi: 10.1029/GL008i002p00140
Mishra, M., Sharma, G., Kothari, R., Vijay, Y., and Sharma, B. (2012). Electronic structure of CaX (X= O, S, Se) compounds using Compton spectroscopy. Comp. Mater. Sci. 51, 340–346. doi: 10.1016/j.commatsci.2011.07.015
Monkhorst, H. J., and Pack, J. D. (1976). Special points for Brillonin-zone integrations. Phys. Rev. B 13, 5188–5192. doi: 10.1103/PhysRevB.13.5188
Murnaghan, F. (1944). The compressibility of media under extreme pressures. Proc. Natl. Acad. Sci. U.S.A. 30, 244–247. doi: 10.1073/pnas.30.9.244
Murtaza, G., Khenata, R., Safeer, A., Alahmed, Z., and Omran, S. B. (2014). Shift of band gap from indirect to direct and optical response of CaO by doping S, Se, Te. Comp. Mater. Sci. 91, 43–49. doi: 10.1016/j.commatsci.2014.04.039
Nejatipour, H., and Dadsetani, M. (2015). Excitonic effects in the optical properties of alkaline earth chalcogenides from first-principles calculations. Phys. Scripta. 90:085802. doi: 10.1088/0031-8949/90/8/085802
Noor, N., Saddique, M. B., Haq, B. U., Laref, A., and Rashid, M. (2018). Investigations of half-metallic ferromagnetism and thermoelectric properties of cubic XCrO3 (X= Ca, Sr, Ba) compounds via first-principles approaches. Phys. Lett. A 382, 3095–3102. doi: 10.1016/j.physleta.2018.07.045
Ohno, H. (1998). Making nonmagnetic semiconductors ferromagnetic. Science 281, 951–956. doi: 10.1126/science.281.5379.951
Özdemir, E. G., and Merdan, Z. (2019). First-principles calculations on half-metal ferromagnetic results of VZrAs and VZrSb half-heusler compounds and Al1−xMxAs (M= Co, Fe and x= 0.0625, 0.125, 0.25) diluted magnetic semiconductors. J. Alloy. Compd. 807:151656. doi: 10.1016/j.jallcom.2019.151656
Perdew, J. P., Burke, K., and Ernzerhof, M. (1996). Generalized gradient approximation made simple. Phys. Rev. Lett. 77, 3865–3868. doi: 10.1103/PhysRevLett.77.3865
Perdew, J. P., Ruzsinszky, A., Csonka, G. I., Vydrov, O. A., Scuseria, G. E., Constantin, L. A., et al. (2008). Restoring the density-gradient expansion for exchange in solids and surfaces. Phys. Rev. Lett. 100:136406. doi: 10.1103/PhysRevLett.100.136406
Perdew, J. P., and Wang, Y. (1992). Accurate and simple analytic representation of the electron-gas correlation energy. Phys. Rev. B 45, 13244–13249. doi: 10.1103/PhysRevB.45.13244
Perdew, J. P., and Zunger, A. (1981). Self-interaction correction to density-functional approximations for many-electron systems. Phys. Rev. B 23, 5048–5079. doi: 10.1103/PhysRevB.23.5048
Raebiger, H., Ayuela, A., and Nieminen, R. (2004). Intrinsic hole localization mechanism in magnetic semiconductors. J. Phys-Condens. Mat. 16, L457–L462. doi: 10.1088/0953-8984/16/41/L05
Rai, D., Laref, A., Shankar, A., Sakhya, A. P., Khenata, R., and Thapa, R. (2018). Spin-induced transition metal (TM) doped SnO2 a dilute magnetic semiconductor (DMS): a first principles study. J. Phys. Chem. Solids 120, 104–108. doi: 10.1016/j.jpcs.2018.04.006
Richet, P., Mao, H. K., and Bell, P. M. (1988). Static compression and equation of state of CaO to 1.35 Mbar. J. Geophys. Res. 93, 15279–15288. doi: 10.1029/JB093iB12p15279
Sajjad, M., Manzoor, S., Zhang, H., Noor, N., Alay-E-Abbas, S., Shaukat, A., et al. (2015). The half-metallic ferromagnetism character in Be1−xVxY (Y= Se and Te) alloys: an ab-initio study. J. Magn. Magn. Mater. 379, 63–73. doi: 10.1016/j.jmmm.2014.11.004
Salam, M. M. A. (2018). Theoretical study of CaO, CaS and CaSe via first-principles calculations. Results Phys. 10, 934–945. doi: 10.1016/j.rinp.2018.07.042
Santana, J. A., Krogel, J. T., Kent, P. R., and Reboredo, F. A. (2016). Cohesive energy and structural parameters of binary oxides of groups IIA and IIIB from diffusion quantum Monte Carlo. J. Chem. Phys. 144:174707. doi: 10.1063/1.4947569
Sanvito, S., Ordejón, P., and Hill, N. A. (2001). First-principles study of the origin and nature of ferromagnetism in Ga1−xMnxAs. Phys. Rev. B 63:165206. doi: 10.1103/PhysRevB.63.165206
Sato, K., Bergqvist, L., Kudrnovsk,ý, J., Dederichs, P. H., Eriksson, O., Turek, I., et al. (2010). First-principles theory of dilute magnetic semiconductors. Rev. Mod. Phys. 82, 1633–1690. doi: 10.1103/RevModPhys.82.1633
Sato, K., Dederichs, P., Araki, K., and Katayama-Yoshida, H. (2003a). Ab initio materials design and Curie temperature of GaN-based ferromagnetic semiconductors. Phys. Status Solidi, 2855–2859. doi: 10.1002/pssc.200303289
Sato, K., and Katayama-Yoshida, H. (2002). First principles materials design for semiconductor spintronics. Semicond. Sci. Tech. 17, 367–376. doi: 10.1088/0268-1242/17/4/309
Sato, K., Katayama-Yoshida, H., and Dederichs, P. (2003b). Curie temperatures of III–V diluted magnetic semiconductors calculated from first-principles in mean field approximation. J. Supercond. 16, 31–35. doi: 10.1023/A:1023264113451
Singh, D. J., and Nordstrom, L. (2006). Planewaves, Pseudopotentials, and the LAPW Method. New York, NY: Springer Science and Business Media.
Soulen, R., Byers, J., Osofsky, M., Nadgorny, B., Ambrose, T., Cheng, S., et al. (1998). Measuring the spin polarization of a metal with a superconducting point contact. Science 282, 85–88. doi: 10.1126/science.282.5386.85
Speziale, S., Shieh, S. R., and Duffy, T. S. (2006). High-pressure elasticity of calcium oxide: a comparison between Brillouin spectroscopy and radial X-ray diffraction. J. Geophys. Res. 111:B02203. doi: 10.1029/2005JB003823
Togo, A., and Tanaka, I. (2015). First principles phonon calculations in materials science. Scripta Mater. 108, 1–5. doi: 10.1016/j.scriptamat.2015.07.021
Tran, F., and Blaha, P. (2009). Accurate band gaps of semiconductors and insulators with a semilocal exchange-correlation potential. Phys. Rev. Lett. 102:226401. doi: 10.1103/PhysRevLett.102.226401
Tran, F., and Blaha, P. (2017). Importance of the kinetic energy density for band gap calculations in solids with density functional theory. J. Phys. Chem. A 121, 3318–3325. doi: 10.1021/acs.jpca.7b02882
Tran, F., Laskowski, R., Blaha, P., and Schwarz, K. (2007). Performance on molecules, surfaces, and solids of the Wu-Cohen GGA exchange-correlation energy functional. Phys. Rev. B 75:115131. doi: 10.1103/PhysRevB.75.115131
Verma, U., Sharma, S., Devi, N., Bisht, P., and Rajaram, P. (2011). Spin-polarized structural, electronic and magnetic properties of diluted magnetic semiconductors Cd1−xMnxTe in zinc blende phase. J. Magn. Magn. Mater. 323, 394–399. doi: 10.1016/j.jmmm.2010.09.016
Wang, X., Khachai, H., Khenata, R., Yuan, H., Wang, L., Wang, W., et al. (2017). Structural, electronic, magnetic, half-metallic, mechanical, and thermodynamic properties of the quaternary Heusler compound FeCrRuSi: a first-principles study. Sci. Rep. Uk 7:16183. doi: 10.1038/s41598-017-16324-2
Whited, R., Flaten, C. J., and Walker, W. (1973). Exciton thermoreflectance of MgO and CaO. Solid. State. Commun. 13, 1903–1905. doi: 10.1016/0038-1098(73)90754-0
Wolf, S., Awschalom, D., Buhrman, R., Daughton, J., Von Molnar, S., Roukes, M., et al. (2001). Spintronics: a spin-based electronics vision for the future. Science 294, 1488–1495. doi: 10.1126/science.1065389
Wu, X., Liu, L., Li, W., Wang, R., and Liu, Q. (2014). Effect of temperature on elastic constants, generalized stacking fault energy and dislocation cores in MgO and CaO. Comp. Cond. Matter. 1, 38–44. doi: 10.1016/j.cocom.2014.10.005
Wu, Z., and Cohen, R. E. (2006). More accurate generalized gradient approximation for solids. Phys. Rev. B 73:235116. doi: 10.1103/PhysRevB.73.235116
Yang, X., Wang, Y., Yan, H., and Chen, Y. (2016). Effects of epitaxial strains on spontaneous polarizations and band gaps of alkaline-earth-metal oxides MO (M= Mg, Ca, Sr, Ba). Comp. Mater. Sci. 121, 61–66. doi: 10.1016/j.commatsci.2016.04.021
Yaqoob, N., Sabir, B., Murtaza, G., Khalil, R. M. A., Muhammad, N., and Laref, A. (2019). Structural, electronic, magnetic, optical and thermoelectric response of half-metallic AMnTe2 (A= Li, Na, K): an ab-initio calculations. Physica B 574:311656. doi: 10.1016/j.physb.2019.08.033
Zhang, L., Wang, X., and Cheng, Z. (2017). Electronic, magnetic, mechanical, half-metallic and highly dispersive zero-gap half-metallic properties of rare-earth-element-based quaternary Heusler compounds. J. Alloy. Compd. 718, 63–74. doi: 10.1016/j.jallcom.2017.05.116
Keywords: first-principles calculations, electronic structures, ferromagnetism, exchange splittings, spintronics
Citation: Doumi B, Mokaddem A, Tadjer A and Sayede A (2020) Recent Insights Into Electronic Performance, Magnetism and Exchange Splittings in the Cr-substituted CaO. Front. Chem. 8:526. doi: 10.3389/fchem.2020.00526
Received: 04 April 2020; Accepted: 22 May 2020;
Published: 10 July 2020.
Edited by:
Xiaotian Wang, Southwest University, ChinaCopyright © 2020 Doumi, Mokaddem, Tadjer and Sayede. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Adlane Sayede, adlane.sayede@univ-artois.fr