diluted magnetic semiconductors ( dms ) @xcite with transition temperatures as high as 150k ( gamnas at @xmath3 mn @xcite ) or even above room temperature ( gamnn , gamnp@xcite ) are attracting much attention lately in part because of possible ` spintronic ' applications @xcite .
the prototypical dms material is ga@xmath4mn@xmath5as with the mn ions substitutionally replacing ga at the cation sites .
it is widely believed that the ferromagnetism in this material is carrier induced with holes donated by mn ions mediating a ferromagnetic interaction between the randomly localized mn spins @xcite .
the coupling of carrier spins ( holes in gamnas ) and localized moments ( mn impurities ) gives rise to the unique magnetic and transport ( as well as optical ) properties in dms .
dms transport properties are influenced by the exchange interaction between the carriers and the localized moments as spin fluctuation scattering contributes to the resistivity .
the experimentally measured dc - resistivity in dms materials @xcite shows interesting behavior strongly depending on the concentration of the magnetic impurity and temperature . in in@xmath4mn@xmath5as @xcite and ga@xmath4mn@xmath5as @xcite with low mn concentration ( @xmath6 ) only an insulating behavior has been observed in transport measurements .
however , near optimal doping @xmath7 , where the highest value of @xmath2 is reported , the non - monotonic behavior ( insulator - metal - insulator ) of the resistivity as a function of temperature is observed .
a resistivity peak appears near the critical temperature ( @xmath2 ) , and the resistivity shows metallic behavior ( @xmath8 ) below @xmath2 and insulating behavior ( @xmath9 ) at higher temperatures .
the peak has been understood as the critical scattering effects at @xmath2 of spin fluctuations @xcite .    in this paper
we present theoretical calculations for dms transport properties and study the role of the carrier - spin coupling which are crucial for ferromagnetic properties .
transport measurements have proven useful in understanding the physics of the colossal magnetoresistance ( cmr ) manganites where carrier - spin coupling is also crucial @xcite . in our calculation
we use a recently developed non - perturbative method , the `` dynamical mean - field theory '' ( dmft),@xcite for calculating the dms transport properties .
a non - perturbative method is needed in dms materials because the crucial physics involves bound - state formation and other aspects of intermediate carrier - spin couplings not accessible to perturbative methods .
( we note that the most interesting phenomena in dms involve intermediate couplings and intermediate temperatures .
this regime is very difficult to treat by standard analytical or numerical methods . )
the dmft has been recently applied to the dms system to calculate the magnetic transition temperature and the optical conductivity @xcite .
dmft is essentially a lattice quantum version of the weiss mean field theory where the appropriate density of states ( including impurity band formation ) along with temporal fluctuations are incorporated within an effective local field theory .
an important ingredient of dmft @xcite is that it reduces to the standard rkky physics in the weak - coupling regime and the double exchange physics in the strong coupling regime .
our dmft results show interesting dependence of resistivity on carrier - spin coupling ( @xmath10 ) , carrier density ( @xmath11 ) , doping of the magnetic material ( @xmath12 ) , and temperature ( @xmath13 ) revealing key features of the underlying physics .
we find that our results show many similarities to colossal magnetoresistance ( cmr ) manganites @xcite , especially , for a large carrier - spin coupling and near half filling of the impurity band ( i.e. in the double exchange regime ) since the spin - disorder scattering dominates in this parameter range .
however , the dms dc - resistivity exhibits novel features not found in the cmr in the weak coupling rkky limit ( @xmath14 , where @xmath15 is the band width of the carrier ) .
these novel features arise mainly from the bound state formation of the carrier - local spins and carrier occupation of the minority spin band .
the formation of the bound state gives rise to an insulating behavior even if the carrier localization effects are not taken into account .
experimental observation of our predictions should lead to crucial information about bound state formation and impurity band physics in this problem .
even though we focus on iii - v compound based dms such as ga@xmath4mn@xmath5as , the results presented in this paper are general for all dms materials .
an important limitation of dmft is that it is a non - perturbative _ local _ theory that can not really incorporate _ spatial _ fluctuations , and therefore resistive charged impurity scattering with its strong momentum dependence is essentially impossible to handle in dmft .
we therefore consider the boltzmann transport theory to calculate the resistivity of the dms systems @xcite .
the boltzmann theory is used to calculate the charged impurity disorder limited dms resistivity ( whereas the dmft is used for obtaining the spin disorder limited dms resistivity ) .
charged impurity disorder arises here both from the ionized mn acceptor in gamnas as well as other charged defects / impurities invariably present in a semiconductor . using relaxation time approximation
we assume that the boltzmann resistivity is due to ionized impurities .
the carriers ( holes ) are scattered by the screened coulomb potential , which we calculate using the linearized thomas - fermi ( tf ) approximation and the random phase approximation ( rpa ) .
we find that the dominant temperature dependence of the resistivity comes from the change in the screening length in the high density metallic samples . in the metallic gamnas
dms samples the change in the measured resistivity , when the temperature goes from @xmath2 to zero , is about 20% in good agreement with our calculation . in the ferromagnetic state ( @xmath1 ) , as the temperature decreases from @xmath16 , we find strong temperature dependence of the resistivity arising from the low temperature screening function and the thermal change of the carrier densities in each spin - split subband . when all the carriers are polarized ( this happens in the low density limit ) the screening function is almost independent of temperature and the resistivity is also temperature independent .
inclusion of both spin disorder and ionized impurity disorder in dms transport is the important ingredient of our theory .
this paper is organized as follows . in sec .
ii , we describe our model and our theoretical approach based on the dynamical mean - field theory ( dmft ) . in sec .
iii , we study , in detail , the density of states and describe the formation of the spin polarized impurity band . in sec .
iv , we provide the results of calculated dms resistivity within dmft . in sec .
v we calculate the transport resistivity within the boltzmann transport theory for ionized impurities . in sec .
vi , we summarize our qualitative findings and providing a critical discussion of the applicability of our results to dms systems . a brief conclusion is given in sec .
our basic model of dms systems is that of magnetic dopants ( `` impurities '' ) interacting through a local exchange coupling with carriers in the host semiconductor material .
the generally accepted hamiltonian of the system is given by @xmath17 where @xmath18 describes carrier propagation in the host semiconductor band .
for simplicity we consider a host material with a single non - degenerate band .
we therefore write @xmath19 \psi _ { \alpha } ( x),\ ] ] where @xmath20 is the spin index and @xmath21 is a random potential arising from non - magnetic defects in the material ( e.g. , as antisite defects , unintentional background charged impurities , etc . ) .
the second ( magnetic ) term in eq . ( [ ham ] ) , @xmath22 , describes coupling of the carriers to an array of impurity ( e.g. mn ) spins at positions @xmath23 , @xmath24 \psi_{\beta } ( r_i),\ ] ] where @xmath10 is the local exchange coupling between the spin of the magnetic impurity and the the spins of the semiconductor carriers , @xmath25 is the ( coulombic ) potential arising from the magnetic dopant , @xmath23 are the positions of the magnetic dopants , and @xmath26 is the pauli matrix . here
we absorb the magnitude of the impurity spin into the coupling @xmath10 ( which we take to be positive ) , and represent the spin direction by the unit vector @xmath27 .
the third term in eq .
( [ ham ] ) , @xmath28 , is the direct mn@xmath29mn short - range antiferromagnetic exchange interaction @xmath30 where @xmath31 is a direct antiferromagnetic exchange coupling between impurity spins .    in this section
we approximate our model by neglecting the nonmagnetic random potential , @xmath21 , and the direct antiferromagnetic exchange interaction , @xmath31 .
lattice defects may be playing an important role in determining magnetic and transport properties of the samples , but we assume here that these defects enter our theory only in determining the basic parameters of the model , namely , the density of magnetically active dopants @xmath32 , the hole density @xmath33 , and perhaps the local effective exchange coupling @xmath10 between the holes and the magnetic impurities , and do not include any defects into our model explicitly .
we include charged impurity scattering through the boltzmann equation in section iv of this paper .
we believe that the effects of the antiferromagnetic coupling between magnetic impurities are either negligibly small or incorporated into the effective parameters of the model . actually ,
in the parameter range of interest to us ( @xmath34 ) , where dms ferromagnetism typically occurs , the magnetic impurities are separated from each other by non - magnetic atoms , and this short - range antiferromagnetic interaction , which rapidly decays with the distance , should be negligible .
these approximations are nonessential and are done in the spirit of identifying the minimal dms magnetic model of interest .
both of these effects , which may be of quantitative importance in some situations , can be included in the theory by adjusting the parameters of the model or perhaps at the cost of introducing more unknown parameters characterizing these interactions .
recently several theories for dms ferromagnetism explicitly including spatial disorder effects and antiferromagnetic coupling have been developed @xcite .    in our dmft dms model
there are two sources of coupling between the carrier and the impurity magnetic moment : a spin - spin coupling ( @xmath10 ) and a potential scattering ( @xmath25 ) .
the crucial physical issues are revealed by the consideration of a ferromagnetic state in which all impurity spins @xmath35 are aligned , say , in the @xmath36 direction .
then the carriers with spin parallel to @xmath35 feel a potential @xmath37 on each magnetic impurity site and anti - parallel carriers feel a potential @xmath38 .
these potentials self - consistently rearrange the electronic structure .
the spin - dependent part of this rearrangement provides the energy gain which stabilizes the ferromagnetic state .
the key physics issue is , evidently , whether the potential @xmath39 is weak ( so its effect on carriers near the lower band edge is simply a scattering phase shift ) or strong ( so only majority spin or perhaps both species of carriers are confined into spin - polarized impurity bands ) .
recent density functional supercell calculations @xcite suggest that in gamnas @xmath37 is close to the critical value for bound state formation for the majority spin systems .
unfortunately , the precise effective values of @xmath10 and @xmath25 are typically unknown in a dms system , and may have to be extracted experimentally .
we assume that magnetic impurities under consideration enter substitutionally at the cation sites ( _ e .
g. _ mn impurities at ga sites ) and the iii@xmath4mn@xmath5v system as a lattice of sites , which are randomly nonmagnetic ( with probability @xmath40 ) or magnetic ( with probability @xmath12 ) , where @xmath12 now indicates the relative concentration ( _ i.e. _ per ga site ) of active mn local moments in iiimnv .
if more complete information about mn locations on the gaas lattice becomes available it will be straightforward to incorporate that in the dmft formalism .
we now introduce the dmft for the hamiltonian given by eq .
( [ ham ] ) . within the general scheme of the dmft
, the local ( momentum independent ) self energy of the system , @xmath41 , can be obtained from the time dependent mean field function .
then the single particle green function is approximated by @xmath42 where @xmath43 is the chemical potential . with the local self energy
all of the relevant physics may be determined from the local ( momentum - integrated ) green function defined by @xmath44 where @xmath45 is the density of states ( dos ) for the noninteracting system .
the information of the lattice geometry is included through the noninteracting dos . in our model
@xmath46 is a matrix in spin index and depends on whether one is considering a magnetic ( @xmath47 or non - magnetic ( @xmath48 site . since @xmath46 is a local function , it is the solution of a local problem specified by a mean - field function @xmath49 , which is related to the partition function @xmath50 with action @xmath51   c_{\beta } ( \tau ^{\prime } ) \nonumber \\ & + & \int d\tau \sum_{\alpha,\beta } c_{\alpha } ^{+}(\tau ) \left [ j \hat { \textbf{s } } \cdot \bm{\sigma}_{\alpha   \beta }   + w \delta_{\alpha,\beta } \right ] c_{\beta } ( \tau)\end{aligned}\ ] ] on the @xmath52 ( magnetic ) site and @xmath53 on the non - magnetic ( @xmath54 ) site . here
@xmath55 ( @xmath56 ) is the destruction ( creation ) operator of a fermion in the spin state @xmath20 and at time @xmath57 .
@xmath58 plays the role of the weiss mean field ( bare green function for the local effective action @xmath59 ) and is a function of time .
@xmath60 depends only on frequency and is therefore the solution of a single - site problem .
the local green function @xmath60 of the effective single - site problem is solely determined by the partition function , @xmath61 , namely , @xmath62 which is identical to the local green function computed by performing the momentum integral using the same self energy .
then the self - energy is defined by @xmath63    the @xmath52-site mean - field function @xmath64 can be written as @xmath65 with @xmath66 the magnetization direction and @xmath67 vanishing in the paramagnetic state ( @xmath0 ) . since the spin axis is chosen parallel to @xmath68 @xmath69
becomes a diagonal matrix with components parallel ( @xmath70 ) and antiparallel ( @xmath71 ) to @xmath68 .
the form of the dispersion given in full hamiltonian eq .
( 1 ) applies only near the band edges .
it is necessary for the method to impose a momentum cutoff , arising physically from the carrier band - width .
we impose the cutoff by assuming a semicircular density of states @xmath72 with @xmath73 .
the parameter @xmath15 is chosen to correctly reproduce the band edge density of states .
other choices of upper cutoff would lead to numerically similar results .
this choice of cutoff corresponds to a bethe lattice in infinite dimensions .
other ( perhaps more realistic ) choices for the density of states would give results qualitatively similar to our results since the band edge density of states has the correct physical behavior in our model . for this @xmath74 the self consistent equation for @xmath75
obeys the equation @xmath76^{-1 } \right \rangle,\end{aligned}\ ] ] where the angular brackets denote averages performed in the ensemble defined by the appropriate @xmath61 , i.e. , @xmath77 , with @xmath78 . with these mean field functions the local green function can be written as @xmath79 the self energies are evaluated using eq .
( [ self ] ) and the full green function from eq .
( [ green ] ) .
physical observables can be obtained from the full green function @xmath80 . in particular , the mean field function @xmath81 can be easily calculated at @xmath82 @xmath83 and @xmath84 @xmath85
the density of states plays crucial roles in determining the physical properties of the dms system .
especially , the formation of the impurity band arising from the impurity doping gives rise to many different aspects from the continuum ( i.e. virtual crystal approximation ) semiconductor band model . in this section
, we calculate the dos for different parameters ( @xmath10 , @xmath12 , @xmath25 , and @xmath13 ) and show how the impurity bands are formed and separated from the main band .
we describe the dos of dynamical mean - field calculations applying to simple semi - circle models .
the dos is given by the imaginary part of the green function @xmath86    in our model for a strong magnetic coupling @xmath87 two spin polarized impurity bands appear at the bottom ( majority spin ) and at the top ( minority spin ) of the main band .
each isolated impurity band has the weight @xmath12 .
however , if the coupling is not strong ( @xmath88 ) the impurity band is not completely separated from the main bend .
all dms samples show that the carrier density is much smaller than the impurity concentration ( @xmath89 ) due to the heavy compensation .
as the system is the partially compensated , the chemical potential @xmath43 is located in the lower impurity band ( if the impurity bands are formed ) or the lower band edge ( if the bands are not formed ) .
thus all physical properties are determined in the lower energy band edge . throughout this paper
we only show the dos near the lower energy band .
temperature dependence of the dos for @xmath90 and for a fixed coupling ( a ) @xmath91 and ( b ) @xmath92 .
the evolutions of majority ( minority ) spin dos show in top ( bottom ) panels for various temperature @xmath930.0 , 0.5 , 0.7 , 0.9 , and 1.0 . , title="fig:",width=249 ]   temperature dependence of the dos for @xmath90 and for a fixed coupling ( a ) @xmath91 and ( b ) @xmath92 .
the evolutions of majority ( minority ) spin dos show in top ( bottom ) panels for various temperature @xmath930.0 , 0.5 , 0.7 , 0.9 , and 1.0 .
, title="fig:",width=249 ]     the calculated majority spin dos at @xmath82 ( upper half ) and @xmath16 ( lower half ) is shown for various doping @xmath94 0.001 , 0.01 , 0.05 , 0.1 and a fixed coupling ( a ) @xmath91 , and ( b ) @xmath95 . , title="fig:",width=249 ]   the calculated majority spin dos at @xmath82 ( upper half ) and @xmath16 ( lower half ) is shown for various doping @xmath94 0.001 , 0.01 , 0.05 , 0.1 and a fixed coupling ( a ) @xmath91 , and ( b ) @xmath95 .
, title="fig:",width=249 ]    in fig .
[ fig_1 ] we show the calculated dos for various temperatures as a function of energy .
the evolutions of majority ( minority ) spin dos are shown in top ( bottom ) panels . in fig .
[ fig_1](a ) the strength of the coupling constant ( @xmath91 ) is not strong enough to form the impurity band .
note that this value of coupling constant ( @xmath91 ) is the critical value for impurity band formation as @xmath96 . at @xmath82 the majority ( minority ) spin band
is shifted to lower ( higher ) energy compared with the noninteracting band which has a band edge at @xmath97 .
thus all carriers are fully polarized when the carrier density is low .
as the carrier density increases they start to occupy the minority band if the chemical potential crosses into the minority spin band .
recently we showed that the optical conductivity of the system is dramatically changed with the occupation of the minority band @xcite .
we show in this paper that the calculated dc transport properties are also very sensitive to the minority band occupation . as the temperature increases
, the minority band occupation grows and the carriers with minority spin increase due to thermal fluctuations . at @xmath16 both spins are equally populated and the bands becomes symmetric . as expected we have a separated impurity band from main band for a strong coupling @xmath92 shown in fig . [
fig_1](b ) .
when @xmath98 the impurity band is fully occupied , no low energy hopping processes to main band are allowed and the system becomes a band insulator . if the impurity band is partially occupied the delocalization energy increases . at the half filling of the impurity band
the system has the highest @xmath2 @xcite .
as the temperature increases spin disorder grows and the band becomes symmetric .
but the impurity band width of the paramagnetic state ( @xmath99 ) is smaller than that of the ferromagnetic state .
this band shrinking occurs because the neighboring spins in the paramagnetic state are uncorrelated .     the calculated majority spin dos at ( a ) @xmath82 and ( b ) @xmath16 is shown for a fixed @xmath90 and for various coupling constant @xmath100 0.5 , 1.0 , 1.25 , 1.5 , 2.0 .
, width=268 ]    in fig .
[ fig_2 ] we show the calculated majority spin dos corresponding to the disordered spin state ( at @xmath16 , bottom panels ) and ferromagnetic state ( at @xmath82 , top panels ) .
the evolutions of the energy ( @xmath101 ) dependent dos are shown for different doping parameter @xmath940.001 , 0.01 , 0.05 , 0.1 and for fixed coupling constant ( a ) @xmath91 and ( b ) @xmath95 . in our model
the impurity level ( acceptor energy level ) and the formation of an impurity band depend on the ferromagnetic coupling @xmath10 . if @xmath102 the impurity level is not isolated from the main band , but if @xmath103 we find an isolated impurity level below the main band .
the small dot in fig . [ fig_2](b ) indicates the isolated impurity level in the dilute limit ( @xmath96 ) .
as @xmath12 increases for @xmath103 an impurity band centered around the impurity level is formed below the main band .
for @xmath88 the impurities give rise to band tailing in the main band edge instead of forming impurity band .
the band width of the impurity band for @xmath103 increases with @xmath12 since the number of states of the impurity band increases with @xmath12 . if @xmath12 is bigger than @xmath104 the impurity band merges into the main band .
for @xmath95 we have @xmath105 at @xmath82k and @xmath106 at @xmath16 .    in fig .
[ fig_3 ] we show the majority spin dos at ( a ) @xmath82 and ( b ) @xmath16 for a fixed @xmath90 and for various coupling constant @xmath100 @xmath107 , @xmath108 , @xmath109 , 1.5 t , 2.0 t .
for @xmath88 we see the expected band shift proportional to @xmath110 . for @xmath111 an impurity band centered at impurity level and containing @xmath12 states
is seen to split off from the main band . due to the widening of the impurity band with @xmath12 we find the separated impurity band when @xmath112 at @xmath82 and @xmath113 at @xmath16 for @xmath90 .
thus the separated impurity band is expected with small @xmath12 and large values of @xmath10 .
the majority ( top panel ) and minority ( bottom panel ) dos at @xmath82 are shown for a fixed value of @xmath95 and @xmath90 for various potential scattering @xmath114 , @xmath115 , 0.0 , 0.25 , 0.5 . ,
width=268 ]    in fig .
[ fig_4 ] we show the calculated dos at @xmath82 for a fixed value of @xmath95 and @xmath90 , and for various values of potential scattering @xmath116@xmath117 , @xmath118 , 0.0 , @xmath119 , @xmath107 .
the dos of majority ( minority ) spin is shown in the top ( bottom ) panel .
when we include the potential scattering ( @xmath25 ) in addition to the spin - spin coupling ( @xmath10 ) , the carriers with spin parallel to local impurity spin feel a potential @xmath37 on each magnetic impurity site and anti - parallel carriers feel a potential @xmath38 .
thus the formation of the impurity band and the corresponding physical properties of the system depend on the combined coupling @xmath120 .
[ fig_4 ] shows that while the band edge of the minority spins is slightly dependent on the potential scattering , the majority spin dos is strongly affected by the potential scattering .
even weak potential scattering can change the extended majority spin band into the spin - polarized impurity bands or the well formed impurity band into the band tail of the main band .    in the following section
we show that our calculated resistivities depend strongly on whether the carriers are within the impurity band or in the band tail of the extended main band .
the conductivity is calculated from the usual kubo formula .
the kubo formula for the conductivity @xmath26 involves the two particle current - current response function .
since the irreducible vertex in the response function is purely local in our approximation of dmft there is no vertex correction @xcite .
thus , within this approximation only the simple bubble survives and the real part of the finite frequency conductivity is given by @xmath121 } { \omega } \nonumber \\ & \times &   a_{\sigma } ( \varepsilon , \omega ) a_{\sigma } ( \varepsilon , \omega + \omega ) , \label{sigxx}\end{aligned}\ ] ] where @xmath122 is the spectral function , @xmath123 is the current vertex for the bethe lattice @xcite , @xmath124 is the fermi distribution function . for
the hypercubic lattice @xcite @xmath125 has been used in eq .
( [ sigxx ] ) , but for the bethe lattice the explicit form is derived in ref . .
the dc - resistivity @xmath126 is then found from eq .
( [ sigxx ] ) in the limit @xmath127 with @xmath128 , which is given by @xmath129     the density dependent resistivity at @xmath16 for a fixed value of @xmath90 and for various coupling constant @xmath130 , 1.0 , 1.2 , 1.4 , 1.5 , 2.0 .
insets show the dos with the corresponding parameters .
, width=288 ]    in fig .
[ fig_5 ] we show the calculated resistivity at @xmath16 as a function of carrier density for a fixed value of @xmath90 and various coupling constants @xmath131=0.7 , 1.0 , 1.2 , 1.4 , 1.5 and 2.0 .
insets show the density of states near the band edge corresponding to the disordered spin state ( @xmath16 ) .
all calculated results show that @xmath132 diverges as @xmath133 due to the absence of carriers . as density increases
we find two different behaviors depending on the formation of the impurity band .
( the critical coupling constant which gives rise to the formation of the well separated impurity band below the main band for @xmath90 , @xmath16 is @xmath134 . ) in the weak coupling limit @xmath135 , where the impurity band formation dose not happen , the resistivity decreases as the density increases monotonically .
however , in the strong coupling limit @xmath136 , where the impurity band is formed , the resistivity diverges again when the impurity bands are fully filled ( i.e. , for @xmath98 ) and the system becomes a band insulator . in the intermediate coupling regime ( @xmath137 )
the resistivity shows non - monotonic behavior .     the calculated resistivity for ( a ) a strong coupling @xmath92 with densities @xmath138 0.01 , 0.015 , 0.02 , 0.025 , 0.03 0.035 , 0.04 , 0.045 ( top to bottom)and ( b ) a weak coupling @xmath92 with densities @xmath138 0.008 , 0.01 , 0.015 , 0.02 , 0.025 , 0.035 , 0.04 , 0.045 ( top to bottom ) .
dots indicate @xmath2 for given densities . , title="fig:",width=240 ]   the calculated resistivity for ( a ) a strong coupling @xmath92 with densities @xmath138 0.01 , 0.015 , 0.02 , 0.025 , 0.03 0.035 , 0.04 , 0.045 ( top to bottom)and ( b ) a weak coupling @xmath92 with densities @xmath138 0.008 , 0.01 , 0.015 , 0.02 , 0.025 , 0.035 , 0.04 , 0.045 ( top to bottom ) .
dots indicate @xmath2 for given densities . , title="fig:",width=240 ]    in fig .
[ fig_6 ] we show the calculated resistivity as a function of temperature for various density . for the strong coupling limit ( @xmath92 ) , in fig .
[ fig_6](a ) , we find a crossover ( resistivity peak ) separating a good metal at low @xmath13 from a semiconductor at higher @xmath13 .
the resistivity peaks are proportional to the energy separation between the chemical potential and the band edge of the minority band . in the high temperature regime above the resistivity peak the decreasing resistance with increasing temperature , characteristic of a semiconductor or insulator , is due to thermal excitation of the carriers from impurity band to the upper minority spin band
. the resistivity decreases with density because the carriers in the main band are scattering off the impurities .
but at low temperatures most carriers in the impurity band contribute to the scattering and give non - monotonic density dependence of the resistivity .
the metallic behavior at low temperature can be understood by the disappearance of the coherent central quasiparticle peak in the dos . for the weak coupling limit ( @xmath91 ) , in fig .
[ fig_6](b ) , the resistivity can be explained by the scattering effects in the main band except the behavior in the ferromagnetic state ( @xmath1 ) . for @xmath1
the resistivity changes from insulating at low densities to metallic at high density .
details of this behavior are given in fig .
[ fig_7 ] .
we see that the resistivity crossover takes place only at low densities since the minority band is occupied by the carriers at high density .
the low temperature dc - resistivity as a function of temperature for a fixed @xmath90 and for various densities . in ( a ) the results for a strong coupling constant ( @xmath92 ) and densities @xmath139 , 0.012 , 0.015 , 0.025 , 0.3 , 0.35 , 0.38
, 0.40 are shown . in ( b )
we use @xmath91 and @xmath140 , 0.016 , 0.0163 , 0.0165 , 0.18 , 0.2 , 0.25 , 0.3 ( from top to bottom ) .
dots indicate the @xmath2 for given parameters .
, title="fig:",width=240 ]   the low temperature dc - resistivity as a function of temperature for a fixed @xmath90 and for various densities . in ( a ) the results for a strong coupling
constant ( @xmath92 ) and densities @xmath139 , 0.012 , 0.015 , 0.025 , 0.3 , 0.35 , 0.38 , 0.40 are shown . in ( b )
we use @xmath91 and @xmath140 , 0.016 , 0.0163 , 0.0165 , 0.18 , 0.2 , 0.25 , 0.3 ( from top to bottom ) .
dots indicate the @xmath2 for given parameters .
, title="fig:",width=240 ]    in fig .
[ fig_7 ] we show the dc - resistivity as a function of temperature for different densities in the low temperature regime ( @xmath141 ) . in fig .
[ fig_7](a ) we use the parameters @xmath90 and a strong coupling @xmath92 . in this case
all carriers ( if @xmath142 ) occupy the impurity band and stay mainly at mn sites . thus , the carriers in the impurity band follow the fluctuation of the localized mn spin . at half filling of the spin polarized impurity band ( @xmath143 ) the resistivity has the lowest value and
its behavior corresponds to that of the double exchange ( de ) model ( see fig .
[ fig_10 ] ) , that is , the resistivity decreases below the critical temperature because of the spin disorder scattering as the temperature decreases @xcite . in the low density limit
( @xmath144 , less than half filling ) the resistivity is dominated by the ` impurity band ' contribution and the resistivity increases as density decreases due to the lack of mobile carriers . as the carrier density is increased above the half filling ( @xmath145 ) the resistivity increases due to the filling of the band .
we also find a very different temperature behavior of the resistivity from the low density case .
as the temperature decreases the resistivity increases just below @xmath2 , then decreases at very low temperatures .
the counter - intuitive increase of resistivity just below @xmath2 ( as @xmath13 is decreased ) arises because , as the carrier spins are aligned to the impurity spins , the binding of the carriers to the impurity spins increases corresponding to an increase in the basic scattering rate . in the very low temperature range , however , the de - like mechanism dominates , which gives rise to decrease of the resistivity .
( when the impurity band is spin - polarized , carriers which are bound to impurity site must have spins parallel to impurity spin .
thus , as the spins order ferromagnetically , the basic ability of carriers to move in the impurity band increases . )
the overall temperature dependence of resistivity is very weak in the strong coupling limit .
this weak @xmath13-dependence of the dc resistivity below @xmath2 occurs because the increase in scattering rate due to the binding is compensated by de - like mechanism .    in fig .
[ fig_7](b ) we show the results for a weak coupling limit , @xmath91 . in this case
the impurities contribute to form the band tail of the main band and the chemical potential lies in the main band . in the low density limit ( @xmath146 , where @xmath147 is the density above which the minority band starts to fill ) the resistivity below @xmath2 increases as the temperature decreases .
the increase in resistivity is due to increased carrier - spin coupling as mentioned above .
but , in the high density limit ( @xmath148 , where at @xmath82 the minority spin - band is occupied ) we find the resistivity decreasing as the temperature decreases .
this metallic behavior in high density regime and for a low coupling constant can be understood by the small scattering rate of the carriers in the minority band . if the carriers are near the edge of the majority spin band , the carriers form a spin - polarized bound state , so the effective scattering rate strongly increases , which gives rise to the insulating behavior as shown in the low density results . on the other hand , if the carriers are in the minority spin band the carriers form an anti - bound state at the top of the band , so that at the physically relevant lower band edge , the effective scattering rate decreases .
in addition in three dimensions the vanishing of the density of states at the band edge further decreases the scattering .
these effects are quite large and dominate as the minority band is occupied , which gives rise to the metallic behavior in the high density limits . in the paramagnetic state
the resistivity is almost temperature independent . in the limit
@xmath0 the local green function @xmath149 is temperature independent and the resistivity depends on temperature only weakly through the fermi function ( thermal smearing around the chemical potential ) .
the energy evolution of optical conductivity for ( a ) a intermediate coupling @xmath91 and ( b ) a strong coupling @xmath92 . solid ( dashed ) lines correspond to the results for @xmath82 ( @xmath16 ) .
insets show the dc conductivities with corresponding parameters .
arrows indicate @xmath2 . , width=288 ]    in fig .
[ fig_8 ] we show the relation between the optical conductivity and dc - resistivity .
the main panels of fig .
[ fig_8 ] show the evolution of the conductivity for two couplings ; weak ( @xmath91 , fig .
[ fig_8](a ) , where the impurity band formation is not accomplished ) , and strong ( @xmath92 , fig .
[ fig_8](b ) , where the impurity band is well formed ) ; the insets show the dc resistivity with the same parameters .
the solid ( dashed ) lines indicate the optical conductivity for @xmath82 ( @xmath16 ) .
when the impurity band is not formed ( for @xmath91 ) we find approximately the drude form for optical conductivity expected for carriers scattering off random impurities ( a closer examination reveals minor differences due to density of states variations near the band edge ) . in this case
( and for @xmath150 ) the dc conductivity shows insulating behavior due to the formation of the bound state .
since the carrier density is low enough not to fill the minority band , the formation of anti - bound state dose not take place , which reduce the scattering rate . in the @xmath92 case
the density of states plot shows the formation of an impurity band and the corresponding conductivity has two structures : a low - frequency quasi - drude peak corresponding to motion within the impurity band and a higher frequency peak corresponding to excitations from the impurity band to the main band . in this case
the dc resistivity shows metallic behavior because the reduction of spin - disorder scattering dominates over the bound - state formation . in in@xmath4mn@xmath5as
@xcite we find the drude - like conductivity is correlated with the insulating behavior , but in ga@xmath4mn@xmath5as @xcite the mid infrared peak in the optical conductivity is closely related to the decrease of the dc - resistivity below the critical temperature .
the calculated dc - resistivity as a function of temperature for a fixed @xmath90 , @xmath139 , @xmath91 , and for various potential scattering @xmath114 , @xmath151 , 0.0 , 0.05 , 0.2 .
arrows indicate the critical temperature @xmath2 .
, width=288 ]    fig .
[ fig_9 ] shows the sensitivity of the predicted behavior to potential scattering . in this figure
we use the parameters : @xmath90 , @xmath139 , @xmath91 , and for various potential scatterings @xmath152@xmath153 , @xmath151 , 0.0 , 0.05 , 0.2 . at zero scalar potential ( @xmath154 , middle panel ) the impurity band
is not formed , and the resistivity shows insulating behavior below @xmath2 due to the bound state formation of the carrier spins with impurity spins . as the potential is made more attractive ( negative ) , the impurity band features become pronounced and the spin - disorder scattering decreases ( the carriers move easily in the impurity band ) . as the potential is made more repulsive ( positive ) , the impurity band rapidly rejoins the main band .
this reduces the energy gap between the fermi energy and the band edge of the minority band , and for large enough repulsive potential the minority band starts to be filled by the carriers .
thus , the anti - bound state formation dominates and it reduces the scattering rate below @xmath2 and shows a metallic behavior .
the dms transport properties are sensitive to the combined coupling @xmath120 , and not solely to the exchange coupling @xmath10 .
resistivity @xmath155 for various magnetic field is given as a function of temperature for ( a ) @xmath156 , @xmath157 , @xmath158 , and ( b ) @xmath156 , @xmath90 , @xmath143 .
insets show the normalized magnetization @xmath159 , where @xmath160 is the saturation magnetization .
, title="fig:",width=211 ]   resistivity @xmath155 for various magnetic field is given as a function of temperature for ( a ) @xmath156 , @xmath157 , @xmath158 , and ( b ) @xmath156 , @xmath90 , @xmath143 .
insets show the normalized magnetization @xmath159 , where @xmath160 is the saturation magnetization .
, title="fig:",width=211 ]    now we compare the transport properties of dms to those of another system with strong carrier - spin couplings , namely cmr manganites @xcite . in cmr , instead of being dilute random impurities as in dms ,
the mn ions form an ordered lattice .
they possess large local moment , to which mobile carriers are very strongly coupled .
thus instead of a spin - polarized impurity band , there is a spin - polarized conduction band , sufficiently well separated from other spin bands .
the periodic arrangement of the mn sites means that ( in the absence of other physics ) the scattering rate decreases as t is lowered .
cmr materials can be understood well by the double exchange ( de ) model @xcite . in our model
this corresponds to @xmath161 , in which all magnetic ions replace ga at the cation sites . in this case
( @xmath161 ) the fully polarized spin band is well separated from the other bands instead of forming impurity band .
the temperature dependences of the resistivity @xmath155 and magnetization @xmath162 are given in fig .
[ fig_10](a ) for @xmath157 with @xmath156 , @xmath158 , and various magnetic fields .
above @xmath2 the resistivity has a small temperature dependence since the local spin fluctuation is saturated above @xmath2 . below @xmath2 resistivity decrease as magnetization increases .
the origin of the resistivity dependence on the magnetization is spin disorder scattering , which gives rise to the scaled behavior of the resistivity , @xmath163 , where @xmath164 is a temperature / field independent constant @xcite .
the origin of the resistivity is qualitatively explained by the carrier scattering due to the thermally fluctuating spin configurations , or the spin - disorder scattering . as the spontaneous or the induced magnetic moment
is developed , the amplitude of the spin fluctuation decreases so that the resistivity also decreases . in fig .
[ fig_10](b ) we show the resistivity for dms system with @xmath90 .
the overall features of temperature and magnetic field dependence look similar to the de model .
however , we find that the negative magnetoresistence at @xmath2 is very weak and the resistivity above @xmath2 is not saturated
. the fast drop of the resistivity just below @xmath2 can be explained by spin - disorder scattering .
these ideas of de model have limited applicability to the dms systems , i.e. only for strong couplings and near half filling of the spin polarized impurity band .
as shown in previous figures these ideas can not explain the resistivity behavior in low coupling and high density regimes .
note that an important ingredient of dms transport properties is missing from our dmft theory which , while accounting well for the non - perturbative effects of spin disorder ( @xmath10 ) and local potential ( @xmath25 ) scattering by the magnetic impurities , leaves out all ionized impurity disorder that may very well be important .
when the dominant scattering mechanism is the scattering by charged impurities the boltzmann transport theory may be used to calculate the electrical resistivity of the carriers since dmft is not well - suited for treating long - range disorder . due to the band splitting in the ferromagnetic state the carrier densities @xmath165 for spin up / down are not equal .
note that the total density @xmath166 .
in this situation the total conductivity can be expressed as a sum of spin up / down contributions @xmath167 where @xmath168 is the conductivity of the ( @xmath169 ) spin subband .
the conductivities @xmath168 are given by @xmath170 where @xmath171 is the carrier effective mass , and the energy averaged transport relaxation time @xmath172 for the ( @xmath169 ) subbands are given in the boltzmann theory by @xmath173 }     { \int d\varepsilon \varepsilon \left [       -\frac{\partial f^{\pm}(\varepsilon)}{\partial \varepsilon } \right    ] } , \label{tau}\ ] ] where @xmath174 is the energy dependent relaxation time for the ( @xmath169 ) subbands , and @xmath175 is the carrier ( fermi ) distribution function @xmath176/k_bt}},\ ] ] where @xmath177 is the chemical potential at finite temperature .
the energy dependent relaxation time is given in the born approximation by @xmath178^{-1 } = \frac{2\pi}{\hbar } \int \frac{d^3k'}{(2\pi)^3}n_i|v^{sc}_{\bf k - k'}|^2 ( 1-\cos\theta_{\bf kk ' } ) \delta(\varepsilon_{\bf k } - \varepsilon_{\bf k'}),\ ] ] where @xmath179 is the charged impurity concentration , and @xmath180 is the screened carrier - impurity coulomb interaction , which can be expressed as @xmath181 ^ 2},\ ] ] where @xmath182 is the charge of impurities , @xmath183 the background ( gaas ) lattice dielectric constant , and @xmath184 the temperature dependent screening function .
wave vector dependence of the effective screening function ( polarizability function @xmath185 normalized its paramagnetic state value @xmath186 ) for various polarizations @xmath187 . , width=211 ]    in this paper we consider two screening approaches : thomas - fermi ( tf ) approximation and random phase approximation ( rpa ) . within tf screening
we have @xmath188,\ ] ] where @xmath189 with @xmath190 being the fermi wave vector of the spin unpolarized state and @xmath191 the effective bohr radius .
@xmath192 is the fermi wave vector of the each spin - split subband .
note that tf screening wave vector , being a long wavelength approximation , is dependent only on the spin polarization , but not on the temperature explicitly .
it is therefore temperature independent above the critical temperature . within tf
screening approximation we have the energy dependent relaxation time ( by integrating eq .
21 ) @xmath178^{-1 } = \frac{2\pi}{3}\frac{n_i}{n}\frac{e_{f}^0}{\hbar } \frac{k k_f^0 } { \tilde{k}_f^2 } f(4k^2/q_{s}^2),\ ] ] where @xmath193 and @xmath194 $ ] . when a system is fully polarized we have @xmath195 and @xmath196 . for the random phase approximation ( rpa )
the temperature dependent screen function is given by @xmath197,\ ] ] where @xmath198 is the temperature dependent static lindhard function for each spin subband . in fig .
[ fig_11 ] we show the screening function @xmath199 , where @xmath186 , for various polarization @xmath187 . for an unpolarized state ( @xmath200 ) there is an inflection point at @xmath201 , which is the usual kohn anomaly . for partially polarized state @xmath202 we have two inflection points at @xmath203 and the value of @xmath204 decreases as @xmath205 increases .
when the system is fully polarized @xmath206 , and @xmath185 has an inflection point at @xmath207 , and @xmath208 .
calculated resistivity as a function of @xmath209 for various value of @xmath210 1.0 , 2.0 , 5.0 , 10.0 ( top to bottom ) . solid ( dashed ) lines show the results calculated with rpa ( tf ) screening function.,width=240 ]    at @xmath82 , @xmath211 where @xmath212 is the fermi energy for the ( @xmath169 ) subbands , and then @xmath213 giving the familiar result @xmath214 .
within tf screening the ratio of the resistivity of the fully polarized state to that of the unpolarized state becomes @xmath215 , where @xmath216 .
( note @xmath217 . ) in fig .
[ fig_12 ] we show the calculated resistivity as a function of @xmath209 for several value of @xmath218 .
solid ( dashed ) lines show the results calculated with rpa ( tf ) screening . in general , for small values of @xmath218 , @xmath219 , but for large @xmath218 , @xmath220 . at @xmath221
we have @xmath222 . for
gaas @xmath223 corresponds to the hole density @xmath224 . in relatively high density limits ( large @xmath218 ) the two approximations agree very well , which indicates that the @xmath225 scattering mostly contributes to the scattering time .
however , in the low density regime ( small @xmath218 ) we that find the two screening theories give very different results for the spin polarized state .
noting that the tf approximation is just the long wavelength limit of rpa , we emphasize that rpa is obviously the better theory .     calculated resistivity , @xmath226 , as a function of temperature . in ( a )
the rpa results are shown with hole densities @xmath227 0.05 , 0.1 , 0.2 , 0.5 . in
( b ) the tf screening results are shown with hole densities @xmath227 0.01 , 0.05 , 0.1 , 0.2 . here
the impurity density @xmath228 is used .
, title="fig:",width=211 ]   calculated resistivity , @xmath226 , as a function of temperature . in ( a ) the rpa results are shown with hole densities @xmath227 0.05 , 0.1 , 0.2 , 0.5 . in
( b ) the tf screening results are shown with hole densities @xmath227 0.01 , 0.05 , 0.1 , 0.2 . here
the impurity density @xmath228 is used .
, title="fig:",width=211 ]     calculated resistivity as a function of temperature for a hole density ( a ) @xmath229 0.1 and ( b ) @xmath2300.2 with finite temperature rpa screening function .
the top insets show the temperature dependence of the mn ions and hole magnetization .
the bottom insets show the calculated temperature dependence of the tf screening wave vector @xmath231 . here
the impurity density @xmath228 is used .
, title="fig:",width=211 ]   calculated resistivity as a function of temperature for a hole density ( a ) @xmath229 0.1 and ( b ) @xmath2300.2 with finite temperature rpa screening function .
the top insets show the temperature dependence of the mn ions and hole magnetization .
the bottom insets show the calculated temperature dependence of the tf screening wave vector @xmath231 . here
the impurity density @xmath228 is used .
, title="fig:",width=211 ]    in figs .
[ fig_13 ] and [ fig_14 ] we show our calculated resistivity for gamnas samples .
we use the parameters corresponding to gamnas : dielectric constant @xmath232 , hole effective mass @xmath233 , impurity density @xmath234 , and the magnetic coupling constant @xmath235120 mev nm@xmath236 .
[ fig_13 ] shows normalized resistivity @xmath226 as a function of temperature @xmath237 for several values of hole density . in fig .
[ fig_13 ] there are two independent sources of temperature dependence in our calculated resistivity
 one source is the energy averaging defined in eq .
