Next Article in Journal
Rare-Earth Activated Nitride Phosphors: Synthesis, Luminescence and Applications
Next Article in Special Issue
Mn-doped Ge and Si: A Review of the Experimental Status
Previous Article in Journal
Electrochemical Light-Emitting Gel
Previous Article in Special Issue
Ti-doped ZnO Thin Films Prepared at Different Ambient Conditions: Electronic Structures and Magnetic Properties
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Carrier States in Ferromagnetic Semiconductors and Diluted Magnetic Semiconductors—Coherent Potential Approach—

Kanagawa Institute of Technology, 1030 Shimo-Ogino, Atsugi, 243-0292, Japan
Submission received: 5 May 2010 / Accepted: 8 June 2010 / Published: 21 June 2010
(This article belongs to the Special Issue Ferromagnetic Semiconductors)

Abstract

:
The theoretical study of magnetic semiconductors using the dynamical coherent potential approximation (dynamical CPA) is briefly reviewed. First, we give the results for ferromagnetic semiconductors (FMSs) such as EuO and EuS by applying the dynamical CPA to the s-f model. Next, applying the dynamical CPA to a simple model for A 1 x Mn x B-type diluted magnetic semiconductors (DMSs), we show the results for three typical cases to clarify the nature and properties of the carrier states in DMSs. On the basis of this model, we discuss the difference in the optical band edges between II-V DMSs and III-V-based DMSs, and show that two types of ferromagnetism can occur in DMSs when carriers are introduced. The carrier-induced ferromagnetism of Ga 1 x Mn x As is ascribed to a double-exchange (DE)-like mechanism realized in the magnetic impurity band/or in the band tail.

1. Introduction

For more than three decades, magnetic semiconductors have attracted much attention because of their combination of magnetic and semiconducting properties.
In the 1960s and 1970s, many papers were published on ferromagnetic semiconductors (FMSs) such as EuO and EuS [1,2]. Pure EuO and EuS are considered as typical Heisenberg ferromagnets. Their magnetic moments originate from the half-filled and highly localized 4 f shell of the Eu ion at regular sites. At the high-temperature limit ( T = ) , the orientation of each localized spin (f spin) is completely random. At paramagnetic temperatures the total magnetization is zero, although short-range order is formed as the temperature decreases to the Curie temperature ( T c ) . Below the Curie temperature, spontaneous magnetic ordering develops, and at T = 0 the orientations of f spins are perfectly arranged in one direction. Eu chalcogenides are also insulators with a NaCl-type structure. When a single electron is introduced into the crystal, the electron enters an otherwise empty ( 5 d , 6 s ) conduction band. The electron (hereafter referred to as the s electron) moves in the crystal while interacting with f spins through the s-f exchange interaction. Thus, the conduction- (s-) electron state in FMSs is strongly affected by the magnetic order of f spins. This causes numerous anomalous phenomena in FMSs, such as the redshift of the optical absorption edge with decreasing temperature, the magnetic polaron effect, spin-disorder scattering around T c [3], and the metal-insulator transition in Eu-rich EuO [4]. Its extreme properties, such as high electron-spin polarization, colossal magnetic resistance (CMR), and the enhancement of T c due to Gd doping, make electron-doped EuO interesting for spintronics applications. Recently, these features have therefore stimulated systematic experimental studies with modern techniques and improved sample quality [5,6,7,8,9,10,11,12,13] as well as theoretical studies [14,15,16,17,18]. The conduction-electron states in FMSs have been discussed on the basis of the s-f (exchange) model [19,20]; the s-f model is sometimes referred to as the Kondo lattice model [21,22].
Diluted magnetic semiconductors (DMSs) are semiconducting alloys whose lattice is partly made up of substitutional magnetic atoms. The most extensively studied DMSs since the 1980s are A 1 x II Mn x B VI -type (II-VI-based) DMSs, in which a fraction of the group II sublattice is replaced at random by magnetic Mn atoms. Mn impurities substituting for 2+ cations act as stable 2+ ions; therefore, there are few carriers, which makes these DMSs insulators. It is widely accepted that in II-VI-based DMSs the carriers (s electrons and p holes) move over many sites while interacting with the localized (d) spins at Mn sites through the sp-d exchange interaction [23,24,25]. The exchange interaction strongly enhances the effect of the magnetic field on band splitting, leading to spectacular magnetooptical effects (e.g., giant Faraday rotation or Zeeman splitting). Since 1996, attention has also been focused on the III-V-based DMSs of Ga 1 x Mn x As and In 1 x Mn x As owing to their high potential for new device applications. It is highly noteworthy that the doping of Mn into GaAs and InAs leads to ferromagnetism and interesting magneto-optical and magnetotransport phenomena. This ferromagnetism is generally referred to as “carrier-induced ferromagnetism” because hole carriers introduced by Mn incorporation mediate the ferromagnetic coupling between Mn ions [26,27]. The main difference between II-VI and III-V-based DMSs is that in the latter, Mn ions act not only as magnetic impurities but also as acceptors. The microscopic mechanism for carrier-induced ferromagnetism is still controversial. To elucidate the origin and mechanism of the carrier-induced ferromagnetism, however, it appears necessary to clarify the nature of the carrier states in DMSs.
The carrier states in FMSs and DMSs have not yet been explained theoretically in a sufficiently clear way. The exchange interaction between the carrier and magnetic moments (localized spin), however, seems to be a clue to solving most of these problems. Since the localized spins thermally fluctuate in FMSs, the theory should properly take into account the effect of the thermally fluctuating spins on carrier states through the exchange interaction. In DMSs, substitutional disorder exists in addition to the thermal fluctuation of localized spins at Mn sites.
The coherent potential approximation (CPA) is a superior mean-field theory which was originally devised for describing the electronic structure and/or the properties of binary substitutional alloys [28,29,30]. Rangette e t a l . first applied the CPA to the s-f model at high temperatures, assuming the orientations of the magnetic moments to be completely random [31]. In 1974 Kubo extended the CPA to thermally fluctuating spin systems; Kubo formulated the CPA in terms of effective locators [32]. In 1996 Takahashi e t a l . formulated the CPA in the t-matrix formula, which has been proved to be equivalent to the locator-formula CPA. Since then, the CPA has been applied to investigate FMSs [33,34,35,36,37] and/or DMSs [38,39,40,41,42,43]. This method is referred to as “dynamical CPA", because the dynamical spin-flip process is properly taken into account in the thermal averaging operation. In the classical spin limit, the numerical results obtained by dynamical CPA are in good agreement with those obtained by dynamical mean-field theory (DMFT) [44,45]. It has also been reported that the result for optical conductivity obtained by the dynamical CPA is in reasonable agreement with that obtained by Monte Carlo (MC) simulation [46]. Note that the dynamical CPA can be applied even for a finite magnitude of localized spin. Furthermore, the analytic formula in limiting cases, such as at high temperatures and at the diluted impurity limit, can be easily deduced.
In this article we briefly review the study of magnetic semiconductors using the dynamical CPA. We first discuss the conduction-electron states in FMSs on the basis of the s-f model. We formulate the dynamical CPA using multiple-scattering theory and discuss the numerical results. Next, we study the carrier states in II-VI and III-V-based DMSs using a simple model. The numerical results and discussion are first given for two typical cases with no nonmagnetic attractive potential: (i) the case of a strong exchange interaction and (ii) the case of a moderate exchange interaction. Since the screened Coulomb attractive potential acts between a hole and a Mn acceptor center in III-V-based DMSs, we also study (iii) the case of a moderate exchange interaction with a strong nonmagnetic attractive potential. Then, based on the Curie temperature T c numerically estimated in a simple way, we investigate the type and properties of the magnetism that may occur when carriers are introduced into DMSs. After that, we focus on (Ga,Mn)As, which has attracted much attention in recent years owing to its so-called carrier-induced ferromagnetism. Throughout the present review, we study the effect of the exchange interaction between the carrier and the localized spins on carrier states in magnetic semiconductors. In the Appendix, the locator formula of the dynamical CPA is briefly summarized.

2. Conduction-electron States in a Ferromagnetic Semiconductor (FMS)

2.1. Coherent Potential Approach to the s-f Model

The s-f (exchange) model is currently accepted as a basis for studying the conduction-electron states in an ordinary FMS such as EuO or EuS. In this model, there are magnetic moments at the regular sites and a well-defined conduction band. A single conduction electron (s electron) moves in the crystal while interacting with localized magnetic moments (f spins) through the s-f exchange interaction. Therefore, the total Hamiltonian H t is expressed as
H t = H s + H f + H s f
where
H s = k , μ ε k a k μ a k μ
H f = m n J m n S m · S n
H s f = I n μ , ν a n μ σ μ ν · S n a n ν
Here, H s represents the translational energy of an s electron; a k μ and a k μ are, respectively, the creation and annihilation operators for the Bloch state of an s electron with spin μ and wave vector k, and ε k is the energy of the Bloch state. H f is the Heisenberg-type exchange interaction between f spins; S n is the operator of the f spin located at the nth lattice site, and J m n is the exchange interaction constant between f spins at the mth and nth sites. The s-f exchange interaction between an s electron and f spins, H s f , is represented by the simplest form of the intra-atomic exchange model; I is the exchange constant and σ μ ν is an element of the Pauli matrix for an s electron; a n μ and a n μ are the creation and annihilation operators for the Wannier states of an s electron with spin μ at site n, respectively.
Generally, magnetic excitation energy is very small compared with the conduction bandwidth and the s-f exchange energy. Thus, we treat the f spins as a quasi-static system, that is, we take the thermal average for fluctuating f spins at the final stage of the derivation of physical quantities. Throughout this paper, we shall confine our discussion to the so-called one-particle picture. Thus, we define the single-particle Green’s function as
G ( ω ) = 1 ω H
with
H = H s + H s f
and write its thermal average for f spins as G ( ω ) . Hereafter, the variable ω will be omitted from the operators for cases where this will cause no confusion.
To apply the multiple-scattering theory [28,29], we divide H into the unperturbed Hamiltonian K and the perturbation term V. When magnetization arises, an s electron in an FMS is subjected to different effective potentials through the s-f exchange interaction depending on the orientation of its spin. Thus, we assume a spin-dependent effective medium in which a carrier is subject to a coherent potential, Σ or Σ , according to the orientation of its spin. The coherent potential Σ ( Σ ) is an energy (ω)-dependent complex potential. Then, an s electron moving in this effective medium is described by the (unperturbed) reference Hamiltonian K:
K = k μ ( ε k + Σ μ ) a k μ a k μ
Thus, the perturbation term V ( H K ) is written as the following sum over each lattice site:
V = n v n
with
v n = μ , ν a n μ ( I σ μ ν · S n Σ μ δ μ ν ) a n μ
Next, using the reference Green’s function P given by
P ( ω ) = 1 ω K
we define the t-matrix of the s-f exchange interaction due to a localized spin S n embedded in the effective medium by
t n = v n [ 1 P v n ] 1
Note that t n represents the complete scattering associated with the isolated potential v n in the effective medium, and that K, and thus P, includes no localized spin operator. According to the multiple-scattering theory, the total scattering operator T, which is related to G 1 / ( ω H ) as
G = P + P T P
is expressed as the multiple-scattering series
T = n t n + n t n P m ( n ) t m + n t n P m ( n ) t m P l ( m ) t l +
Within the single-site approximation, the condition
t n = 0 , for   any   site   n
leads to T 0 , and thus G P . This is the CPA.
Here, we introduce the diagonal component of P in the Wannier representation, F μ ( ω ) n μ | P ( ω ) | n μ (independent of n), as
F μ ( ω ) n μ | P ( ω ) | n μ = 1 N k 1 ω ε k Σ μ
In the CPA, the coherent potential Σ μ is set to satisfy the condition Equation (14) so that G P . Thus, the density of states (DOS) of the electron with spin μ ( = , or ) is calculated as
D μ ( ω ) = 1 N Tr μ | δ ( ω H ) | μ = 1 N π Im Tr μ | 1 ω + i 0 H | μ
1 N π Im Tr μ | 1 ω K | μ = 1 N π Im k 1 ω ε k Σ μ = 1 π Im F μ ( ω )

2.2. t-Matrix Elements and Their Thermal Average

Here we show that the t-matrix elements of the s-f exchange interaction for the f spin embedded in the effective medium, Σ ( Σ ), can be derived without further approximation [33,43]. Because the exchange term v has four spin-matrix elements, the t-matrix also has four components. For simplicity, we will omit the site suffix n in the Wannier representation; t μ ν n μ | t | n ν , where μ , ν = and/or ↓. In accordance with the definition of the t-matrix in Equation (11), we have
t [ 1 P v ] = v
Equation (17) is written in the spin-matrix-element expression as
t t F v t F v = v
t t F v t F v = v
where F μ μ | P | μ . Then, Equation (18a) × ( F v ) 1 + Equation (18b) × ( 1 F v ) 1 leads to an equation including t only ( t is canceled):
t [ ( 1 F v ) ( F v ) 1 F v ( 1 F v ) 1 ] = v ( F v ) 1 + v ( 1 F v ) 1
Recalling the commutation relationships between the components of S,
S S z = ( S z + 1 ) S
( S + ) 1 = [ S ( S + 1 ) ( S z ) 2 S z ] 1 ( S )
we obtain an explicit expression for t using no further approximations. Here, S is the magnitude of the localized spin S , and the z component of the f spin is S z : S + = S x + i S y and S = S x i S y . Other t-matrix elements are also obtained by a similar procedure.
In order to show the resulting expressions simply, it is convenient to introduce the following symbols:
V v = I S z Σ
V v = + I S z Σ
v = I S
v = I S +
U I ( S z 1 ) Σ
U + I ( S z + 1 ) Σ
W I 2 S S + = I 2 [ S ( S + 1 ) S z 2 S z ]
W I 2 S + S = I 2 [ S ( S + 1 ) S z 2 + S z ]
The physical meanings of the above symbols can be easily explained. V ( V ) is the spin-diagonal component of the s-f exchange interaction, wherein an s electron with ( ) spin interacts with an f spin embedded in the medium of Σ ( Σ ). U ( U ) is the spin-diagonal component of the s-f exchange interaction, wherein an s electron with ( ) spin interacts with an f spin that has already flipped in the previous scattering; thus, the f spin operator S z is replaced by S z 1 ( S z + 1 ). Both V ( V ) and U ( U ) describe the scattering process without a spin flip. On the other hand, W ( W ) represents the interaction energy required by an s electron with ↑ (↓) spin to flip and then reverse its spin after intermediate propagation with a flipped spin. The resulting expressions are
t = V + F ( W V U ) 1 F U F [ V + F ( W V U ) ]
t = V + F ( W V U ) 1 F U F [ V + F ( W V U ) ]
t = 1 1 F U F [ V + F ( W V U ) ] ( I S ) = ( I S ) 1 1 F U F [ V + F ( W V U ) ]
t = 1 1 F U F [ V + F ( W V U ) ] ( I S + ) = ( I S + ) 1 1 F U F [ V + F ( W V U ) ]
Note that the diagonal-matrix element t μ μ involves only S z as an operator [i.e., t μ μ t μ μ ( S z ) ]. Thus, the thermal average over the fluctuation of the localized spin is taken as
t μ μ = S z = S S t μ μ ( S z ) exp h S z k B T / S z = S S exp h S z k B T
where h denotes the effective field acting on the localized f spins. Since there is a one-to-one correspondence between S z and the parameter λ h / k B T through the relationship
S z = S z = S S S z exp h S z k B T / S z = S S exp h S z k B T
we can describe the electron states in terms of S z instead of λ. Note that the thermal average of the off-diagonal elements t = t is 0 because the magnetization is assumed to be along the z axis.
The conditions for the CPA are expressed as
t = 0
t = 0
It is worth noting that the spin-flip processes are properly taken into account in the above expression for the t-matrix elements. As a result, each t matrix element t μ μ depends on both Σ and Σ Equations (26a) and (26b) simultaneously.
For an undisturbed DOS, we assume the model DOS to have a semicircular form with a half-bandwidth Δ:
ρ ( ε ) = 2 π Δ 1 ε Δ 2
Then, F μ ( ω ) is calculated as
F μ ( ω ) = 1 N k 1 ω ε k Σ μ = Δ Δ ρ ( ε ) ω ε Σ μ d ε
= 2 Δ ω Σ μ Δ ω Σ μ Δ 2 1
The procedure of the numerical calculation is as follows. For a given ω, by assigning a suitable complex value of Σ μ (for μ = or ↓), F μ is calculated by Equation (29). Taking the thermal average for fluctuating f spins using Equation (24), Σ and Σ are simultaneously determined by Equation (26). Then, we can calculate F μ again. This procedure is repeated until the calculation converges (see Reference [33] for details). In all of the present numerical results, we have numerically verified that
D ( ω ) d ω = D ( ω ) d ω = 1

2.3. Results for the Conduction-electron States in an FMS

The parameters that are necessary to describe the present model are the conduction bandwidth 2 Δ , the exchange energy I S ( = I × S ) , the magnetization S z / S , and the quantum spin factor 1 / S . Here, we take S = 7 / 2 for the f spin. We first confirm that the exchange interaction term I σ · S has two eigenstates (i.e., the parallel-coupling state and antiparallel-coupling state) according to the manner of coupling between the s electron’s spin and the localized f spin. The parallel-coupling state (denoted by p) has an energy eigenvalue of ε p = I S with the degeneracy of 2 S + 2 , while the antiparallel-coupling state (denoted by a) has an energy eigenvalue of ε a = + I ( S + 1 ) with the degeneracy of 2 S . Therefore, the band splits into two subbands when the exchange energy I S is large compared with the bandwidth 2 Δ .
In Figure 1(a), we show the DOS for various exchange energies I S / Δ in the paramagnetic states. At high temperatures the orientation of f spins is completely random. Therefore, the present result is equivalent to that obtained by Rangette e t a l . [31], who applied the CPA to an electron propagating in a disordered binary alloy in which two kinds of atoms with atomic energies of I S and + I ( S + 1 ) are randomly distributed with concentrations of ( S + 1 ) / ( 2 S + 1 ) and S / ( 2 S + 1 ) , respectively. When I S / Δ is small, the band remains as a single band, while it broadens with the increase of I S / Δ . For I ( 2 S + 1 ) Δ (or I S / Δ 0 . 44 ) the band splits into two subbands, the parallel-coupling and antiparallel-coupling subbands, whose centers are at I S and + I ( S + 1 ) , and whose total numbers of states are 2 ( S + 1 ) / ( 2 S + 1 ) = 1 . 125 and 2 S / ( 2 S + 1 ) = 0 . 875 , respectively.
In Figure 1(b), we show the DOS in the completely ferromagnetic states. In this case, only the value S z = S is realized in the thermal average over the f spin states. Hence, from Equation (26) we obtain
Σ = I S
Σ = + I S ( 1 + I F ) ( 1 I F )
Figure 1. The DOS shown as a function of ω/Δ for I S /Δ = 0.0, 0.2, 0.4, 0.6, 0.8, and 1.0: (a) paramagnetic states ( S z = 0; left panel) and (b) completely ferromagnetic states ( S z = S; right panel). The solid curves show up-spin states and the dotted curves show down-spin states. The points on the energy axes indicate I S and + I ( S + 1 ) , and the straight dotted lines indicate ( Δ + I S ) and + ( Δ + I S ) for reference.
Figure 1. The DOS shown as a function of ω/Δ for I S /Δ = 0.0, 0.2, 0.4, 0.6, 0.8, and 1.0: (a) paramagnetic states ( S z = 0; left panel) and (b) completely ferromagnetic states ( S z = S; right panel). The solid curves show up-spin states and the dotted curves show down-spin states. The points on the energy axes indicate I S and + I ( S + 1 ) , and the straight dotted lines indicate ( Δ + I S ) and + ( Δ + I S ) for reference.
Materials 03 03740 g001
The results can be interpreted as follows. For the completely ferromagnetic case (i.e., T = 0 ), the states of an s electron with up spin only shift by I S with no damping. On the other hand, the s-electron states with down spin are damped because they can flip their spin under the condition that the total spin ( = S 1 / 2 ) is conserved if the DOS with up spin is not zero therein. This is because Σ involves F [ F ( ω ) ] . This spin-flip process of the s electron is a quantum effect due to the finiteness of the magnitude of the f spin. Thus, in the classical spin limit [i.e., setting S 1 and S z 1 in Equation (23)] , Equation (31b) is replaced by
Σ = I S
In Figure 2, we show the DOS for two typical cases: (a) weak exchange interaction (left panel; I S / Δ = 0 . 2 ) and (b) strong exchange interaction (right panel; I S / Δ = 0 . 8 ). In the weak exchange interaction limit ( I S / Δ 0 ), the present results agree with those obtained using first-order perturbation theory. Substituting Equation (23a) into Equation (26a) and Equation (23b) into Equation (26b), and expanding them in I to the first order, we obtain
V = 0 then Σ = I S z
V = 0 then Σ = + I S z
This means that the ferromagnetic ordering of f spins gives rise to the I S z shift in the up-spin band and the + I S z shift in the down-spin band. However, even when I S / Δ = 0 . 2 , this is not the case, as is shown in Figure 2(a). The bands are broadened owing to the fluctuation of the f spins, and the down-spin band has a tail which reaches the bottom of the up-spin band even in the case of completely ferromagnetic states ( S z = S ). This explains why the electron-spin polarization cannot reach 100%; the origin is the quantum effect coming from the finiteness of the magnitude of the f spin, as already discussed [34].
Figure 2. The DOS shown for magnetizations S z / S = 0.0, 0.2, 0.4, 0.6, 0.8, and 1.0: (a) I S / Δ = 0 . 2 (left panel) and (b) I S / Δ = 0 . 8 (right panel). The solid curves show up-spin states and the dotted curves show down-spin states.
Figure 2. The DOS shown for magnetizations S z / S = 0.0, 0.2, 0.4, 0.6, 0.8, and 1.0: (a) I S / Δ = 0 . 2 (left panel) and (b) I S / Δ = 0 . 8 (right panel). The solid curves show up-spin states and the dotted curves show down-spin states.
Materials 03 03740 g002
As shown in Figure 2(b), in a case of the strong exchange interaction, the band splits into two subbands depending on the coupling of the s electron spin and f spins: the parallel-coupling subband (lower-energy side) and the antiparallel-coupling subband (higher-energy side). The total number of states in the parallel-coupling subband per site is 2 ( S + 1 ) / ( 2 S + 1 ) = 1 . 125 and that in the antiparallel-coupling subband is 2 S / ( 2 S + 1 ) = 0 . 875 , irrespective of the value of S z . When S z = 0 , both subbands are composed of the same number of up- and down-spin states. When S z = S , on the contrary, all states in the antiparallel-coupling subband have down spin, while the states in the parallel-coupling subband are composed of all the up-spin states and part of the down-spin states; the number of states with down spin per site is 0 . 125 in the parallel-coupling subband.
The energy of the bottom of the band ω b is shifted by the s-f exchange interaction from that of the undisturbed (model) band ( ω = Δ ), even when S z = 0 . In Figure 3(a), we show the energy shift of the bottom of the band in the paramagnetic states normalized by I S , ( ω b + Δ ) / I S , as a function of I S / Δ ; the exact solution for ω b is given in the Appendix of Reference [38]. Using the result in Figure 3(a), we can explain why Eu chalcogenide FMSs exhibit different redshifts despite the fact that they all have the same lattice structure, the same f spin value of S = 7 / 2 , and almost the same exchange interaction energy I S [47]. Since the optical absorption band is assigned to the 4 f 7 4 f 6 5 d ( t 2 g ) band transition, the redshift is ascribed to the lowering of the energy of the bottom of the (d-like) conduction band with the decrease in temperature due to the (d-f) exchange interaction [1]. In Figure 3(a), we indicate the magnetic redshifts experimentally observed for EuO (0.27 eV), EuS (0.18 eV), and EuSe (0.13 eV) with arrows. In this approach, ω b is a function of I S / Δ at paramagnetic temperatures, while ω b = I S (independent of Δ) at T = 0 . Since the exchange interaction has an intra-atomic character due to the strong localization of f electrons within the Eu 2 + ion, the value of I S does not greatly differ among these Eu chalcogenides. Thus, the difference in the total redshift can be ascribed to the difference in the bandwidth 2 Δ . From Figure 3(a), using I S = 0 . 35 eV we estimated the values of the bandwidth ( 2 Δ ) as 7 eV (EuO), 2.5 eV (EuS), and 1.6 eV (EuSe). Though uncertainty remains due to the experimental error in measuring the position of the absorption edge, the broad conduction-band picture for EuO is consistent with the result recently obtained by spin-resolved spectroscopy [12]. The tendency that the width of the conduction band of EuX decreases with the change of chalcogenides from X = O to Se is consistent with the electronic structure obtained by optical measurement [1].
Figure 3. (a) The shift of the bottom of the band in paramagnetic states normalized by the exchange energy, ( ω b + Δ ) / I S , shown as a function of I S / Δ . The arrows indicate the experimentally obtained magnetic redshifts of EuO (0.27 eV), EuS (0.18 eV), and EuSe (0.13 eV) with I S = 0 . 35 eV. (b) The energy of the bottom of the band depicted as a function of T / T c for I S / Δ = 0 . 2 , 0 . 4 , 0 . 6 , and 1 . 0 . The value of S z / S , which corresponds to the result of the weak exchange interaction limit, is also shown.
Figure 3. (a) The shift of the bottom of the band in paramagnetic states normalized by the exchange energy, ( ω b + Δ ) / I S , shown as a function of I S / Δ . The arrows indicate the experimentally obtained magnetic redshifts of EuO (0.27 eV), EuS (0.18 eV), and EuSe (0.13 eV) with I S = 0 . 35 eV. (b) The energy of the bottom of the band depicted as a function of T / T c for I S / Δ = 0 . 2 , 0 . 4 , 0 . 6 , and 1 . 0 . The value of S z / S , which corresponds to the result of the weak exchange interaction limit, is also shown.
Materials 03 03740 g003
In Figure 3(b), the energy shift of the bottom of the band is shown as a function of normalized temperature T / T c . From Figure 3(b), we can also explain the reason for the apparent success of first-order perturbation theory in explaining the magnetic redshift. The temperature dependence of the shift in the energy of the bottom of the band is approximately I eff S z . Thus, I eff can be regarded to be the effective exchange constant in which the effect of multiple scattering has already been renormalized.
The results obtained by applying the dynamical CPA for the s-f model of FMSs are consistent with the experimental observation of the electron-spin polarization [34], and they explain the temperature dependence of the electrical resistivity of degenerate FMSs [35]. We should also add that the coherent potential approach can explain not only the redshift of FMSs but also the blue shift of the antiferromagnetic semiconductor EuTe [36,37].