( [ tau ] ) and the other is the explicit temperature dependence of the finite temperature rpa screening function @xmath238 . since the fermi temperature is much higher than the critical temperature ( @xmath239 ) the screening function is weakly temperature dependent above the critical temperature ( unpolarized state ) .
thus the decrease of the resistivity above @xmath2 with increasing temperature arises from the thermal energy averaging .
this is a well - known high - t effect of ionized impurity scattering in semiconductors : @xmath240 decreases with increasing @xmath13 due to the thermal averaging . as density decreases this effect
is enhanced . in the ferromagnetic polarized state ( @xmath1 ) , as the temperature decreases from @xmath16 , the screening length increases until all the holes are polarized . in this temperature range
we find a strong temperature dependent resistivity .
this strong temperature dependence arises from the low temperature screening function and the change of the carrier densities of the each subband . when all the carriers are polarized ( i.e. at low density and temperature ) the screening function is almost independent of temperature and the resistivity is also therefore temperature independent .
the interplay between the screening length and the down spin (  ) carrier density gives rise to the minimum of the resistivity in the low density regime .
the decrease of the scattering time due to the fast increase of the down spin density overwhelms the increase of the scattering time due to the decrease of the screening length , producing the non - monotonic behavior for @xmath241 and @xmath242 in fig . 13 .
but in the high density regime , where the spins are not fully polarized even at low @xmath13 , screening is the dominant effect on the temperature dependent resistivity .
as the density increases the relative low temperature resistivity , @xmath226 , decreases until the holes are partially polarized .
the change of screening wave vector is larger in this case , leading to the monotonically increasing @xmath240 in the @xmath243 regime .    in fig .
[ fig_14 ] we show our calculated resistivity as a function of temperature for two hole densities ( a ) @xmath227 0.1 and ( b ) @xmath2300.2 with finite temperature rpa screening function . the temperature dependence of the impurity ( i.e. mn ) moment magnetization as well as the hole spin polarization is given in the insets , which is calculated using the weiss molecular mean - field theory for delocalized carriers @xcite .
note that the hole gas is almost fully spin polarized upto @xmath244 at low density , but at higher density the holes are partially polarized even at @xmath82 .
the bottom insets show the calculated temperature dependence of the finite temperature screening wave vector @xmath231 . at high density
the screening wave vector changes by 6% when the temperature increases from zero to @xmath2 due to the partial polarization of the holes at @xmath82 .
but at low density the decrease of the screening wave vector is about 20% . in the metallic gamnas samples
the change of the resistivity when the temperature goes from @xmath2 to zero is about 20% in good agreement with our calculation .
similarly the observed decrease of @xmath240 for @xmath0 also arises naturally in our theory as a consequence of thermal averaging .
thus , the temperature dependence of the resistivity in the metallic dms samples may be arising almost entirely from the temperature dependence of screening and thermal averaging in the charged impurity scattering .
we have developed in this work two complementary theories for understanding dms transport properties .
our work establishes that dms transport , even in ideal intrinsic circumstances , is rather complex , and depends sensitively on many system parameters including the exchange coupling , the magnetic impurity density , the carrier density , the temperature , the band parameters of the parent semiconductor ( e.g. effective mass , etc . ) , and the details of the charged impurity distribution ( and therefore compensation ) .
given such a complex parameter space , it is not meaningful to try to develop quantitative transport theories at this early stage of the subject since all the intrinsic parameters may not be known .
we therefore focus in this work on developing a comprehensive qualitative theoretical description which emphasizes general broad features of how various parameters affect dms transport behavior . as such
we have concentrated in this work on understanding temperature and carrier density dependence of dc resistivity in a model dms system , keeping primarily the extensively studied ga@xmath4mn@xmath5as system in mind . even for gamnas ,
the transport data for various values of mn concentration ( @xmath245 ) cover much too broad a range of behavior for a unified and coherent theoretical description .
for example , low ( and sometimes large ) mn concentrations ( @xmath246 and sometimes @xmath247 ) are known to lead to insulating transport behavior usually attributed to a disorder - driven metal - insulator ( anderson ) localization .
we neglect _ all _ localization effects in our theory .
the localized gamnas regime in all likelihood requires its own characteristic theoretical description such as the bound magnetic polaron percolation theory @xcite whose transport properties @xcite have recently been theoretically analyzed .    even without the disorder - induced strong localization complications , neglected completely in this work , we face the formidable difficulty of using the semiconductor valence / conduction band ( the valence band for gamnas , where the carriers mediating the ferromagnetic interaction are holes ) or the magnetic impurity induced impurity band ( i.e. mn induced d - band in the fundamental band gap of gaas ) picture for describing the carrier dynamics . the precise nature (
i.e. valence band versus impurity band ) of the dms carriers is still a controversial issue although it is likely that at the very high doping densities ( e.g. mn density @xmath248 in gamnas ) of dms interest the impurity band overlaps strongly with the valence band ( i.e. forms the tail of the valence band ) , and therefore , the distinction between the valence and the impurity band picture is not a real qualitative difference .
our dmft theory , presented in sections ii - iv of this paper , clearly shows that in the strong exchange coupling ( @xmath249 ) regime the impurity band picture applies whereas in the weak - coupling regime ( @xmath250 ) , the valence band picture applies .
it is possible that gamnas belongs to the intermediate coupling regime ( @xmath251 ) , where it may be more appropriate to think of the holes to be residing in the extended tail of the valence band , presumably with an enhanced effective mass compared with the gaas valence band hole mass .
such a coupled impurity - valence band picture of gamnas is consistent with recent optical spectroscopy measurements @xcite , but more experimental work is needed to settle this question .
the theoretical strength of our dmft description is that , being a nonperturbative technique , it can handle both strong - coupling and weak - coupling regimes , and our results presented in figs .
6  10 of this paper show qualitative difference between the strong - coupling regime ( @xmath252 ) with an impurity band well - separated from the semiconductor band and the intermediate - coupling regime ( @xmath253 ) with only band - tailing and no separate impurity band formation .
temperature , carrier density , and impurity concentration all play qualitatively important roles in determining the dc transport properties within dmft , and sorting out the details with respect to experimental results may be extremely difficult .
the weakness of dmft is that it can only include spin - disorder scattering ( controlled by @xmath10 ) and mn impurity induced local potential scattering ( controlled by @xmath25 ) effects on transport properties . as such it leaves out the most important scattering mechanism which may be operational in real samples , namely scattering by charged impurities which is often the most important resistive scattering process in heavily doped semiconductos below the room temperature ( or the optical phonon scattering regime which may well be above the room temperature ) .
the reason dmft is unable to account for charged impurity scattering is that the coulombic impurity potential is long ranged , and dmft by construction is a local theory .
thus , rewriting our starting hamiltonian ( eq .
1 ) more completely we have @xmath254,\ ] ] where @xmath255 is the carrier - charged impurity interaction and @xmath256 is the carrier - carrier ( i.e. hole - hole in gamnas ) interaction . in principle , the terms ( i.e. @xmath255 , @xmath256 ) within the square bracket are parts of the @xmath18 , but it is important to appreciate their considerable ( perhaps even dominant ) importance in determining the dc transport properties . to include the charged impurity scattering effects on transport , we use the highly successful and robust semiclassical boltzmann transport theory to dms systems assuming a mean - field approach where the long - range coulomb impurity potential arising from @xmath255 ( assuming random impurity scattering ) is screened by the polarization bubble diagrams arising from @xmath256 .
this type of boltzmann transport theory is extremely successful in describing semiconductor transport properties @xcite .
we note that our dmft and boltzmann transport theories are complementary
dmft treats the spin disorder and the local potential scattering associated with mn dopants and boltzmann theory treats the scattering by screened charged impurity scattering .
the magnetic dms properties enter the boltzmann theory only indirectly through the carrier spin polarization calculations .
spin disorder scattering is _ not _ explicitly included in the boltzmann theory although it is straightforward to do so .
our boltzmann theory manifests nontrivial interplay among temperature dependent screening , temperature dependent spin polarization ( i.e. spin up - down carrier densities ) , and thermal energy averaging , leading to temperature dependent resistivity ( figs .
13 and 14 ) that are rather similar to experimental observations @xcite in gamnas . based on this qualitative similarity
we conclude that much of gamnas dms transport is dominated by screened charged impurity scattering with spin disorder scattering playing only a rather minor quantitative role .
recently , lopez - sancho and brey @xcite have come to a very similar conclusion for gamnas transport properties .
schematic diagram for the experimental gamnas dc - resistivity for decreasing hole density ( from bottom to top ) ( see ref . [ 3 ] ) .
, width=211 ]    before concluding we now discuss our theoretical results in light of the existing dms transport data . although there are some experimental transport data in other dms systems ( most notably inmnas with qualitatively similar behavior to gamnas ) , truly extensive and reliable transport data@xcite are available only for ga@xmath4mn@xmath5as ( in the @xmath257 regime ) dms samples . even for this well - studied system ,
experimental transport results are problematic , and have considerable spread in the sense that nominally `` identical '' gamnas samples ( i.e. same nominal carrier density and mn concentration made in the same growth run ) may have different @xmath2 and transport properties .
this situation is improving as sample quality and processing ( e.g. annealing ) techniques improve , but experimental dms transport properties are still not robust in a quantitative way .
this is of course very understandable given the very large parameter space ( i.e. exchange coupling , hole density , mn concentration , defects and impurities , compensation , band structure parameters , etc ) that dms transport depends on . with these serious caveats we show in fig .
15 a schematic depiction of the generic experimental observation for @xmath240 in gamnas for various hole densities . at high hole density
the system shows `` metallic '' behavior for @xmath1 with @xmath240 increasing somewhat with @xmath13 upto @xmath2 , and then decreasing slowly for @xmath0 .
even the optimally doped most `` metallic '' gamnas is , however , at best a bad metal with mobilities of the order of 10 @xmath258 with @xmath259 where @xmath260 is the transport mean free path .
it is important to realize that dms transport is always extremely highly resistive due to the very large amount of impurities and defects invariably present in the low temperature mbe process needed for producing homogeneous gamnas samples .
thus , from the perspective of heavily doped semiconductors , although the dms systems may be above the mott limit ( i.e. the carrier density in the mott metallic range ) , they are close to being anderson insulators due to strong disorder effects .
as hole density decreases the system eventually becomes an insulator at low enough carrier densities ( the top curve in fig .
15 ) , with @xmath240 decreasing monotonically increasing @xmath13 .
two generic features of fig .
15 are : ( a ) a peak ( or a kink ) in @xmath240 near @xmath261 , and ( 2 ) the slow decrease of @xmath240 for @xmath0 . both of these generic features of @xmath240 as well as the `` metallic '' high - density behavior are qualitatively ( perhaps even semi - quantitatively ) well - explained by our boltzmann theory approach including only scattering by ionized impurity scattering .
this is apparent by comparing fig .
15 with figs .
14 and 13 where our boltzmann theory results are shown .
physically , the increasing @xmath240 with @xmath13 ( @xmath262 ) arises from the decreasing strength of screening due to the interplay of two independent and competing physical effects : temperature induced suppression of screening and spin polarization ( i.e. the spin polarization decreasing with increasing temperature ) induced enhancement of screening with increasing temperature . as explained in sec .
v the competition between these two effects depend on the carrier density , leading to some weak non - monotonicity in @xmath240 for @xmath1 . in this screening picture ,
the resistivity peak or cusp at @xmath263 arises from the ferromagnetic to paramagnetic transition which affects screening as the carrier density of states ( which is inversely proportional to the screening length ) changes from one in the fully spin polarized state to two in the fully paramagnetic phase .
thus , even without any spin disorder scattering effects , just screened ionized impurity scattering by itself will give rise to the peak or the kink in @xmath240 at @xmath16 .
we believe that spin disorder scattering , which also produces a kink at @xmath263 ( see , for example , figs . 7 and 9 ) , play only a minor role in the resistivity `` maximum '' at @xmath2 in optimally metallic gamnas with most of the peak structure arising from the screening properties of ionized impurity scattering .
the second generic experimental feature in fig .
15 , the slow decrease of @xmath240 for @xmath0 , can not be explained _ at all _ by spin disorder scattering since spin disorder should remain large in the paramagnetic system ( @xmath99 ) and certainly should not decrease with increasing @xmath13 .
our boltzmann theory provides a natural explanation ( figs .
13 and 14 ) for the decreasing @xmath264 as arising from the energy averaging of the relaxation time ( eq .
19 ) , i.e. , @xmath240 decreases with increasing @xmath13 simply because the holes move `` faster '' at higher temperatures ( i.e. increasing kinetic energy with increasing @xmath13 ) .
this decreasing @xmath240 with increasing @xmath13 ( @xmath265 ) also shows that our neglect of phonon scattering in the transport theory is a valid approximation since phonon effects will always increase @xmath240 with increasing @xmath13 .
@xmath240 will increase again when phonon scattering starts dominating over ionized impurity scattering at much higher temperatures .
the relative lack of importance of phonon scattering in dms systems arises from their very strong charged impurity resistive scattering effects as reflected in very small sample mobilities .
in this paper we investigate the transport properties of the diluted magnetic semiconductors using dynamical mean field theory and boltzmann transport theory .
we have shown that the resistivity depends strongly on the system parameters , i.e. , exchange coupling , carrier density , doping , and temperature . the resistivity drop with decreasing temperature in the ferromagnetic state
can be partially explained by the screening theory for metallic samples .
the parameter dependence of the resistivity contains important information about the physics of diluted magnetic semiconductors .
we find that in the strong exchange coupling regime the spin disorder scattering and the formation of the bound state in the impurity band compete to produce an unusual behavior in the temperature dependent resistivity .
we also show that in the weak coupling regime the occupation of the minority spin band is critical to the scattering mechanisms , and substantially reduces the resistivity because the repulsive interaction between local moments and `` wrong - spin '' carriers suppresses the carrier amplitude at the impurity site , reducing the effective carrier - spin coupling .
our boltzmann transport theory for charged impurity scattering is good qualitative agreement with the existing dms experimental data , showing that transport in dms gamnas system may very well be dominated primarily by screened ionized impurity scattering ( with spin disorder scattering playing only a minor secondary role ) , at least for the optimally doped metallic gamnas samples .
we have completely neglected detailed band structure complications ( e.g. spin - orbit coupling in the valence band ) in our theory .
these effects are certainly very important , but our interest in this paper is the development of a conceptually coherent qualitative theory for dms transport identifying the main transport mechanisms , and as such we have neglected all nonessential complications .
we thank darpa and us - onr for support .
f. matsukura _
et al . _ , ,
r2037(1998 ) ; s.j .
potashnik , k.c .
chun , j.j .
berry , n. samarth , and p. schiffer , appl .
. lett . * 79 * , 1495 ( 2001 ) ; sh.u .
yuldashev , h. i m , v.sh .
yalishev , c.s .
park , t.w .
kang , s. lee , y. sasaki , x. liu , and j.k .
furdina , appl .
phys . lett . * 82 * , 1206 ( 2003 ) .
et al . _ ,
lett . * 79 * , 3473 ( 2001 ) ; n. theodoropoulou , a. f. hebard , m. e. overberg , c. r. abernathy , s. j. pearton , s. n. g. chu , and r. g. wilson phys . rev . lett . * 89 * , 107203 ( 2002 )