3. Carrier States in Diluted Magnetic Semiconductors (DMSs)

3.1. Model Hamiltonian for a Carrier in a DMS and the Application of the Dynamical CPA

In order to study the effect of the sp-d exchange interaction between a carrier (an s conduction electron or p hole) and localized magnetic moments (d spins) together with magnetic and chemical disorder in DMSs, we introduce a simple model for A 1 x Mn x B -type DMSs. In this model, the local potentials of nonmagnetic (A) ions in a semiconducting compound ( A B ) are substituted randomly, with mole fraction ( x ) , by local potentials that include the exchange interaction between a carrier and the localized spin moment on a Mn (denoted by M) ion. Thus, the potential to which a carrier is subjected at a site depends on whether the site is occupied by an A ion or M ion. The Hamiltonian H is given by
H = m , n , μ ε m n a m μ a n μ + n u n
where u n is either u n A (at the A site) or u n M (at the M site), depending on the ion species occupying the n site:
u n A = E A μ a n μ a n μ
u n M = E M μ a n μ a n μ I μ , ν a n μ σ μ ν · S n a n ν
Here, a n μ and a n μ are, respectively, the creation and annihilation operators for a carrier with spin μ at the n site. The transfer-matrix element between m and n, ε m n , is assumed to be independent of the types of constituent atoms that occupy the m and n sites. In II-VI-based DMSs of the A 1 x II Mn x B VI -type, E A ( E M ) represents a nonmagnetic local potential at the A 2 + (Mn 2 + ) sites. In III-V-based DMSs such as Ga 1 x Mn x As, the spin-independent potential E M ( < 0 ) can be regarded as a screened Coulomb attractive potential between a carrier (hole) and the Mn 2 + ion (acceptor center). The exchange interaction between the carrier and localized spin S n of the Mn site n is expressed by I a n μ σ μ ν · S n a n ν , where σ μ ν represents the element of the Pauli spin matrices. We disregard the electron-electron, hole-hole, and/or electron-hole interactions.
The application of the dynamical CPA to the Hamiltonian in Equation (34) is straightforward [41,42,43]. Since the present system includes both substitutional disorder and the thermal fluctuation of the localized spin at an M site, the average of the t-matrix is written as
t av = ( 1 x ) t A + x t M
Here, we express the average of t n over the disorder in the system as t n av ; ( 1 x ) and x are the mole fractions of A and M atoms, respectively. t A is the t-matrix that represents the multiple scattering of a carrier due to the A ion potential u A [in Equation (35)] embedded in the effective medium. Assuming the spin-dependent coherent potential Σ μ ( μ = or ↓), the t-matrix elements are given as [29]
t μ μ A = E A Σ μ 1 ( E A Σ μ ) F μ
Here, F μ is a diagonal matrix element of a propagator with respect to the effective medium and is given by Equations (15) and/or (29). t M is the t-matrix that represents the multiple scattering of a carrier due to the M ion potential u M [in Equation (36)] embedded in the effective medium; t M is the thermal average of t M over the fluctuation of the localized spin. Explicit expressions for t μ ν are obtained after minor substitutions in Equation (23); Σ μ Σ μ E M . In the dynamical CPA, the coherent potential Σ μ is set such that the effective scattering of a carrier at the chosen site embedded in the effective medium is zero on average. Thus, the dynamical CPA conditions are given by
( 1 x ) t A + x t M = 0
( 1 x ) t A + x t M = 0
In Appendix A, we outline the dynamical CPA using locator formalism, which has been proved to be equivalent with the dynamical CPA using t-matrix formalism. The advantage of the locator formalism is that it can be easily used to obtain the species-resolved DOSs D μ A ( ω ) and D μ M ( ω ) , which are the DOSs with μ spin associated with an A ion and M ion, respectively.

3.2. General Consideration for the Carrier States in a DMS

In the subsections below, we treat the localized spins classically, although S = 5 / 2 for the Mn 2 + ion. One of the advantages of the classical spin approximation is that it reduces the number of physical parameters. It is sufficient to assign a value to the exchange energy I S = I × S instead of assigning the values of I and S separately. In the classical spin approximation, we let 1 / S approach 0 (i.e., S ) while keeping I S constant. The quantum fluctuation of the localized spin is ignored. As a consequence, both eigenstates have the degeneracy of 2 S , and the eigenenergies become symmetric: ε p = I S and ε a = + I S . The s-d exchange interaction between a conduction electron (s electron) and a localized spin (d spin) favors parallel coupling. On the other hand, in most DMSs, the p-d exchange interaction between a hole (p hole) and a localized spin favors antiparallel coupling, and the magnitude of the p-d exchange interaction is several times larger than that of the s-d exchange interaction. The p-d exchange interaction plays an important role in magneto-optical effects in II-VI-based DMSs and is related to the carrier-induced ferromagnetism in III-V-based DMSs. Hence, keeping the p-d exchange interaction in mind, we assume I S < 0 hereafter. Note that the present model requires only two parameters, I S / Δ and E M / Δ , after we set E A 0 .
In Section 3.3, Section 3.4 and Section 3.5, we will discuss some typical cases of A 1 x M x B -type DMSs in which 5% of the nonmagnetic (A) ions are randomly substituted by magnetic (M) ions. The results reveal the nature of the magnetic impurity bands and how the carrier states behave with changing magnetization. In Section 3.6, based on the Curie temperature T c calculated in a simple way, we will discuss the type and properties of ferromagnetism which may occur when carriers introduced into A 1 x M x B -type DMSs. In Section 3.7, we will focus on the case of (Ga, Mn)As.

3.3. The Case of Strong Exchange Interaction

In Figure 4, we show the numerical results with I S = Δ and E M = 0 . 0 , referred to as the case of strong exchange interaction hereafter. In the left panel of Figure 4, the spin-polarized DOSs, D ( ω ) and D ( ω ) , are depicted for various values of S z / S . In the dilute impurity limit ( x 0 ), impurity levels appear at the energies of E a Δ = E M I S Δ + 1 4 Δ E M I S = ± 1.25. When x = 0.05, impurity bands appear around the impurity levels. The total number of states of each impurity band is x, irrespective of S z . The lower (higher)-energy impurity band corresponds to the antiparallel (parallel)-coupling state.
Figure 4. Left panel: DOS as a function of ω / Δ for various values of S z / S . The solid lines represent the down-spin carrier and dotted lines represent the up-spin carrier. The impurity levels E a = ± 1 . 25 Δ are indicated by dots on the abscissa. Right panel: (a) A-site component of the DOS, ( 1 x ) D A ( ω ) Δ and ( 1 x ) D A ( ω ) Δ . (b) M-site component of the DOS, x D M ( ω ) Δ and x D M ( ω ) Δ . The thick, thin, and dotted lines represent the cases of S z / S = 1 . 0 , 0 . 5 , and 0 . 0 , respectively. Note the difference in the scale of the vertical axes of (a) and (b). From Takahashi and Kubo [42].
Figure 4. Left panel: DOS as a function of ω / Δ for various values of S z / S . The solid lines represent the down-spin carrier and dotted lines represent the up-spin carrier. The impurity levels E a = ± 1 . 25 Δ are indicated by dots on the abscissa. Right panel: (a) A-site component of the DOS, ( 1 x ) D A ( ω ) Δ and ( 1 x ) D A ( ω ) Δ . (b) M-site component of the DOS, x D M ( ω ) Δ and x D M ( ω ) Δ . The thick, thin, and dotted lines represent the cases of S z / S = 1 . 0 , 0 . 5 , and 0 . 0 , respectively. Note the difference in the scale of the vertical axes of (a) and (b). From Takahashi and Kubo [42].
Materials 03 03740 g004
The impurity bands are strongly affected by changes in S z . On the other hand, the host band is negligibly affected. To elucidate the origin of the change in the DOS, we calculate the species-resolved DOS. In the left panel of Figure 4, we depict the A- and M-site components of the DOS, ( 1 x ) D μ A ( ω ) and x D μ M ( ω ) , respectively. D μ A ( ω ) [ D μ M ( ω ) ] represents the local DOS with μ spin ( μ = or ↓) associated with the A (M) ion. Note that
D μ ( ω ) = ( 1 x ) D μ A ( ω ) + x D μ M ( ω )
Since D A ( ω ) and D M ( ω ) are normalized, the total numbers of A-site states and M-site states are 1 x and x, respectively. The numerical result shown in the right panel reveals that the impurity state is mainly composed of M-site states and that the change in the impurity band is mainly ascribed to the change in D μ M ( ω ) .

3.4. The Case of Moderate Exchange Strength (II-VI-based DMSs)

To the best of our knowledge, no magnetic impurity band has been observed in A 1 x II Mn x B VI -type DMSs [23]. Thus, II-VI-based DMSs correspond to the present model with | E M ± I S | < 0 . 5 Δ . For convenience, we take the exchange energy I S = 0 . 4 Δ and the band offset energy E M = 0 . 0 for II-VI-based DMSs. Numerical results for the parameters are presented in Figure 5, Figure 6 and Figure 7. In the left panel of Figure 5, the spin-polarized DOS is shown for various values of S z . As can be seen in Figure 5(b), the carrier states at the M site spread over the whole range of band energy, in contrast with the case of strong exchange interaction. The weak S z dependence of x D M ( ω ) suggests that the coupling between the carrier spin and the localized spin is not strong, except at the band edges.
Figure 5. The results for II-VI-based DMSs. Left panel: DOS as a function of ω / Δ for various values of S z / S . The solid lines represent the down-spin carrier and the dotted lines represent the up-spin carrier. The arrow indicates the Fermi level ε F / Δ for n = x ( = 0 . 05 ) . Right panel: (a) A-site component of the DOS, ( 1 x ) D A ( ω ) Δ and ( 1 x ) D A ( ω ) Δ . (b) M-site component of the DOS, x D M ( ω ) Δ and x D M ( ω ) Δ . The thick, thin, and dotted lines represent the cases of S z / S = 1 . 0 , 0 . 5 , and 0 . 0 , respectively. From Takahashi and Kubo [42].
Figure 5. The results for II-VI-based DMSs. Left panel: DOS as a function of ω / Δ for various values of S z / S . The solid lines represent the down-spin carrier and the dotted lines represent the up-spin carrier. The arrow indicates the Fermi level ε F / Δ for n = x ( = 0 . 05 ) . Right panel: (a) A-site component of the DOS, ( 1 x ) D A ( ω ) Δ and ( 1 x ) D A ( ω ) Δ . (b) M-site component of the DOS, x D M ( ω ) Δ and x D M ( ω ) Δ . The thick, thin, and dotted lines represent the cases of S z / S = 1 . 0 , 0 . 5 , and 0 . 0 , respectively. From Takahashi and Kubo [42].
Materials 03 03740 g005
The most powerful tool for studying the exchange interaction between a carrier and localized spins in II-VI-based DMSs is optical measurement. In most II-VI-based DMSs the energy splitting between σ + and σ A excitons in a DMS is well described as [23,24]
δ E = N 0 ( α β ) x S z
The virtual crystal approximation (VCA) is a first-order perturbation theory with respect to the exchange interaction. When applying the VCA to the present model, we obtain the energy shift in the band edge due to the M-site local potential:
Σ b ( ± ) = ( E M I S z ) x
depending on the carrier’s spin ( ± ) . Therefore, N 0 α and N 0 β are regarded as the exchange constants for conduction electrons and valence electrons, respectively; N 0 α = 2 I for an electron and N 0 β = 2 I for a hole. However, some experimental results indicate that the application of the VCA is limited [38], although the VCA has been widely accepted as describing the field-induced exchange splitting of extended states in II-VI-based DMSs [23,24,25]. Therefore, it is important to clarify the limit of application of the VCA and to devise a theoretical treatment beyond the VCA.
Here we consider the optical absorption spectrum on the basis of the dynamical CPA treatment in which the multiple-scattering effect is considered [43]. In calculating the optical absorption spectrum, we assume that the transition dipole moments of the A and M ions are the same. Under this assumption, the optical absorption spectrum is given by the k = 0 components of the DOS. Since the explicit k dependence of ε k is not used in the present framework, we assume that k = 0 corresponds to the minimum point of the model band. Therefore, taking ε 0 = Δ , we define the optical absorption spectrum as [30]
A μ ( ω ) = 1 π Im 1 ω + Δ Σ μ ( ω )
Figure 6(a) shows how the band tail is modified and spin-polarized with the development of magnetization. Note that even with S z = 0 , the band is not the same as the model band, owing to the disorder of the random distribution of M ions and the fluctuation of localized spins. With increasing S z , the bottom of the down-spin band shifts to a lower energy, accompanied by an energy shift of the bottom of the up-spin band. The two band edges agree with each other except in the case of S z = S , although the down-spin band is strongly suppressed in the band tail. The agreement of the band edges is a consequence of the spin-flip of a carrier. Thus, the present result for the band-edge shift is very different from that obtained by the VCA. On the other hand, exchange band splitting is observed in magneto-optical measurements such as magneto-absorption and magnetoreflectivity spectra. Hence, we assume that the peak of the optical absorption spectrum A ( ω ) , shown in Figure 6(b), corresponds to the band edge observed in optical measurements. In the right panel of Figure 7, the optical band-edge energy, ω p , at which the optical absorption spectrum shows a peak, is presented for up- and down-spin bands as a function of S z . The behavior of ω p convincingly explains the asymmetrical splitting of the Zeeman energy component; when a magnetic field is applied, the pattern of spin splitting of the A exciton term is asymmetric relative to the position of S z = 0 [40]. In the right panel of Figure 7, the spin-splitting energy ω p ( up ) - ω p ( down ) is displayed as a function of x S z for various values of x. The data for each x are fit by a straight line. With increasing x, the slope of the line, corresponding to N 0 β , decreases. The same behavior has been experimentally observed in Zn 1 x Mn x Te [48] and Cd 1 x Mn x Te [49].
Figure 6. The results for II-VI-based DMSs. (a) Low-energy part of the DOS D ( ω ) . (b) Optical absorption spectrum A ( ω ) in arbitrary units (arb. units). The solid lines represent the down-spin carrier and the dotted lines represent the up-spin carrier. The values of the band-edge energy, ω b / Δ , obtained by a simple approximation are plotted as dots on the lines of S z = S and S z = 0 (see text). Note that the energy of the bottom of the model band is ω = Δ . From Takahashi [43].
Figure 6. The results for II-VI-based DMSs. (a) Low-energy part of the DOS D ( ω ) . (b) Optical absorption spectrum A ( ω ) in arbitrary units (arb. units). The solid lines represent the down-spin carrier and the dotted lines represent the up-spin carrier. The values of the band-edge energy, ω b / Δ , obtained by a simple approximation are plotted as dots on the lines of S z = S and S z = 0 (see text). Note that the energy of the bottom of the model band is ω = Δ . From Takahashi [43].
Materials 03 03740 g006
Assuming that ω b = Δ + Σ ( ω b ) under the dynamical CPA condition [Equation (39)], we can obtain equations for the band-edge energy ω b ; the equations are cubic when S z = 0 and quadratic when S z = S . The approximate values for the band-edge energy are plotted as dots on the lines of S z = S and S z = 0 in Figure 6(a), (b), and in the left panel of Figure 7. Using the approximate values for the band-edge shift with up and down spins, Σ b ( + ) and Σ b ( ) when S z = S , we calculate N 0 β from N 0 β = [ Σ b ( ) Σ b ( + ) ] / x S . The result shows that N 0 β is a function of x, and [39,40]
N 0 β 2 I 1 + 2 E M Δ 2 2 I S Δ 2 when   x 0
N 0 β = 2 I when   x = 1
The present treatment reveals that the apparent enhancement of | N 0 β | with decreasing x observed in Cd 1 x Mn x S is a consequence of the multiple-scattering effect, which is significant for small x in a disordered system [39]. The VCA is applicable when | E M | Δ and | I S | Δ [50].
Figure 7. Left panel: Absorption peak energy ω p / Δ as a function of S z / S . The solid line represents the down-spin carrier and the dotted line represents the up-spin carrier. The error bar represents the half-peak width. The values of ω b / Δ obtained by a simple approximation are plotted as dots on the lines of S z = S and S z = 0 (see text). Right panel: Exchange splitting [ ω p ( up ) ω p ( down ) ] / Δ as a function of x S z / S for various values of x. The straight lines are adjusted to obtain the best fit with each set of x data. From Takahashi [43].
Figure 7. Left panel: Absorption peak energy ω p / Δ as a function of S z / S . The solid line represents the down-spin carrier and the dotted line represents the up-spin carrier. The error bar represents the half-peak width. The values of ω b / Δ obtained by a simple approximation are plotted as dots on the lines of S z = S and S z = 0 (see text). Right panel: Exchange splitting [ ω p ( up ) ω p ( down ) ] / Δ as a function of x S z / S for various values of x. The straight lines are adjusted to obtain the best fit with each set of x data. From Takahashi [43].
Materials 03 03740 g007

3.5. The Case of Strong Attractive Potential

Figure 8 shows the results for I S = 0 . 4 Δ and E M = 0 . 6 Δ . Although the same exchange energy, I S = 0 . 4 Δ , as that for the II-VI-based DMSs, is assumed, the nonmagnetic local potential E M generates an impurity band. Furthermore, comparing the DOSs shown in Figure 4 (strong exchange interaction) and Figure 8 (strong attractive potential), it is easy to see a strong similarity. The lower-energy part of the DOS for I S = 0 . 4 Δ and E M = 0 . 6 Δ is almost the same as that for I S = Δ and E M = 0 . 0 . The reason for the similarity can be explained as follows. First, the impurity level is the same, E a = 1 . 25 Δ , because it is determined by the value of I S + E M ( = Δ ) . When S z = 0 , a magnetic impurity band forms around the impurity level; the number of down- and up-spin states is x / 2 each. The impurity band forms in imitation of the model band. When S z = S , the down-spin bands in the two cases agree with each other because the DOSs were calculated for the same value of E M + I S . When S z = S , the up-spin bands shift toward higher energies and merge into the host band in both cases, although the up-spin bands do not coincide because different values of E M I S were used. Therefore, even in the up-spin bands, we observe a similar tendency in the behavior of the two cases. The above similarities of the carrier state in the two cases can be explained as follows. In DMSs with strong attractive potential, the carrier is so strongly attracted to the M site due to E M ( < 0 ) that the exchange interaction operates effectively in comparison with the case of E M = 0 . Consequently, the exchange interaction in the case of strong attractive potential yields a very similar effect to that in DMSs with a strong exchange interaction.
Figure 8. The results for the DMSs with strong attractive potential. Left panel: DOS as a function of ω / Δ for various values of S z / S . The solid lines represent the down-spin carrier and the dotted lines represent the up-spin carrier. The impurity level E a = 1 . 25 Δ is indicated by the dot on the abscissa. Right panel: (a) A-site component of the DOS, ( 1 x ) D A ( ω ) Δ and ( 1 x ) D A ( ω ) Δ , (b) M-site component of the DOS, x D M ( ω ) Δ and x D M ( ω ) Δ . The thick, thin, and dotted lines represent the cases of S z / S = 1 . 0 , 0 . 5 , and 0 . 0 , respectively. From Takahashi and Kubo [42].
Figure 8. The results for the DMSs with strong attractive potential. Left panel: DOS as a function of ω / Δ for various values of S z / S . The solid lines represent the down-spin carrier and the dotted lines represent the up-spin carrier. The impurity level E a = 1 . 25 Δ is indicated by the dot on the abscissa. Right panel: (a) A-site component of the DOS, ( 1 x ) D A ( ω ) Δ and ( 1 x ) D A ( ω ) Δ , (b) M-site component of the DOS, x D M ( ω ) Δ and x D M ( ω ) Δ . The thick, thin, and dotted lines represent the cases of S z / S = 1 . 0 , 0 . 5 , and 0 . 0 , respectively. From Takahashi and Kubo [42].
Materials 03 03740 g008

3.6. Mechanism of Carrier-induced Ferromagnetism in DMSs

In order to study the mechanism of carrier-induced ferromagnetism that may occur when carriers are introduced into DMSs, we calculate the Curie temperature T c in a very simple way. Throughout this article we assume that the carriers are degenerate. Then we obtain the carrier density with μ spin n μ and the total energy E ( S z ) by
n μ = ε F D μ ( ω ) d ω
and
E ( S z ) = ε F ω [ D ( ω ) + D ( ω ) ] d ω
respectively, as functions of the Fermi level ε F . Note that E ( S z ) is the sum of the kinetic and exchange energies. For a fixed value of S z / S , the total carrier density n ( n + n ) has a one-to-one correspondence with ε F and therefore E ( S z ) can be expressed as a function of n. Thus, we can estimate T c as a function of n using
k B T c = 2 3 x [ E ( 0 ) E ( S ) ]
where E ( 0 ) and E ( S ) are the energies of the paramagnetic state and the completely ferromagnetic state, respectively.
First we investigate the case of E M = 0 . In Figure 9 the result for T c / Δ is presented as a function of n for various values of I S / Δ . We immediately notice that there are two different types of behavior of T c as a function of n depending on the size of | I S | / Δ . When | I S | / Δ is small ( | I S | / Δ 0 . 3 ), ferromagnetism occurs over a wide range of n. The Curie temperature gradually increases with the increase in n and reaches a broad maximum. Then it gently decreases and vanishes at a critical value n c . The maximum T c stays at a low value, and n c is much larger than x ( = 0 . 05 ) . On the other hand, when | I S | / Δ is large ( | I S | / Δ 0 . 7 ), ferromagnetism occurs in a narrow range of n ( x ). T c rises steeply and reaches a maximum at n x x / 2 , and then it decreases rapidly. The maximum T c is high and n c is slightly less than x. These two different features can also be seen clearly in Figure 10, where n c and the maximum T c are depicted as functions of n. The carrier density n x at which T c reaches the maximum is also shown. Two different characteristic features were recognized in the phase diagrams obtained in earlier studies by Chattopadhyay et al. [51], Yagi and Kayanuma [52], and Calderón et al. [53].
Figure 9. The result for Curie temperature T c / Δ as a function of carrier density n for various values of I S / Δ with x = 0 . 05 and E M = 0 . The result based on the assumption that an impurity band has the same shape as the model DOS is drawn as ‘LIMIT’ (see text). From Takahashi and Kubo [42].
Figure 9. The result for Curie temperature T c / Δ as a function of carrier density n for various values of I S / Δ with x = 0 . 05 and E M = 0 . The result based on the assumption that an impurity band has the same shape as the model DOS is drawn as ‘LIMIT’ (see text). From Takahashi and Kubo [42].
Materials 03 03740 g009
When | I S | / Δ 0 . 3 , the maximum T c is approximately proportional to ( I S / Δ ) 2 . This suggests that in the range of | I S | / Δ 0 . 3 the perturbative treatment of I S / Δ is available and that a Ruderman-Kittel-Kasuya-Yosida (RKKY)-like mechanism is expected to operate for a moderate carrier density. In order to clarify the properties of ferromagnetism that occurs when | I S | / Δ is small enough, we show the explicit expression for T c using the mean field approach. Assuming spin-split ± x I S z for the model band, we estimate the gain in total energy E ( 0 ) E ( S z ) to be ρ ( ε F ) ( x I S z ) 2 , where ρ ( ε F ) is the DOS at the Fermi level ε F . Consequently, the T c is obtained as
k B T c = 2 3 x ρ ( ε F ) ( I S ) 2
With a further increase in | I S | / Δ , the maximum T c rises rapidly ( 0 . 3 | I S | / Δ 0 . 7 ) and then tends to saturate. For | I S | / Δ 0 . 7 , ferromagnetism is induced only when n x and the maximum T c occurs at n x x / 2 . The case with I S = Δ shown in Figure 4 belongs to this region.
Figure 10. Phase diagram for E M = 0 and x = 0 . 05 . The critical value n c (solid line; left scale) and the maximum T c (dotted line; right scale) are presented as functions of | I S | / Δ . The carrier density n x at which T c reaches the maximum is included (solid line; left scale). From Takahashi and Kubo [42].
Figure 10. Phase diagram for E M = 0 and x = 0 . 05 . The critical value n c (solid line; left scale) and the maximum T c (dotted line; right scale) are presented as functions of | I S | / Δ . The carrier density n x at which T c reaches the maximum is included (solid line; left scale). From Takahashi and Kubo [42].
Materials 03 03740 g010
Here, we consider the mechanism of ferromagnetism that occurs in the magnetic impurity band of DMSs. In Figure 11, we extract the lower-energy part of the DOS from Figure 4. A magnetic impurity band forms around impurity level “A” and imitates the model band. The total number of states in the impurity band per site is equal to x, irrespective of the value of S z . When S z = S , all states in the impurity band are down-spin states, whereas when S z = 0 , the impurity band is composed of the same number of up- and down-spin states. The impurity band has a larger bandwidth in the ferromagnetic state than in the paramagnetic state. Hence, when the carrier concentration n is small, the ferromagnetic state has lower energy than the paramagnetic state. The energy gain initially increases with increasing n and reaches a maximum at n x / 2 . Then it gradually decreases and finally vanishes at n x . The gain of kinetic energy results in ferromagnetism below a certain temperature. This implies that a double-exchange (DE)-like mechanism for ferromagnetism is operative in the impurity band. Assuming that the impurity band has the same shape as the model DOS [defined by Equation (27)] in the limit of the strong exchange interaction, the bandwidth of the magnetic impurity band is estimated to be 2 x Δ when S z = S and 2 x Δ when S z = 0 . The results for T c based on this assumption are inserted in Figure 9 and Figure 12 as “LIMIT”. The maximum T c is estimated to be
k B T c = 2 ( 2 2 ) 9 π x Δ
at n = x / 2 . The maximum T c ( = 0 . 0093 Δ ) for x = 0 . 05 is indicated by an arrow in Figure 10.
Figure 11. The DOS of the impurity band in the case of I S / Δ = 1 . 0 and E M = 0 . The thick, thin, and dotted lines represent the cases of S z / S = 1 . 0 , 0 . 5 , and 0 . 0 , respectively. Dot A indicates the impurity level for x 0 . From Takahashi and Kubo [41].
Figure 11. The DOS of the impurity band in the case of I S / Δ = 1 . 0 and E M = 0 . The thick, thin, and dotted lines represent the cases of S z / S = 1 . 0 , 0 . 5 , and 0 . 0 , respectively. Dot A indicates the impurity level for x 0 . From Takahashi and Kubo [41].
Materials 03 03740 g011
It is worth noting that the Zener double-exchange mechanism for ferromagnetism is usually understood to be effective only when the exchange energy is larger than the width of the carrier band (or | I S | 2 Δ ) in the case where magnetic ions sit on every site [54]. In the present case, the exchange energy is not greater than the width of the model band. Nevertheless, a DE-like mechanism functions because the magnetic impurity bandwidth is smaller than the exchange energy (or | I S | 2 x Δ ).
From the n dependence of T c shown in Figure 9, we conclude that the DE-like mechanism in a magnetic impurity band becomes dominant when | I S | / Δ 0 . 7 if E M = 0 . On the other hand, the DE-like mechanism is not relevant to the ferromagnetism in II-VI DMSs because the magnetic impurity level does not appear as illustrated in Figure 5 and Figure 6, in which the parameters | I S | / Δ = 0 . 4 and E M = 0 were employed.
Next, we study the role of the attractive potential in order to elucidate the origin of the carrier-induced ferromagnetism in III-V-based DMSs. We have already pointed out the similarity of the lower-energy part of the DOSs between Figure 4 (strong exchange interaction) and Figure 8 (strong attractive potential). The similarity is due to the fact that the impurity level has the same energy, which is determined by the effective attractive potential E M + I S . From the strong similarity in the low-energy part of the DOS, we may expect that ferromagnetism occurs through the same mechanism in both cases. In Figure 12 the effect of the nonmagnetic potential E M on T c is presented for I S fixed at 0 . 4 Δ . The impurity level appears for E M < 0 . 1 Δ in this case. When E M 0 . 0 the T c stays low and n c is much larger than x, while for E M 0 . 2 Δ a high T c is realized and n c is less than x. In the latter region the DE-like mechanism becomes operative. The criterion for the DE-like mechanism to operate is roughly estimated to be I S + 0 . 4 E M 0 . 6 Δ . The result suggests that the DE-like mechanism can be operative when an attractive potential assists in the production of an impurity band even though the exchange interaction is not particularly strong.
Figure 12. The result for T c / Δ as a function of n for various values of E M / Δ with x = 0 . 05 and I S = 0 . 4 Δ . The result denoted by `LIMIT’ is included (see text). From Takahashi and Kubo [42].
Figure 12. The result for T c / Δ as a function of n for various values of E M / Δ with x = 0 . 05 and I S = 0 . 4 Δ . The result denoted by `LIMIT’ is included (see text). From Takahashi and Kubo [42].
Materials 03 03740 g012

3.7. Specific Results for (Ga,Mn)As

Here we discuss (Ga,Mn)As, which has attracted much attention in recent years owing to its so-called carrier-induced ferromagnetism. Though the microscopic mechanism for carrier-induced ferromagnetism is still controversial, the following properties seem to be generally accepted for (Ga,Mn)As: (i) Mn ions substitute randomly for Ga cations in the zincblende structure [26]. (ii) A Mn ion in GaAs gives rise to an acceptor level at about 0.113 eV above the valence band [55]. (iii) The Mn ion has highly localized d states with a magnetic moment of 5 μ B (or S = 5 / 2 ) [55,56,57]. (iv) The Mn-induced states near the Fermi energy play a key role in the origin of ferromagnetism. According to photoemission studies [58,59,60], X-ray absorption spectroscopy [61], and band calculations [62,63], these states are mainly created in As 4p orbits. (v) The p-d exchange interaction between the As 4p hole and the localized d spin is antiferromagnetic [64,65], and its amplitude is not very different from that in II-VI-based DMSs [60,66]. (vi) As antisite defects (As ions sitting on Ga lattice site) and Mn interstitials are common in semiconductor samples grown by low-temperature molecular beam epitaxy [27,67,68,69,70].
Mn interstitials act as double donors. Many holes may be trapped not at Mn acceptors but at such defects, although we may expect that one hole is donated by a Mn atom. The density of the holes and that of the Mn ions are therefore regarded as separate sample-dependent quantities that are to be determined experimentally.
For Ga 1 x Mn x As, we set the (p-) bandwidth 2 Δ = 4 eV from the band calculations [62,63], and take I S = 0 . 4 Δ , which corresponds to N 0 β = 0 . 64 eV. Then we determine E M to be 0 . 3 Δ so as to yield the acceptor energy of 0.113 eV ( = 0 . 057 Δ ) [55]. The x dependence of the DOS is shown in Figure 13; S z = 0 (a) and S z = S (b). The present model parameters lead to an impurity level at the energy of E a = 1 . 057 Δ in the dilute limit ( x 0). With an increase in x, an impurity band forms, and for x 0 . 02 it merges into the host valence band. The results consistently explain the experimental observation of impurity-band-like states [59,71] and the insulator-metal transition at x 0 . 03 [72,73].
Figure 13. The lower-energy part of the DOS of Ga 1 x Mn x As for various values of x: (a) S z = 0 and (b) S z = 0 . The solid lines represent the down-spin carrier and the dotted lines represent the up-spin carrier. The impurity level E a = 1 . 057 Δ is indicated by a dot on the line x = 0 . From Takahashi and Kubo [42].
Figure 13. The lower-energy part of the DOS of Ga 1 x Mn x As for various values of x: (a) S z = 0 and (b) S z = 0 . The solid lines represent the down-spin carrier and the dotted lines represent the up-spin carrier. The impurity level E a = 1 . 057 Δ is indicated by a dot on the line x = 0 . From Takahashi and Kubo [42].
Materials 03 03740 g013
In Figure 14 and Figure 15, the results for x = 0 . 05 are presented. The lower band tail shown in the left panel of Figure 14 is strongly affected by the change in S z . The results of species-resolved analysis shown in Figure 14a,b reveal that the change in the band tail is mainly ascribed to the change in the local DOS at the Mn site. This result indicates that a carrier at the band tail usually stays at Mn sites in spite of the small x ( = 0 . 05 ) . The present result is in sharp contrast with the free-carrier picture, which is the premise for the application of the RKKY model [74], but is consistent with the nearly bound hole picture deduced on the basis of infrared optical absorption measurement [75,76]. The carrier’s spin is tightly coupled with the localized spin. In the left panel of Figure 15, we plot the Curie temperature T c as a function of n. It has been reported that ferromagnetism with T c = 110 K is realized in Ga 1 x Mn x As with x = 0 . 053 when n is 30% of the nominal concentration (x) of Mn [26]. The agreement of the present result with the experimental observation is satisfactory. The T dependence of the magnetization is presented in the right panel of Figure 15 for various values of n. The present result is consistent with the experimentally obtained magnetization [26]. We have also verified that the result obtained by applying DMFT to the present model is almost the same as the present result [44].
Figure 14. The results for Ga 1 x Mn x As with x = 0 . 05 . Left panel: Low-energy part of the DOS shown for various values of S z / S . The solid lines represent the down-spin carrier and the dotted lines represent the up-spin carrier. The arrow indicates the Fermi level for n = x ( = 0 . 05 ) . The impurity level E a = 1 . 057 Δ is indicated by a dot on the horizontal line S z = 0 . Right panel: (a) Ga-site component of the DOS, ( 1 x ) D A ( ω ) Δ and ( 1 x ) D A ( ω ) Δ and (b) Mn-site component of the DOS, x D M ( ω ) Δ and x D M ( ω ) Δ . The thick, thin, and dotted lines represent the cases of S z / S = 1 . 0 , 0 . 5 , and 0 . 0 , respectively. From Takahashi and Kubo [42].
Figure 14. The results for Ga 1 x Mn x As with x = 0 . 05 . Left panel: Low-energy part of the DOS shown for various values of S z / S . The solid lines represent the down-spin carrier and the dotted lines represent the up-spin carrier. The arrow indicates the Fermi level for n = x ( = 0 . 05 ) . The impurity level E a = 1 . 057 Δ is indicated by a dot on the horizontal line S z = 0 . Right panel: (a) Ga-site component of the DOS, ( 1 x ) D A ( ω ) Δ and ( 1 x ) D A ( ω ) Δ and (b) Mn-site component of the DOS, x D M ( ω ) Δ and x D M ( ω ) Δ . The thick, thin, and dotted lines represent the cases of S z / S = 1 . 0 , 0 . 5 , and 0 . 0 , respectively. From Takahashi and Kubo [42].
Materials 03 03740 g014
Here, we briefly discuss the origin and the mechanism of the carrier-induced ferromagnetism of (Ga,Mn)As. The results for the n dependence of T c shown in the left panel of Figure 15 are well explained by the DE-like mechanism described above. Although the impurity band is not separated in the case of x 0 . 05 , the carriers in the band tail of Ga 1 x Mn x As have such a high local carrier density at Mn sites that the carrier spins couple strongly to the fluctuating localized spins. Thus, the hopping of carriers among Mn sites causes the ferromagnetic ordering of the localized spins through the DE-like mechanism.
Since Zener originally proposed the DE interaction for (La,Ca)MnO 3 , where 3d holes hop among the magnetic ions located at the regular lattice sites [54], it might be understood that the DE mechanism is only relevant to the hopping of 3d holes in (Ga,Mn)As [77]. The DE mechanism for ferromagnetism, in fact, works quite generally. The only condition required for the mechanism is very strong spin coupling between carrier spins and localized spins. If this is satisfied, carriers may have any character and the localized spins can be arranged randomly. In the case of III-V-based DMSs, carriers are considered to have 4p character [60,61] and the strength of the p-d exchange interaction is not very different from that in II-VI-based DMSs [66]. The Coulomb interaction between the carrier (hole) and a Mn 2 + ion (acceptor center), however, promotes the formation of a magnetic impurity band, and strong spin coupling is realized in the magnetic impurity band and/or in the band tail. Therefore, the DE mechanism induces ferromagnetism. We call this mechanism the DE-like mechanism to avoid confusion with the argument assuming d holes [77]. Note that no Mn 3 + ( d 4 configuration) states have been experimentally detected in (Ga,Mn)As [55,57,60,64,78]. All these experimental observations suggest that the fixed valence state Mn 2 + ( S = 5 / 2 ) is realized in (Ga,Mn)As.
Figure 15. The results for Ga 1 x Mn x As with x = 0 . 05 . Left panel: Curie temperature T c as a function of n. Right panel: Magnetization as a function of the temperature T for various values of n. From Takahashi and Kubo [42] and Takahashi et al. [44].
Figure 15. The results for Ga 1 x Mn x As with x = 0 . 05 . Left panel: Curie temperature T c as a function of n. Right panel: Magnetization as a function of the temperature T for various values of n. From Takahashi and Kubo [42] and Takahashi et al. [44].
Materials 03 03740 g015
Although the exchange interaction between p holes and d spins has been experimentally proved to be antiferromagnetic [60,64,65], in the early stage of research, the ferromagnetic exchange interaction was reported on the basis of polarized magnetoreflection measurement [79]. When an impurity band exists, the optically observed band edge is not necessarily the band edge of the impurity band. The present result for x = 0 . 005 is shown in Figure 16. When x = 0 . 5 % , the impurity band is separate from the host band, as shown in Figure 16(a). The optical absorption spectrum corresponds to the k = 0 component of the DOS. Therefore, the optical absorption spectrum takes negligible values in the impurity band, because the impurity band is constituted from the wide range of k space. The optical band edge, ω p , lies almost at the bottom of the host band. Thus, although a negative I S is assumed, the optical band edge with up (down) spin shifts to the low (high)-energy side with increasing S z , as shown in Figure 16(b). Hence, the direction of the shift in the optical band edge is opposite to the direction obtained by the VCA. Note that ferromagnetic spin coupling is realized near the bottom of the host band edge. In magnetoreflection measurement, not the shift of the impurity band edge but the shift of the host band edge was detected. A simple calculation based on the present approach with N 0 β = 0 . 64 eV predicts that the magnetoreflection measurement will deduce the apparent value of N 0 β to be + 1 . 3 eV, which is consistent with the experimental observation [79].
Figure 16. The result for Ga 1 x Mn x As with x = 0 . 005 : (a) low-energy part of the DOS D ( ω ) , (b) optical absorption band A ( ω ) in arbitrary units (arb. units). The solid lines represent the down-spin carrier and the dotted lines represent the up-spin carrier. Note that the energy of the bottom of the model band is ω = Δ . Along the upper horizontal axis of (a), energies for Δ = 2 eV are graduated in eV, where the origin of the energy, 0.0 eV, is taken at ω = Δ . The impurity level E a = 1 . 057 Δ (or 0 . 113 eV) is indicated by a dot on the line S z = 0 . From Takahashi [43].
Figure 16. The result for Ga 1 x Mn x As with x = 0 . 005 : (a) low-energy part of the DOS D ( ω ) , (b) optical absorption band A ( ω ) in arbitrary units (arb. units). The solid lines represent the down-spin carrier and the dotted lines represent the up-spin carrier. Note that the energy of the bottom of the model band is ω = Δ . Along the upper horizontal axis of (a), energies for Δ = 2 eV are graduated in eV, where the origin of the energy, 0.0 eV, is taken at ω = Δ . The impurity level E a = 1 . 057 Δ (or 0 . 113 eV) is indicated by a dot on the line S z = 0 . From Takahashi [43].
Materials 03 03740 g016
The difference between the characters of a hole in a II-VI-based DMS and a hole in a III-VI-based DMS is illustrated in Figure 17. The hole in a II-VI-based DMS can move freely over many sites while undergoing exchange interactions with the d spin on Mn sites. On the contrary, the hole in a III-V-based DMS moves while hopping from a Mn site to another Mn site because the attractive Coulomb potential makes it favorable for the hole to remain at Mn sites.
Figure 17. An illustration of the difference between the characters of a hole in a II-VI-based DMS and in a III-VI-based DMS.
Figure 17. An illustration of the difference between the characters of a hole in a II-VI-based DMS and in a III-VI-based DMS.
Materials 03 03740 g017

4. Summary and Concluding Comments

Throughout this review article, we considered a simple model with fewer physical parameters [80]. The advantage of this approach is that the basic physics of systems can be explained in a simple way, regardless of the strong material dependence, although realistic electronic-structure calculations such as first-principles calculations based on density functional theory [16,17,81,82] or numerical methods such as the quantum Monte Carlo method [83,84] may give more realistic results.
First, applying the dynamical CPA to the s-f model, we showed the numerical results for FMSs such as EuO and EuS. The results for the DOS and the energy of the bottom of the band were given for various exchange energies and temperatures. Not only the dependence of the band edge on the temperature but also the magnitude of the redshift can be explained satisfactorily. We pointed out that the electron-spin polarization cannot reach 100% due to the quantum effect coming from the finiteness of the magnitude of the f spin [34]. We also add that the temperature dependence of the electrical resistivity of a degenerate FMS can be consistently explained except around T c [35]; the dynamical CPA becomes inefficient when the short-range order of f spins develops near T c [85,86,87,88] .
Next, we considered a simple model for carriers in the A 1 x Mn x B-type DMS. The model includes not only the exchange interaction but also the nonmagnetic local potential at the magnetic Mn ion site. Based on the results obtained by applying the dynamical CPA to the simple model, we discussed the carrier states of three typical cases: cases with strong and moderate exchange interactions in the absence of nonmagnetic potentials, and the case with strong attractive nonmagnetic potentials in addition to moderate exchange interaction. We showed that the mechanism of carrier-induced ferromagnetism changes from the RKKY-type mechanism to the DE-like mechanism with the increase in the amplitude of the exchange interaction. Furthermore, the DE-like mechanism can be operative when an attractive potential assists in the production of an impurity band even though the exchange interaction is not particularly strong.
Carrier states in II-VI-based DMSs are well described by the present model with moderate exchange interaction. The results obtained by applying the dynamical CPA to the model explain the anomalous x dependence (i.e., bowing) of the energy gap of wide-band-gap II-VI-based DMSs [38], the apparent enhancement of the p-d exchange interaction, and the asymmetric splitting of exciton states in Cd 1 x Mn x S [39,40]; these effects cannot be explained by the VCA.
In the case of III-VI-based DMSs such as Ga 1 x Mn x As and In 1 x Mn x As, Coulomb attractive potentials assist the formation of the impurity band and/or band tail and promote the DE-like mechanism although the exchange interaction is not sufficiently strong. Setting the parameters so as to yield an experimentally observed impurity level, we calculated the DOS for various x and magnetizations, as well as the dependence of T c on n. The results for the local DOS suggest that the exchange coupling between a carrier and localized spins at Mn sites is very strong [41,42]. Thus, we conclude that the carrier-induced ferromagnetism of Ga 1 x Mn x As is ascribed to a DE-like mechanism realized in the magnetic impurity band/or in the band tail. The present results also reveal the difference between the optical band edges in II-VI-based DMSs and in III-V-based DMSs [43].
However, we have to point out that the many-body effect [89] and clustering effect [90] are not considered in our CPA. With the increase in Mn fraction x, the effects of the direct antiferromagnetic (AF) superexchange interaction between neighboring Mn impurities become important [90]. It should also be considered that the effect of the Coulomb attractive potential E M becomes less important as the carrier density increases [45]. As noted in Section 3.7, in real samples of (Ga,Mn)As, Mn ions do not always substitute for Ga cations. Some fraction of Mn may reside in the interstitial lattice sites and act as double donors [68]; many holes may be trapped at such defects, reducing the concentration of holes. Furthermore, As antisites may induce a so-called disordered local moment configuration, where only part of the randomly distributed Mn atoms are ferromagnetically aligned while the rest of the Mn atoms have magnetic moments oriented antiparallel to each other [69]. These effects are beyond the scope of our simple model and approach.
Recently, Gd-doped EuO has been highlighted for potential use in spintronics devices owing to its with high T c [5,6,8]. The mechanism of the enhancement of T c of Gd-doped EuO is a controversial problem. This problem can be treated by applying the dynamical CPA to a simple model adapted for Eu 1 x Gd x O. The metal-insulator transition in Eu-rich EuO is still under discussion [15]. The model for Eu 1 y V y O, where V denotes a vacancy of an O ion, is the same as that for A 1 x Mn x B-type DMSs, and can act as a substitute for EuV y O 1 y , which is the most realistic model of Eu-rich EuO. An extension of the dynamical CPA to these problems is now in progress and will be published separately.

Acknowledgement

The author appreciates the financial support from his elder sister Hisako Sasaki during his childhood, and would like to sincerely thank her for her continuous warm encouragement throughout her life.

References

  1. Wacter, P. Europium chalcogenides: EuO, EuS, EuSe and EuTe. In Handbook on Physics and Chemistry of Rare Earths; Gschneidner, K.A., Eyring, L., Eds.; North-Holland Publishing Company: Amsterdam, The Netherlands, 1979; pp. 507–576. [Google Scholar]
  2. Mauger, A.; Godart, C. The magnetic, optical, and transport properties of representatives of a class of magnetic semiconductors: The europium chalcogenides. Phys. Rep. 1986, 141, 51–176. [Google Scholar] [CrossRef]
  3. von Molnar, S.; Kasuya, T. Evidence of band conduction and critical scattering in dilute Eu-chalcogenide alloys. Phys. Rev. Lett. 1968, 21, 1757–1761. [Google Scholar] [CrossRef]
  4. Oliver, M.R.; Dimmock, J.O.; McWhorter, A.L.; Reed, T.B. Conductivity studies in europium oxide. Phys. Rev. B 1971, 5, 1078–1098. [Google Scholar] [CrossRef]
  5. Matsumoto, T.; Yamaguchi, K.; Yuri, M.; Kawaguchi, K.; Koshizaki, N.; Yamada, K. Preparation of Gd-doped EuO1−x thin films and the magnetic and magneto-transport properties. J. Phys-Condens. Matter 2004, 16, 6017–6028. [Google Scholar] [CrossRef]
  6. Ott, H.; Heise, S.J.; Sutarto, R.; Hu, Z.; Chang, C.F.; Hsieh, H.H.; Lin, H.-J.; Chen, C.T.; Tjeng, L.H. Soft x-ray magnetic circular dichroism study on Gd-doped EuO thin films. Phys. Rev. B 2006, 73, 094407. [Google Scholar] [CrossRef]
  7. Santos, T.S.; Moodera, J.S.; Raman, K.V.; Negusse, E.; Holroyd, J.; Dvorak, J.; Liberati, M.; Idzerda, Y.U.; Arenholz, E. Determining exchange splitting in a magnetic semiconductor by spin-filter tunneling. Phys. Rev. Lett. 2008, 101, 147201. [Google Scholar] [CrossRef] [PubMed]
  8. Sutarto, R.; Altendorf, S.G.; Coloru, B.; Moretti Sala, M.; Haupricht, T.; Chang, C.F.; Hu, Z.; Schusler-Langeheine, C.; Hollmann, N.; Kierspel, H.; Mydosh, J.A.; Hsieh, H.H.; Lin, H.-J.; Chen, C.T.; Tjeng, L.H. Epitaxy, stoichiometry, and magnetic properties of Gd-doped EuO films on YSZ (001). Phys. Rev. B 2009, 80, 085308. [Google Scholar] [CrossRef]
  9. Kimura, S.; Ito, T.; Miyazaki, H.; Mizuno, T.; Iizuka, T.; Takahashi, T. Electronic inhomogeneity EuO: Possibility of magnetic polaron states. Phys. Rev. B 2008, 78, 052409. [Google Scholar] [CrossRef]
  10. Comment, A.; Ansermet, J.-P.; Slichter, C.P.; Rho, H.; Snow, C.S.; Cooper, S.L. Magnetic properties of pure and Gd-doped EuO probed by NMR. Phys. Rev. B 2005, 72, 014428. [Google Scholar] [CrossRef]
  11. Miyazaki, H.; Ito, T.; Im, H.J.; Yagi, S.; Kato, M.; Soda, K.; Kimura, S. Direct observation of momentum-dependent exchange interaction in a Heisenberg ferromagnet. Phys. Rev. Lett. 2009, 102, 227203. [Google Scholar] [CrossRef] [PubMed]
  12. Steeneken, P.G.; Tjeng, L.H.; Elfimov, I.; Sawatzky, G.A.; Ghiringhelli, G.; Brookes, N.B.; Huang, D.-J. Exchange splitting and charge carrier spin polarization in EuO. Phys. Rev. Lett. 2002, 88, 047201. [Google Scholar] [CrossRef] [PubMed]
  13. Souza-Neto, N.M.; Haskel, D.; Tseng, Y.-C.; Lapertot, G. Pressure-induced electronic mixing and enhancement of ferromagnetic ordering in EuX (X=Te, Se, S, O) magnetic semiconductors. Phys. Rev. Lett. 2009, 102, 057206. [Google Scholar] [CrossRef] [PubMed]
  14. Schiller, R.; Müller, W.; Nolting, W. Kondo-lattice model: Application to the temperature-dependent electronic structure of EuO(100) films. Phys. Rev. B 2001, 64, 134409. [Google Scholar] [CrossRef]
  15. Sinjukow, P.; Nolting, W. Fully self-consistent determination of transport properties in Eu-rich EuO. Phys. Rev. B 2004, 69, 214432. [Google Scholar] [CrossRef]
  16. Ghosh, D.B.; De, M.; De, S.K. Electronic structure and magneto-optical properties of magnetic semiconductors: Europium monochalcogenides. Phys. Rev. B 2004, 70, 115211. [Google Scholar] [CrossRef]
  17. Ingle, N.J.C.; Elfimov, I.S. Influence of epitaxial strain on the ferromagnetic semiconductor EuO: First-principles calculations. Phys. Rev. B 2008, 77, 121202. [Google Scholar] [CrossRef]
  18. Arnold, M.; Kroha, J. Simultaneous ferromagnetic metal-semiconductor transition in electron-doped EuO. Phys. Rev. Lett. 2008, 100, 046404. [Google Scholar] [CrossRef] [PubMed]
  19. Nolting, W.; Java, M.S.; Rex, S. Magnetic polaron in ferro- and antiferromagnetic semiconductors. Phys. Rev. B 1996, 54, 14455–14466. [Google Scholar] [CrossRef]
  20. Nolting, W.; Rex, S.; Java, M.S. Magnetism and electronic structure of local moment ferromagnet. J. Phys-Condens. Matter 1997, 9, 1301–1330. [Google Scholar] [CrossRef]
  21. Bryksa, V.; Nolting, W. Disordered Kondo-lattice model: Extension of coherent potential approximation. Phys. Rev. B 2008, 78, 064417. [Google Scholar] [CrossRef]
  22. Tang, G.; Nolting, W. Carrier-induced ferromagnetism in diluted local moment systems. Phys. Rev. B 2007, 75, 024426. [Google Scholar] [CrossRef]
  23. Furdya, J.K. Diluted magnetic semiconductors. J. Appl. Phys. 1988, 64, R29–R64. [Google Scholar] [CrossRef]
  24. Haas, K.C. Band structure and theory of magnetic interactions. In Semimagnetic Semiconductors and Diluted Magnetic Semiconductors; Averous, M., Balkanski, M., Eds.; Plenum: New York, NY, USA, 1991; pp. 59–82. [Google Scholar]
  25. Lascaray, J.P. II1−xMnxVI semimagnetic semiconductors. In Semimagnetic Semiconductors and Diluted Magnetic Semiconductors; Averous, M., Balkanski, M., Eds.; Plenum: New York, NY, USA, 1991; pp. 169–190. [Google Scholar]
  26. Ohno, H. Properties of ferromagnetic III-V semiconductors. J. Magn. Magn. Mater. 1999, 200, 110–129. [Google Scholar] [CrossRef]
  27. Jungwirth, T.; Masek, J.; Kucera, J.; MacDonald, H. Theory of ferromagnetic (III,Mn)V semiconductors. Rev. Mod. Phys. 2006, 78, 809–864. [Google Scholar] [CrossRef]
  28. Ehrenreich, H.; Schwartz, L.M. The electronic structure of alloys. In Solid State Physics, Advances in Research and Application; Ehrenreich, H., Seitz, F., Turnbull, D., Eds.; Academic Press: New York, NY, USA, 1976; Volume 31, pp. 149–286. [Google Scholar]
  29. Gonis, A. Green functions for ordered and disordered systems. In Studies in Mathematical Physics; van Groesen, E., de Jager, E.M., Eds.; North-Holland: Amsterdam, The Netherlands, 1992; Volume 4. [Google Scholar]
  30. Onodera, Y.; Toyozawa, Y. Persistence and amalgamation types in the electronic structure of mixed crystals. J. Phys. Soc. Jpn. 1968, 24, 341–355. [Google Scholar] [CrossRef]
  31. Rangette, A.; Yanase, A.; Kübker, J. CPA treatment of the s-d model at high temperatures. Solid State Commun. 1973, 12, 171–174. [Google Scholar] [CrossRef]
  32. Kubo, K. Electronic states in magnetic semiconductors–An extension of CPA to random spin systems–. J. Phys. Soc. Jpn. 1974, 36, 32–38. [Google Scholar] [CrossRef]
  33. Takahashi, M.; Mitsui, K. Single-site approximation for the s-f model in ferromagnetic semiconductors. Phys. Rev. B 1996, 54, 11298–11304. [Google Scholar] [CrossRef]
  34. Takahashi, M. Electron-spin polarization in ferromagnetic semiconductors. Phys. Rev. B 1997, 56, 7389–7394. [Google Scholar] [CrossRef]
  35. Takahashi, M.; Mitsui, K. Electron scattering due to fluctuating localized spins in degenerate ferromagnetic semiconductors. J. Magn. Magn. Mater. 1998, 182, 329–334. [Google Scholar] [CrossRef]
  36. Takahashi, M. Conduction electron band in antiferromagnetic semiconductors. Phys. Rev. B 1997, 55, 6950–6956. [Google Scholar] [CrossRef]
  37. Takahashi, M.; Nolting, W. Single-site approximation for the s-f model of antiferromagnetic semiconductors. J. Magn. Magn. Mater. 2001, 237, 225–230. [Google Scholar] [CrossRef]
  38. Takahashi, M. Coherent-potential approach to magnetic and chemical disorder in diluted magnetic semiconductors. Phys. Rev. B 1999, 60, 15858–15865. [Google Scholar] [CrossRef]
  39. Takahashi, M. Coherent potential approach to exchange-induced band splitting in diluted magnetic semiconductors under a saturating magnetic field. J. Phys-Condens. Matter 2001, 13, 3433–3442. [Google Scholar] [CrossRef]
  40. Takahashi, M. Asymmetric splitting of exciton states in diluted magnetic semiconductor Cd1−xMnxS. J. Phys. Soc. Jpn. 2001, 70, 2224–2225. [Google Scholar] [CrossRef]
  41. Takahashi, M.; Kubo, K. Mechanism of carrier-induced ferromagnetism in magnetic semiconductors. Phys. Rev. B 2002, 66, 153202. [Google Scholar] [CrossRef]
  42. Takahashi, M.; Kubo, K. Carrier states and ferromagnetism in diluted magnetic semiconductors. J. Phys. Soc. Jpn. 2003, 72, 2866–2879. [Google Scholar] [CrossRef]
  43. Takahashi, M. Optical band edge of diluted magnetic semiconductors. Phys. Rev. B 2004, 70, 035207. [Google Scholar] [CrossRef]
  44. Takahashi, M.; Furukawa, N.; Kubo, K. Mechanism of carrier-induced ferromagnetism in diluted magnetic semiconductors. J. Magn. Magn. Mater. 2003, 272, 2021–2022. [Google Scholar] [CrossRef]
  45. Popescu, F.; Şen, C.; Dagotto, E.; Moreo, A. Crossover from impurity to valence band in diluted magnetic semiconductors: Role of Coulomb attraction by acceptors. Phys. Rev. B 2007, 76, 085206. [Google Scholar] [CrossRef]
  46. Hoang, A.T. Optical properties of diluted magnetic semiconductors in coherent potential approximation. Physica B 2008, 403, 1803–1807. [Google Scholar] [CrossRef]
  47. Nolting, W.; Oreś, A.M. Conduction-band structure of a ferromagnetic semiconductor. Phys. Rev. B 1980, 22, 6184–6195. [Google Scholar] [CrossRef]
  48. Lascaray, J.P.; Deruelle, M.C.D.; Coquillat, D. Magnetization and magnetoreflectivity measurements in Zn1−xMnxTe with 0.25 ≤ x ≤ 0.71. Phys. Rev. B 1987, 35, 675–677. [Google Scholar] [CrossRef]
  49. Lascaray, J.P.; Coquillat, D.; Deportes, J.; Bhattacharjee, A.K. Zeeman splitting of exciton and magnetization in Cd1−xMnxTe: Anomalous behavior at high x. Phys. Rev. B 1988, 38, 7602–7606. [Google Scholar] [CrossRef]
  50. Takahashi, M.; Kubo, K. Limitation of mean field approximation for carrier states and ferromagnetism in diluted magnetic semiconductors. J. Phys. Soc. Jpn. 2005, 74, 1642–1643. [Google Scholar] [CrossRef]
  51. Chattopadhyay, A.; Das Sarma, S.; Millis, A.J. Transition temperature of ferromagnetic semiconductors: A dynamical mean field study. Phys. Rev. Lett. 2001, 87, 227202. [Google Scholar] [CrossRef] [PubMed]
  52. Yagi, M.; Kayanuma, Y. Theory for carrier-induced ferromagnetism in diluted magnetic semiconductors. J. Phys. Soc. Jpn. 2002, 71, 2010–2018. [Google Scholar] [CrossRef]
  53. Calderón, M.J.; Gomez-Santos, G.; Brey, L. Impurity-semiconductor band hybridization effects on the critical temperature of diluted magnetic semiconductors. Phys. Rev. B 2002, 66, 075218. [Google Scholar] [CrossRef]
  54. Zener, C. Interaction between the d-shells in the transition metals. II. Ferromagnetic compounds of manganese with perovskite structure. Phys. Rev. 1951, 82, 403–405. [Google Scholar] [CrossRef]
  55. Linnarsson, M.; Janzen, E.; Monemar, B.; Kleverman, M.; Thilderkvist, A. Electronic structure of the GaAs:MnGa center. Phys. Rev. B 1997, 55, 6938–6944. [Google Scholar] [CrossRef]
  56. Iye, Y.; Oiwa, A.; Endo, A.; Katsumoto, S.; Matsukura, F.; Shen, A.; Ohno, H.; Munekata, H. Metal-insulator transition and magnetotransport in III-V compound diluted magnetic semiconductors. Mater. Sci. Eng. B 1999, 63, 88–95. [Google Scholar] [CrossRef]
  57. Ohldag, H.; Solinus, V.; Hillebrecht, F.U.; Goedkoop, J.B.; Finazzi, M.; Matsukura, F.; Ohno, H. Magnetic moment of Mn in the ferromagnetic semiconductor (Ga0.98Mn0.02)As. Appl. Phys. Lett. 2000, 76, 2928–2930. [Google Scholar] [CrossRef]
  58. Okabayashi, J.; Kimura, A.; Mizokawa, T.; Fujimori, A.; Hayashi, T.; Tanaka, M. Mn 3d partial density of states in Ga1−xMnxAs studied by resonant photoemission spectroscopy. Phys. Rev. B 1999, 59, R2486–R2489. [Google Scholar] [CrossRef]
  59. Okabayashi, J.; Kimura, A.; Rader, O.; Mizokawa, T.; Fujimori, A.; Hayashi, T.; Tanaka, M. Angle-resolved photoemission study of Ga1−xMnxAs. Phys. Rev. B 2001, 64, 125304. [Google Scholar] [CrossRef]
  60. Okabayashi, J.; Kimura, A.; Rader, O.; Mizokawa, T.; Fujimori, A.; Hayashi, T.; Tanaka, M. Core-level photoemission study of Ga1−xMnxAs. Phys. Rev. B 1998, 58, R4211–R4214. [Google Scholar] [CrossRef]
  61. Ishiwata, Y.; Watanabe, M.; Eguchi, R.; Takeuchi, T.; Harada, Y.; Chainani, A.; Shin, S.; Hayashi, T.; Hashimoto, Y.; Katsumoto, S.; Iye, Y. Manganese concentration and low-temperature annealing dependence of Ga1−xMnxAs by x-ray absorption spectroscopy. Phys. Rev. B 2002, 65, 233201. [Google Scholar] [CrossRef]
  62. Park, J.H.; Kwon, S.K.; Min, B.I. Electronic structures of III-V based ferromagnetic semiconductors: half-metallic phase. Physica B 2001, 281, 703–704. [Google Scholar] [CrossRef]
  63. Shirai, M.; Ogawa, T.; Kitagawa, I.; Suzuki, N. Band structures of zinc-blende-type MnAs and (MnAs)1(GaAs)1 superlattice. J. Magn. Magn. Mater. 1998, 177, 1383–1384. [Google Scholar] [CrossRef]
  64. Szczytko, J.; Mac, W.; Twardowski, A.; Matsukura, F.; Ohno, H. Antiferromagnetic p-d exchange in ferromagnetic Ga1−xMnxAs epilayers. Phys. Rev. B 1999, 59, 12935–12939. [Google Scholar] [CrossRef]
  65. Ando, K.; Hayashi, T.; Tanaka, M.; Twardowski, A. Magneto-optic effect of the ferromagnetic diluted magnetic semiconductor Ga1−xMnxAs. J. Appl. Phys. 1998, 83, 6548–6550. [Google Scholar] [CrossRef]
  66. Okabayashi, J.; Mizokawa, T.; Sarma, D.D.; Fujimori, A.; Slupinski, T.; Oiwa, A.; Munekata, H. Electronic structure of In1−xMnxAs studied by photoemission spectroscopy: Comparison with Ga1−xMnxAs. Phys. Rev. B 2002, 65, 161203. [Google Scholar] [CrossRef]
  67. Hayashi, T.; Hashimoto, Y.; Katsumoto, S.; Iye, Y. Effect of low-temperature annealing on transport and magnetism of diluted magnetic semiconductor (Ga,Mn)As. Appl. Phys. Lett. 2001, 78, 1691–1693. [Google Scholar] [CrossRef]
  68. Yu, K.M.; Walukiewicz, W.; Wojtowicz, T.; Kuryliszyn, I.; Liu, X.; Sasaki, Y.; Furdyna, J.K. Effect of the location of Mn sites in ferromagnetic Ga1−xMnxAs on its Curie temperature. Phys. Rev. B 2002, 65, 201303. [Google Scholar] [CrossRef]
  69. Korzhavyi, P.A.; Abrikosov, I.A.; Smirnova, E.A.; Bergqvist, L.; Mohn, P.; Mathieu, R.; Svedlindh, P.; Sadowski, J.; Isaev, E.I.; Vekilov, Yu.Kh.; Eriksson, O. Defect-induced magnetic structure in (Ga1−xMnx)As. Phys. Rev. Lett. 2002, 88, 187202. [Google Scholar] [CrossRef] [PubMed]
  70. Timm, C. Disorder effects in diluted magnetic semiconductors. J. Phys-Condens. Matter 2003, 15, R1865–R1896. [Google Scholar] [CrossRef]
  71. Okabayashi, J.; Kimura, A.; Rader, O.; Mizokawa, T.; Fujimori, A.; Hayashi, T.; Tanaka, M. Electronic structure of GaMnAs studied by angle-resolved photoemission spectroscopy. Physica E 2001, 10, 192–195. [Google Scholar] [CrossRef]
  72. Oiwa, A.; Katsumoto, S.; Endo, A.; Hirakawa, M.; Iye, M.; Ohno, H.; Matsukura, F.; Shen, A.; Sugawara, Y. Nonmetal-metal-nonmetal transition and large negative magnetoresistance in (Ga,Mn)As/GaAs. Solid State Commun. 1997, 103, 209–213. [Google Scholar] [CrossRef]
  73. Oiwa, A.; Katsumoto, S.; Endo, A.; Hirakawa, M.; Iye, M.; Matsukura, F.; Shen, A.; Sugawara, Y.; Ohno, H. Low-temperature conduction and giant negative magnetoresistance in III-V-based diluted magnetic semiconductor. Physica B 1998, 249, 775–779. [Google Scholar] [CrossRef]
  74. Matsukura, F.; Ohno, H.; Shen, A.; Sugawara, Y. Transport properties and origin of ferromagnetism in (Ga,Mn)As. Phys. Rev. B 1998, 57, R2037–R2040. [Google Scholar] [CrossRef]
  75. Hirakawa, K.; Katsumoto, S.; Hayashi, T.; Hashimoto, Y.; Iye, Y. Double-exchange-like interaction in Ga1−xMnxAs investigated by infrared absorption spectroscopy. Phys. Rev. B 2002, 65, 193312. [Google Scholar] [CrossRef]
  76. Katsumoto, S.; Hayashi, T.; Hashimoto, Y.; Iye, Y.; Ishiwata, Y.; Watanabe, M.; Eguchi, R.; Takeuchi, T.; Harada, Y.; Shin, S.; Hirakawa, K. Magnetism and metal-insulator transition in III-V based diluted magnetic semiconductors. Mater. Sci. Eng. 2001, 84, 88–95. [Google Scholar] [CrossRef]
  77. Akai, H. Ferromagnetism and its stability in the diluted magnetic semiconductor (In, Mn)As. Phys. Rev. Lett. 1998, 81, 3002–3005. [Google Scholar] [CrossRef]
  78. Szczytko, J.; Twardowski, A.; Świa̧tek, K.; Palczewska, M.; Tanaka, M.; Hayashi, T.; Ando, K. Mn impurity in Ga1−xMnxAs epilayers. Phys. Rev. B 1999, 60, 8304–8308. [Google Scholar] [CrossRef]
  79. Szczytko, J.; Mac, M.; Stachow, A.; Twardowski, A.; Becla, P.; Tworzydlo, J. The s,p-d exchange interaction in GaAs heavily doped with Mn. Solid State Commun. 1996, 99, 927–931. [Google Scholar] [CrossRef]
  80. Kagami, K.; Takahashi, M.; Yasuda, C.; Kubo, K. Theory of diluted magnetic semiconductors: A minimum model. Sci. Technol. Adv. Mater. 2006, 7, 31–41. [Google Scholar] [CrossRef]
  81. Nazir, S.; Ikram, N.; Tanveer, M.; Shaukat, A.; Saeed, Y.; Reshak, A.H. Spin-polarized structural, electronic, and magnetic properties of diluted magnetic semiconductors Cd1−xMnxS and Cd1−xMnxSe in zinc blende phase. J. Phys. Chem. A 2009, 113, 6022–6027. [Google Scholar] [CrossRef] [PubMed]
  82. Sato, K.; Katayama-Yoshida, H. First-principles materials design for semiconductor spintronics. Semicond. Sci. Technol. 2002, 17, 367–376. [Google Scholar] [CrossRef]
  83. Sakai, O.; Suzuki, S.; Nishizawa, K. Study on the magnetic and transport properties of low density carrier ferromagnetic semiconductors. J. Phys. Soc. Jpn. 2001, 70, 1105–1116. [Google Scholar] [CrossRef]
  84. Ohe, J.; Tomoda, Y.; Bulut, N.; Arita, R.; Nakamura, K.; Maekawa, S. Combined approach of density function theory and quantum Monte Carlo method to electron correlation in diluted magnetic semiconductors. J. Phys. Soc. Jpn. 2009, 78, 083703. [Google Scholar] [CrossRef]
  85. Takahashi, M.; Mitsui, K.; Umehara, M. Conduction-electron states in ferromagnetic semiconductors above the Curie temperature. Phys. Rev. B 1993, 48, 17053–17063, Erratum: Phys. Rev. B 1994, 50, 17759–17760. [Google Scholar] [CrossRef]
  86. Takahashi, M.; Mitsui, K.; Umehara, M. Multiple-scattering approach to the s-f model in ferromagnetic semiconductors above the Curie temperature. Phys. Rev. B 1995, 52, 16313–16316. [Google Scholar] [CrossRef]
  87. Takahashi, M.; Kasuya, T. A theory for a self-trapped molecular magnetic polaron in ferromagnetic semiconductor. J. Phys. Soc. Jpn. 1983, 52, 3127–3137. [Google Scholar] [CrossRef]
  88. Takahashi, M.; Kasuya, T. Effect of electron-phonon interaction in the magnetic polaron. J. Phys. Soc. Jpn. 1983, 52, 3138–3143. [Google Scholar] [CrossRef]
  89. Edwards, D.M.; Green, A.C.M.; Kubo, K. Electronic structure and resistivity of the double exchange model. J. Phys-Condens. Matter 1999, 11, 2791. [Google Scholar] [CrossRef]
  90. Matsunaka, D.; Kasai, H.; Dino, W.A.; Nakanishi, H. Dynamical cluster approximation in disordered systems with magnetic impurities. J. Phys. Soc. Jpn. 2004, 73, 3448–3452. [Google Scholar] [CrossRef]

Appendix

Dynamical CPA–Locator Formalism

In the Appendix, we briefly outline an alternative but equivalent condition of the CPA, which is an extension of the CPA using the locator formalism [32]. Assuming that the spin-dependent effective medium surrounds an arbitrary site n, we consider the transfer of the carrier with spin μ between site n and the effective medium (i.e., Σ and Σ ) by the site-renormalized interactor J μ . Then, the propagator G A ( G M ) associated with the real potential of u n A ( u n M ) embedded at site n in the medium is defined by
G A = 1 ω u n A μ J μ a n μ a n μ
G M = 1 ω u n M μ J μ a n μ a n μ
When we set the coherent potential Σ μ at site n in the effective medium, the reference propagator,
P = 1 ω μ Σ μ a n μ a n μ μ J μ a n μ a n μ
is equivalent to the Green function for the effective medium. Thus, the diagonal matrix element of P is equal to F μ defined by Equation (15):
F μ ( ω ) = n μ | P | n μ = 1 ω Σ μ J μ
Equation (52) gives the relationship between J μ and F μ ; L μ 1 / ( ω Σ μ ) is called a locator. Hereafter, for the sake of simplicity, the site-diagonal elements in the Wannier representation n μ | G A | n ν are written as G μ ν A ( μ , ν = , or ↓).
Then, the spin-diagonal element of G A is given by
F μ A ( ω ) = G μ μ A = 1 ω E A J μ
and the spin-off-diagonal elements are G A = G A = 0 .
In order to obtain the explicit expression for the site-diagonal elements of G M , we rewrite Equation (50b) as
G M ( ω u n M μ J μ a n μ a n μ ) = 1
Equation (54) is written in the spin-matrix-element expression as
G M ( ω E M + I S z J ) + G M ( I S + ) = 1
G M ( I S ) + G M ( ω E M I S z J ) = 0
Equation (55) is the simultaneous equations concerning G M and G M , which can be solved after a somewhat complicated calculation using the commutation relationships between the components of S (but with no further approximations). The resulting expressions are
G M = ω E M I ( S z + 1 ) J [ ω ( E M I S z ) J ] [ ω E M I ( S z + 1 ) J ] I 2 [ S ( S + 1 ) S z 2 S z ]
G M = ω E M + I ( S z 1 ) J [ ω ( E M + I S z ) J ] [ ω E M + I ( S z 1 ) J ] I 2 [ S ( S + 1 ) S z 2 + S z ]
G M = 1 [ ω ( E M I S z ) J ] [ ω E M I ( S z + 1 ) J ] I 2 [ S ( S + 1 ) S z 2 S z ] ( I S )
G M = 1 [ ω ( E M + I S z ) J ] [ ω E M + I ( S z 1 ) J ] I 2 [ S ( S + 1 ) S z 2 + S z ] ( I S + )
Note that the site-diagonal elements of G M involve spin operators. Thus, F μ M ( ω ) is defined as the thermal average of the spin-diagonal element G μ μ M by
F μ M ( ω ) = G μ μ M = S z = S S G μ μ M ( S z ) exp ( λ S z ) / S z = S S exp ( λ S z )
where λ ( h / k B T ) is determined so as to reproduce a given value of S z [see Equation (24)]. Note that the spin-off-diagonal elements G M = G M = 0 , because G M ( G M ) includes S ( S + ) in its final form. Finally, the CPA condition in the locator formula is given by
F μ ( ω ) = ( 1 x ) F μ A ( ω ) + x F μ M ( ω )
When F μ is given, J μ is calculated by Equations (29) and (52) [i.e., J μ = ( Δ 2 / 4 ) F μ for the model band defined by Equation (27)]. Then, F μ A and F μ M are calculated by Equations (53) and (57), respectively, and consequently, F μ is again obtained by Equation (58). Therefore, F μ and J μ are determined self-consistently.
The advantage of the locator formula of the CPA is that it is straightforward to determine the species-resolved DOS, i.e., the DOS associated with each type of ion in the alloy. The local DOS at the A ( M ) site is obtained by
D μ A ( ω ) = 1 π Im F μ A ( ω )
D μ M ( ω ) = 1 π Im F μ M ( ω )
In actual calculations, we first determined F μ and Σ μ by the t-matrix formula of the CPA, then calculated J μ , and finally calculated D μ A ( ω ) and D μ M ( ω ) . We verified the relations
D μ A ( ω ) d ω = D μ M ( ω ) d ω = 1
and
D μ ( ω ) = ( 1 x ) D μ A ( ω ) + x D μ M ( ω )
which are a consequence of the CPA condition in the locator formula, Equation (58).

Share and Cite

MDPI and ACS Style

Takahashi, M. Carrier States in Ferromagnetic Semiconductors and Diluted Magnetic Semiconductors—Coherent Potential Approach—. Materials 2010, 3, 3740-3776. https://doi.org/10.3390/ma3063740

AMA Style

Takahashi M. Carrier States in Ferromagnetic Semiconductors and Diluted Magnetic Semiconductors—Coherent Potential Approach—. Materials. 2010; 3(6):3740-3776. https://doi.org/10.3390/ma3063740

Chicago/Turabian Style

Takahashi, Masao. 2010. "Carrier States in Ferromagnetic Semiconductors and Diluted Magnetic Semiconductors—Coherent Potential Approach—" Materials 3, no. 6: 3740-3776. https://doi.org/10.3390/ma3063740

APA Style

Takahashi, M. (2010). Carrier States in Ferromagnetic Semiconductors and Diluted Magnetic Semiconductors—Coherent Potential Approach—. Materials, 3(6), 3740-3776. https://doi.org/10.3390/ma3063740

Article Metrics

Back to TopTop