diluted magnetic semiconductors ( dms ) based on iii - v semiconductors doped with mn have attracted a lot of interest recently , after critical temperatures for the onset of ferromagnetism of the order of 110 k have been found in ga@xmath0mn@xmath1as , for @xmath2.@xcite more recently , critical temperatures larger than room temperatures have been reported in mn - doped gan , enhancing the hope for extensive technological applications of these materials.@xcite    while there is general agreement that the magnetism is due to charge - carrier mediated , effectively ferromagnetic interactions between the mn spins , there are various theoretical models attempting to understand their detailed behavior .
the dms are alloy systems , with inherent positional disorder of the mn atoms ; further , spin - orbit coupling may play a significant role for hole doping . a theoretical treatment which takes into account all these factors , and their effects on the magnetic and transport properties of dms , is not yet available . instead , theoretical models have tended to concentrate on different aspects of the problem .
one class of models , where much work has been done , neglects both the disorder and spin - orbit coupling effects . at large carrier densities , where coulomb potentials of the impurity ( mn ) atoms are effectively screened out , and other disorder effects ( such as those coming from as antisite defects which are believed to be the cause of the large compensation observed experimentally ) can be neglected ,
the former at least can be justified .
in such a case , holes occupy a fermi sea at the top of the valence band.@xcite studies of the spin - wave spectra@xcite , monte - carlo studies@xcite as well as dynamical mean - field theory studies@xcite have been performed .
overall , their results are rather similar to the physics one would expect to hold in a conventional ferromagnet .
more recently , it has been proposed that spin - orbit coupling of the valence band , in the high carrier density ( strongly metallic ) regime leads to an rkky coupling between mn moments that is _ anisotropic _ in spin space , and where the positional disorder effectively leads to random anisotropy.@xcite this would lead to frustration effects on the magnetic ordering , and therefore unusual ferromagnetism .
a number of very useful ab - initio studies have also been published.@xcite    in parallel , we have studied effects of positional disorder on the magnetic ordering ( in the absence of spin - orbit coupling ) of an impurity band model,@xcite developed from previous work done in the context of insulating ii - vi dms compounds.@xcite such an approach should be of relevance at low doping concentrations x , at and below the metal - insulator transition , and possibly even above .
evidence for the existence of impurity - like states is provided by ab - initio studies,@xcite which found that occupied hole states near the fermi level have wave - functions mostly concentrated on and near the mn impurity .
more recently , angle resolved photoemission spectroscopy has revealed the existence of the impurity band in ga@xmath0mn@xmath1as with @xmath3 ( very close to the metal - insulator transition).@xcite a scanning tunneling microscope study @xcite demonstrates the existence of an impurity band in ( ga , mn)as samples with x=0.005 - 0.06 .
optical spectroscopy@xcite also identifies the impurity band for x=0.0001 and x=0.05 samples .
electron densities in localized states , as well as states close to the mobility edge , can be far from homogeneous , unlike bloch waves . at low mn concentrations and high degree of compensation ( which comes from charged centers ) seen in dms , short length - scale density fluctuations could be significant .
this can lead to inhomogeneities in the local charge densities at different mn sites , which in turn leads to ( microscopically ) inhomogeneous magnetizations of the magnetic ions in the ordered phase .
such inhomogeneity would be expected to alter the nature of the collective magnetic excitations ( spin waves ) of the magnetically ordered system , which in turn would affect charge transport through the magnetic excitation processes which give rise to spin - flip scattering .    in this paper
, we develop a numerical scheme to calculate the spin wave spectrum for a finite but large size model of coupled fermions and spins , in the presence of quenched disorder .
we show that the method accurately reproduces the results for a lattice model of small size obtained from a standard treatment of spin waves via the random phase approximation.@xcite we then apply the method to the model of ref . , and study in detail its collective magnetic excitations .
the plan of the paper is as follows .
the model hamiltonian for dms is described in section ii .
section iii is devoted to developing the numerical scheme .
the accuracy of the scheme is tested by comparing it with the standard rpa method in section iv .
in particular , we demonstrate that we can use our scheme to study large size systems , beyond the capability of the standard rpa approach .
section v presents the results obtained by applying the scheme to the model of section ii , including the density of states and nature of the eigenvectors of the magnetic excitations of the coupled fermion - spin system .
we summarize our conclusions in section vi .
in the following we develop the formalism for an impurity band hamiltonian that we believe to be appropriate for describing ( at least qualitatively ) the diluted magnetic semiconductor ga@xmath0mn@xmath4as for low doping @xmath5 ( in the insulating phase , or not too far from the metal - insulator transition ) .
however , generalizations of this method for other types of hamiltonians , crystal structures , parameters etc . can be done in a straightforward manner .    the iii - v host semiconductor is assumed to have the zinc - blende structure appropriate for gaas .
experimental measurements@xcite suggest that valence - ii mn substitutes for valence - iii ga . as a result ,
the @xmath6 mn dopants are randomly distributed at positions @xmath7 , @xmath8 on the @xmath9 fcc ga sublattice , corresponding to @xmath10 .
all throughout this paper , we assume periodic boundary conditions .    due to the valence mismatch , each mn introduces a charge carrier ( hole ) in the system .
however , experimentally it is found that the hole concentration is considerably suppressed through compensation processes . as a result ,
the number of holes is only a fraction @xmath11 of the number of mn atoms @xmath12 .
each mn atom also has a @xmath13spin @xmath14 from its half - filled @xmath15 level .
the magnetic properties of the doped semiconductor are related to the exchange interaction between the mn spins and the charge carrier spins .
this interaction is known to be antiferromagnetic ( afm ) , and proportional to the probability of finding the charge carrier at the corresponding mn site .
recently , we proposed a simple impurity - band model to describe the behavior of the charge carriers for the low - doping regime.@xcite the main justification is that near and below the metal - insulator transition ( @xmath16 ) the density of charge carriers is not large enough to effectively screen the attractive coulomb potential of the mn dopants . as a result , bound
, hydrogen - like impurity states are created about each mn site , at an energy @xmath17 ry above the top of the valence band . due to interactions ,
these impurity states lead to the appearance of an impurity band , and the holes first occupy states in this band . only if the concentration of holes ( or the temperature ) is large enough , are states in the valence band itself occupied by holes .
thus , it is reasonable to attempt a description of the charge carriers behavior only in terms of this impurity band .    as in previous work,@xcite in the following we will use the electron formalism to treat the problem .
this is equivalent to performing a particle - hole transformation which leads to an inversion @xmath18 of the charge carrier spectrum .
thus , instead of emptying the top of a valence - like impurity band ( i.e. , introducing holes in the system ) we instead occupy the bottom of a conduction - like impurity band .
we also make the simplifying assumption that the isolated impurity wave - function for a charge carrier trapped near the mn at site @xmath7 is the @xmath19 hydrogen orbital @xmath20 , where @xmath21   is the appropriate bohr radius related to the effective mass of the heavy hole ( @xmath22 for gaas ) and the binding energy ( @xmath23 mev for mn in gaas).@xcite our approach neglects both the complicated orbital form of the acceptor wave - function ( ref . ) and spin - orbit coupling .
the former is not expected to lead to any qualitative changes .
it has recently been proposed that spin - orbit coupling leads to frustration in the magnetic ordering.@xcite these effects are left out in the present study , which concentrates on the effect of disorder .
a study combining both effects is indicated for future work .
we consider the hamiltonian @xmath24 here , @xmath25 creates a charge carrier with spin @xmath26 in the impurity state centered at site @xmath7 .
the first term describes the hopping of charge carriers between impurity states .
we use the parameterization @xmath27 ry , where @xmath28 , of magnitude and form appropriate for hopping between two isolated 1@xmath29 impurities which are not too close to one another .
@xcite this hopping term has been shown to lead to the appearance of an impurity band which has a mobility edge , as well as a characteristic energy for the occupied states in agreement with physical expectation.@xcite    the second term of the hamiltonian ( [ 1.1 ] ) describes the afm exchange between the mn spin @xmath30 and the charge carrier spin @xmath31 [ @xmath32 are the pauli spin matrices ] .
this afm exchange is proportional to the probability of finding the charge carrier trapped at @xmath33 near the mn spin at @xmath34 , and therefore @xmath35 .
based on calculations of the isolated mn impurity in gaas , we estimate the exchange coupling between a hole and the trapping mn ( @xmath36 ) to be @xmath37 mev.@xcite as already stated , the number of mn atoms is @xmath6 , and therefore there are a total of @xmath38 states in the impurity band .
the number of charge carriers is fixed to @xmath39 , where we take @xmath40 .    in ref .
we studied the relevance of various other terms , such as an on - site potential ( associated with the coulomb potential of the compensation centers ) , hubbard - like on - site repulsion ( describing interactions between charge - carriers ) , external magnetic field etc .
while they lead to various _ quantitative _ changes , we believe that the _ qualitative _ picture we present in the following sections is not changed by their absence .
the random phase approximation ( rpa ) describes the collective excitations of a system about its self - consistent @xmath41 mean - field ground state . in order to clarify the notation used , we begin this section with a very short review of the derivation of the relevant equations for the mean - field ground state .
we use the equation - of - motion approach of ref . .
at the mean - field level , we are trying to find non - interacting fermionic quasiparticles through a unitary transformation of the on - site charge carrier operators : @xmath42 such that @xmath43 \approx e_{n\sigma } a^{\dagger}_{n\sigma}\ ] ] and @xmath44 \approx - \hbar h_i s_i^{+}\ ] ] for a non - interacting system , eqs .
( [ 22 ] ) and ( [ 23 ] ) would be exact . for an interacting system at the mean - field level , the exact hamiltonian is approximated by the diagonalized , non - interacting hartree - fock hamiltonian @xmath45 for which eqs .
( [ 22 ] ) and ( [ 23 ] ) become exact .
it is straightforward to show that for hamiltonian ( [ 1.1 ] ) @xmath46 = \sum_{ij}^ { } t_{ij } \psi_{n\sigma}(j ) c^{\dagger}_{i\sigma } + \sum_{ij}^ { } j_{ij } \vec{s}_j { \vec{\sigma}_{\alpha\sigma } \over 2 } \psi_{n\sigma}(i ) c^{\dagger}_{i\alpha}\ ] ] the right - hand - side of eq .
( [ 25 ] ) contains spin operators , so it can not be put in the form of eq .
( [ 22 ] ) unless we linearize it , by replacing the spin operators with their mean - field ground state expectation value @xmath47 .
while the choice of a collinear ground - state , with all mn spins aligned along the @xmath48-axis , is not the most general possibility , we have actually shown in ref .
that the ground - state of this model , for the range of parameters we are interested in , is indeed collinear .
this is why we choose to apriori make this assumption in this case .
after the linearization , we can now equate the right - hand - side of eq .
( [ 25 ] ) with @xmath49 [ see eq .
( [ 22 ] ) ] to obtain the hartree - fock equation for the one - particle orbitals : @xmath50 the mean - field ground - state of the charge carriers is obtained by occupying the @xmath51 lowest - energy orbitals , @xmath52    the effective magnetic fields @xmath53 [ see eq .
[ 24 ] ) ] are obtained in a similar way .
we evaluate the commutator of the exact hamiltonian with the spin raising operator @xmath54 @xmath55 = -\hbar \sum_{j}^{}j_{ij } s_i^z s_i^{+ } + \hbar \sum_{j}^ { } j_{ij } s_i^{+ } s_i^{z}\ ] ] the right - hand - side now contains charge carrier operators ( in the charge carrier spin operator ) , beside the spin operators .
we again use linearization @xmath56 [ see eq .
( [ 27 ] ) ] , @xmath57 , and by comparison with eq .
( [ 23 ] ) we find @xmath58    from eq .
( [ 24 ] ) we see that if @xmath59 , then in the mean - field ground state the spin @xmath14 is in the `` up '' state @xmath60 and therefore @xmath61 , while if @xmath62 the spin @xmath14 is in the `` down '' state @xmath63 and @xmath64 ( @xmath65 ) . as a result
, the spin - part of the mean - field ground - state can also be easily found if the expectation value of the charge carrier spins @xmath66 are known from eqs .
( [ 26 ] ) , ( [ 27 ] ) .
thus , we obtain the usual self - consistent hartree - fock equations .
diluted magnetic semiconductors exhibit ferromagnetism at low temperatures . as a result
, we assume that in the mean - field ground state , all mn spins are fully polarized @xmath61 . then , from eq .
( [ 26 ] ) we find the lowest - lying @xmath51 states .
if the charge carriers are also fully polarized at @xmath41 , all the occupied @xmath51 state have @xmath67 , and therefore the mean - field ground - state is of the form : @xmath68    of course , one must check that self - consistency is obeyed by verifying that indeed the first @xmath51 lowest energy one - electron states are all spin - down .
we find that this always holds true for the parameters we study in this paper .
as discussed later , the spin - wave spectrum of the collective excitations also confirms that this ground - state is indeed stable .
however , for higher charge carrier concentrations ( or hamiltonians with other types of interactions ) it is likely that the hf ground - state is only partially spin polarized . in that case , one must do the full iterational search for the self - consistent hf ground state .
the computation for the spin - wave spectrum in the partially - polarized case ( for instance due to spin - orbit coupling ) , can be derived in a similar way to the one we present in the following for the fully polarized case .
we will use the same equation - of - motion approach to derive the rpa equations .
we are interested in spin - waves , which are related to spin - flip ( spin - lowering ) processes .
therefore , the rpa operators for spin - wave collective modes must be constructed in terms of spin - flip operators @xmath69{c } b^{\dagger}_{hp } = a^{\dagger}_{h\uparrow } a_{p\downarrow } \\
b^{\dagger}_i = { 1 \over \sqrt{2s } } s^-_i \\ \end{array } \right.\ ] ] the index @xmath70 runs over the empty ( `` hole '' ) states of the hf ground state , while the index @xmath71 runs over the occupied ( `` particle '' ) states of the hf ground state .
the index @xmath72 runs over all the mn spin positions .
the hf ground - state @xmath73 is fully polarized , and given by eq .
( [ 30 ] ) .
the spin - flip operators have `` bosonic '' nature , as expected for collective modes operators .
indeed , it is straightforward to verify that @xmath74\rangle = \delta_{pp ' } \delta_{hh'}$ ] , @xmath75\rangle = \delta_{ij}$ ] , with all the other commutators identically zero .
here , the notation @xmath76 signifies an average over the hf ground state @xmath77 .
the hamiltonian can be rewritten in terms of the @xmath78 operators as [ see eqs .
( [ 1.1 ] ) , ( [ 21 ] ) ] @xmath79 where @xmath80 and @xmath81    one can easily derive the equations of motion for the spin - flip operators . for instance @xmath82= { 2 \over \sqrt{2s } } \sum_{nm}^ { } j_{n\downarrow , m\uparrow}(i ) a^{\dagger}_{n\downarrow } a_{m\uparrow } s^z_i\ ] ] @xmath83 since rpa is an approximation , not an exact solution , we again must linearize this commutator about the hf ground - state , and keep only the meaningful terms .
in other words , we replace @xmath84 here , the first transformation is the linearization about the hf ground state @xmath73 .
the second transformation results if one uses the hf expectation values @xmath85 , @xmath86 .
the restriction to relevant terms means that in any sum over @xmath87 , we restrict @xmath88 to being an occupied orbital in the hf ground - state @xmath89 , because otherwise @xmath90 , i.e. a spin - flip process from the hf ground - state is impossible [ see eq .
( [ a.1 ] ) ] . performing a similar linearization for @xmath91 and collecting the various terms ,
we find the linearized equation of motion @xmath92\approx h_i b_i^{\dagger } + \sqrt{2s } \sum_{p=1}^{n_h } \sum_{h=1}^{n_d } j_{p\downarrow , h\uparrow}(i ) b_{hp}\ ] ] after similar steps , the linearized equation of motion for the other spin - flip operator @xmath93 is found to be @xmath94 \approx -\left ( e_{h\uparrow } - e_{p\downarrow}\right ) b_{hp } - \sqrt{2s } \sum_{i}^ { } j^*_{p\downarrow , h\uparrow}(i ) b_i^{\dagger}\ ] ]    eqs .
( [ a4 ] ) and ( [ a5 ] ) allow us to `` guess '' the rpa part of the hamiltonian , for which eqs .
( [ a4 ] ) and ( [ a5 ] ) become exact , to be given by @xmath95 @xmath96 thus , the rpa hamiltonian is quartic in electron operators and quadratic in spin operators , showing that it is the next order term after the hartree - fock hamiltonian ( which is quadratic in electron operators and linear in spin operators ) @xmath97 higher order approximation terms describe interactions between the collective modes , which are neglected at the rpa level .
we want to diagonalize the rpa hamiltonian to the canonical form @xmath98 where @xmath99 is the creation operator for a spin - wave with energy @xmath100 , and @xmath101 is an index ( in homogeneous systems , it would be a wave - vector ) .
since they create collective mode excitations , these operators must have bosonic nature @xmath102=0 $ ] . however , since they do not describe exact , but only approximative solutions of the exact hamiltonian , in fact only the weaker condition @xmath103\rangle=\delta_{\alpha\beta}$ ] is obeyed .
the most general form for these operators [ see eq .
( [ a.6 ] ) ] is @xmath104 using the equation of motion @xmath105= \hbar \omega_{\alpha } q^{\dagger}_{\alpha } $ ] and computing the commutator using eqs .
( [ a.6 ] ) , ( [ a.8 ] ) , we find the diagonalization condition to be @xmath106 x^{(\alpha)}_{hp } = \sqrt{2s } \sum_{i}^{}j_{p\downarrow , h\uparrow}(i ) y^{(\alpha)}_i \\ \label{a.9b } & & \left [ \hbar \omega_{\alpha } - h_i\right ] y^{(\alpha)}_i = -   \sqrt{2s } \sum_{p=1}^{n_h}\sum_{h=1}^{n_d }   j_{p\downarrow , h\uparrow}^*(i ) x^{(\alpha)}_{hp}\end{aligned}\ ] ]    these are the rpa equations and they can be recast in the customary standard eigenvalue rpa equation form@xcite @xmath107{c } { \bf x } \\ { - \bf y } \\ \end{array } \right ) = \left ( \begin{array}[c]{cc } { \bf e } & { \bf j}\\ { \bf j}^ * & { \bf h } \\ \end{array } \right ) \cdot \left ( \begin{array}[c]{c } { \bf x } \\ { \bf y } \\
\end{array } \right)\ ] ] where the vectors @xmath108 and @xmath109 contain all the unknowns @xmath110 and @xmath111 , and the matrices @xmath112 and @xmath113 have the elements @xmath114 , @xmath115 and @xmath116 .
the dimension of the rpa matrix ( [ 5.2 ] ) , and therefore the number of normal modes , is @xmath117 .
of these solutions , @xmath6 are proper spin - wave collective modes , while the rest of @xmath118 are spin - flip processes associated with particle - hole excitations .
if there were no interactions ( j=0 ) , the eigenenergies of these spin - flip processes would be @xmath53 for the lowering of the mn spin @xmath119 , respectively @xmath120 for a hole spin - flip [ see the left - hand side of eqs .
( [ a.9]),([a.9b ] ) ] .
interactions renormalize these values [ as described by the right - hand - side of eqs .
( [ a.9]),([a.9b ] ) ] , but one still expects @xmath6 spin - wave collective solutions at low energies , coming from the mn spin lowering , and @xmath121 spin - flip particle - hole excitations at energies comparable or larger than the spin - flip gap @xmath122 .
we now comment on the sign of the rpa spectrum frequencies @xmath100 . since @xmath59
( [ 29 ] ) ] and @xmath123 ( by definition of the hf ground - state ) , it is apparent from eqs .
( [ a.9]),([a.9b ] ) that @xmath6 spin - wave solutions will have positive energies @xmath100 , while the @xmath121 spin - flip particle - hole excitations will have negative energies .
this _ does not _ imply that the system is unstable , but simply that we have not chosen the proper spin - flip creation operators . instead of the choice of eq .
( [ a.8 ] ) , we can also denote @xmath124 in terms of these new operators , we have now @xmath125= \hbar \omega_{\alpha } p_{\alpha}$ ] , i.e. @xmath126 in other words , a negative solution for @xmath127 simply means that we chose the corresponding annihilation operator instead of the creation operator when we wrote eq .
( [ a.8 ] ) . for spin - flip processes , it is obvious that this is the case from the dependence of @xmath99 on @xmath93 , instead of @xmath128
.    an unstable mean - field ground - state is signaled by _
complex _ values of the spectrum frequencies @xmath129.@xcite excited states about the mean - field ground - state @xmath73 are of the general form @xmath130    @xmath131 , leading to a time dependence : @xmath132 |\psi\rangle,\ ] ] where @xmath133 is the hartree - fock energy of the system .
if any of the frequencies @xmath134 has a non - trivial imaginary part , it follows that the expectation value of any operator @xmath135 will move in time exponentially away from its mean - field ground - state expectation value @xmath136 , i.e. the mean - field is unstable to small perturbations .
the advantage of the standard rpa approach is that by solving the eigenvalue problem ( [ 5.2 ] ) , one finds the normal mode spectrum @xmath127 and the spatial distribution of the spin - waves ( related to the @xmath111 coefficients ) at once .
the obvious disadvantage is of a numerical nature : for a disordered system the rpa matrix must be diagonalized numerically .
the size of the matrix is @xmath137 .
the typical sizes we consider are systems with around @xmath138 mn spins , and @xmath139 holes ( although concentrations up to @xmath140 might have to be considered , depending on the doping @xmath5 ) , leading to rpa matrices of typical size @xmath141 . as a result
, one has to either consider much smaller systems ( in which case , finite size effects may be overwhelmingly important ) , or to try sparse matrix techniques to obtain the first few low - energy modes of the rpa matrix ( although the off - diagonal @xmath113 matrix is not necessarily sparse )
.    there is , however , an alternative way of finding the collective mode spectrum and spatial distributions . we can directly solve for the @xmath110 coefficients from eq .
( [ a.9 ] ) and rewrite eq .
( [ a.9b ] ) as @xmath142 where @xmath143 the advantage of this formulation is that one has to deal with much smaller matrices ( the dimension of the @xmath144 matrix is @xmath145 ) .
( [ 5.3 ] ) tells us that for frequencies @xmath146 belonging to the spin - wave spectrum , the matrix @xmath147 has at least one eigenvector corresponding to an eigenvalue @xmath148 ( if the mode is degenerate , there are several such eigenvectors ) .
the corresponding eigenvectors give us the desired spatial mode distribution @xmath149 .
thus , the problem is reduced to sweeping the range of frequencies of interest ( for low - energy collective excitations , this is usually a fairly small range of frequencies near @xmath150 ) , and monitoring the eigenvalues of the @xmath151 matrix . the dependence on @xmath152 of the matrix elements
@xmath153 , and therefore of the eigenvalues @xmath154 is monotonic for small @xmath155 [ see eq . ( [ 5.4 ] ) ] . as a result , each equation @xmath156 has a single solution , and we expect to find @xmath6 collective mode eigenenergies , one for each eigenvalue .
the monotonic dependence of @xmath154 on @xmath152 also simplifies the search for the collective spectrum , since it is enough to evaluate the eigenvalues @xmath154 in the range of interest on a grid of step @xmath157 , and use linear interpolation to find the solutions of the equations @xmath156 .
as we show in the following by comparison with standard rpa ( for small system sizes ) we can very easily achieve relative errors less than @xmath158 .
the rpa spectrum has a second type of solutions , the spin - flip particle - hole continuum .
these solutions are at frequencies of the order of the gap between the last occupied and first empty hf states @xmath159 , and one expects @xmath121 such solutions . in principle , one can use the same method to find these energies as well , although the fact that they extend over a larger range of frequencies and that the @xmath152 dependence is very strong makes the search more difficult .
typically , one is only interested in the lower and upper limits of the spin - flip continuum , which can be found easily with the method just described .
finally , we would like to mention that it is always possible to reformulate the standard rpa equation ( [ 5.2 ] ) in a `` self - consistent '' form of the type shown in eq .
( [ 5.3 ] ) , with a matrix @xmath151 of much smaller dimension than the rpa matrix .
for instance , for hamiltonians describing interacting electron systems , the rpa matrix has a dimension @xmath160 , where @xmath161 is the number of occupied , and @xmath162 is the number of empty orbitals in the mean - field ground - state.@xcite for a system with a finite concentration of electrons @xmath5 on a lattice of linear dimension @xmath163 in a @xmath164-dimensional space , we have @xmath165 and @xmath166 , and the rpa matrix scales as @xmath167 . for such systems , one can always find an equivalent reformulation of the rpa equation in terms of a @xmath151 matrix of size @xmath168.@xcite this allows numerical handling of considerably larger systems , without having to resort to sparse matrix techniques , which often have issues related to instability .
in this section we show , by direct comparison against known calculations , as well as by direct comparison against solving the rpa matrix equation , that the formulation of eqs .
( [ 5.3 ] ) , ( [ 5.4 ] ) gives correct and numerically very accurate results . in the second part of the section we describe in more detail the numerical implementation as well as efficiency of our method .
we first apply our approach by solving eqs .
( [ 5.3 ] ) , ( [ 5.4 ] ) analytically for a simplified case .
we assume that the mn spins are arranged on an ordered superlattice , instead of having random positions .
then , the charge - carrier part of the mean - field hamiltonian is easily diagonalized in the @xmath169-space : @xmath170 where @xmath171 here , @xmath172 is the kinetic energy of the non - interacting electrons , where @xmath173 indexes all the neighboring sites and @xmath174 for which @xmath175 .
also , @xmath176 , where again @xmath177 for which @xmath175 . in the hf ground state , the states @xmath178 given by the condition @xmath179 are occupied . since
the @xmath51 electrons are equally distributed among the @xmath6 mn sites , and are fully polarized , the z - component of the total spin created by electrons at each site , within hfa , is @xmath180 , and therefore we find @xmath181 [ see eq .
( [ 29 ] ) ] .
using the fact that the solutions of eq . ( [ 5.3 ] ) must also have plane - wave structure @xmath182 , it follows that the index @xmath183 becomes a wave - vector , as expected , and we can reduce eqs .
( [ 5.3 ] ) and ( [ 5.4 ] ) to a self - consistent equation for each collective mode @xmath184 : @xmath185 here , @xmath186 .
the occupation number @xmath187 obviously comes from the sum over occupied states we had in eq .
( [ 5.4 ] ) . for a finite - size ordered lattice ,
one expects considerable degeneracy for each mode @xmath188 , due to invariance to various symmetry transformations of the @xmath169-vector . in most cases ,
the number of charge carriers is such that the last orbital ( of degeneracy @xmath189 ) is only partially occupied by @xmath190 charge carriers . in this case
, one must obviously choose @xmath191 for all these orbitals , and @xmath192 for lower , fully occupied orbitals , and @xmath193 for higher , empty orbitals .
otherwise , the translational invariance is broken and the collective spectrum is not indexed by a wave - vector .    the sum in eq .
( [ 55 ] ) can be performed if we assume that the dispersion at the bottom of the band is of quadratic form @xmath194 ( this is a reasonable approximation since we are interested in low filling fraction @xmath195 ) .
then , the sum over occupied states can be transformed to an integral which is straightforward to evaluate , leading to the solution : @xmath196 here , @xmath197 is the fermi velocity , where the fermi vector is given by @xmath198 the function @xmath199 is given by @xmath200 @xmath201 \ln { { \hbar\omega + \epsilon_{q } + j_0s- v_f q } \over { \hbar\omega + \epsilon_{q } + j_0s+ v_f q}}\ ] ] where @xmath202 .    for comparison purposes , we will assume a local afm interaction @xmath203 , leading to @xmath204 . in this case , the hamiltonian becomes identical to the hamiltonian used by konig .
al _ in ref . , provided that the effective mass @xmath205 in the dispersion relation is assumed to be the band effective mass .
we can easily solve eq .
( [ 56 ] ) in the asymptotic limit @xmath206 , to find @xmath207\ ] ] this is indeed identical with the asymptotic long - wavelength spin - wave spectrum obtained in ref . , and has the typical quadratic dependence of spin - waves in conventional ferromagnetic systems .      in this section
we briefly illustrate the accuracy and speed of our formulation of the rpa problem [ eqs .
( [ 5.3 ] ) and ( [ 5.4 ] ) ] , in comparison with the standard rpa approach [ eq .
( [ 5.2 ] ) ] .
we use a rather small system , with @xmath208 mn spins randomly distributed on a fcc lattice with linear size @xmath209 , corresponding to a mn concentration @xmath210 .
the number of charge carriers is @xmath211 , and the other parameters are as defined in section [ sec2 ] .
thus , the standard rpa involves the diagonalization of a non - hermitian matrix of dimension @xmath212 .
the largest eight eigenvalues @xmath154 of the @xmath151 matrix ( circles ) , evaluated on a grid with a step @xmath213 mev .
the full line is @xmath214 .
the spin - wave spectrum is given by the condition @xmath215 . ]
accuracy @xmath216 for all @xmath217 spin - wave frequencies @xmath218 , between the exact rpa ( standard ) values and values obtained in our formulation using a grid with @xmath213 [ see fig . [ fig1 ] ] . ]    in fig .
[ fig1 ] we show the @xmath152-dependence of the largest eight eigenvalues @xmath154 of the @xmath144 matrix ( circles ) , evaluated on a grid with a step @xmath213 mev . [ for real @xmath152 and @xmath219 , the matrix @xmath151 is hermitian and all eigenvalues @xmath154 are real , see eq .
( [ 5.4 ] ) ] .
the full line is @xmath214 , and the spin - wave spectrum is given by the condition @xmath156 .
one can see that the eigenvalues @xmath154 have monotonically increasing dependence on @xmath152 for small @xmath220 mev , and therefore the equation @xmath221 can have at most one solution for each eigenvalue , for small @xmath152 .
if each eigenvalue yields a solution , we have found all @xmath6 spin - wave frequencies .
if one or more eigenvalues do not intersect the @xmath222 line , this means that the mean - field ground - state is unstable ( there are _ complex _ spin - wave frequencies ) .
for all the cases and parameters we investigate , we found that the fully - polarized ground - state is stable for our model .
largest relative error @xmath216 among @xmath223 spin - wave frequencies @xmath218 , between the exact rpa ( standard ) values and values obtained in our formulation using a grid with @xmath224 and 1 mev .
the line is drawn to guide the eye . even for a very large grid ,
the relative error is still reasonably small . ]    to find the spectrum frequencies @xmath100 , we use linear interpolation over the interval @xmath157 where each @xmath154 curve intersects the @xmath214 curve .
this method avoids the need of evaluating the eigenvalues at too many @xmath152 points ( the most time - consuming part of the computation is the evaluation of the matrix elements of @xmath144 )
. the high degree of accuracy obtained with this method is demonstrated in fig .
[ fig2 ] , where we compare the exact rpa spectrum @xmath218 obtained through direct diagonalization of the rpa matrix [ see eq .
( [ 5.2 ] ) ] with the the values @xmath225 obtained through our formulation .
in fact , we plot the relative error @xmath226 for each spectrum frequency @xmath218 . as one can see from fig .
[ fig2 ] , the largest error is of the order @xmath227 , for the grid @xmath213 .
this suggests that one could use a much larger grid @xmath157 and still have a good accuracy .
indeed , we have computed the spectrum through both methods for five different realization of disorder , using various grid values @xmath157 in our formulation , and selected the largest relative error for each case ( from a total of @xmath223 values ) .
the results are shown in fig .
the error is found to scale as the square of the grid size .
even with a very large grid @xmath228 mev , ( which implies the evaluation of the eigenvalues @xmath154 in only very few points ) , we still get a relative error smaller than @xmath229 .
thus , one can use the grid step @xmath157 to optimize and considerably speed up the computation .
however , as the number of spins ( and eigenvalues ) @xmath6 increases , one must take into consideration other complications , as discussed below .
cpu time for the two rpa methods , for systems with @xmath230 and 216 spins .
the standard rpa ( full line ) is significantly more time consuming than our formulation ( dashed line ) . ]    in fig .
[ fig4 ] we compare cpu times for finding the rpa spin - wave spectrum using the standard vs. our rpa formulations .
all simulations were done on the same processor .
we used a small grid @xmath213 mev for our approach , and verified that all relative errors were less than @xmath158 .
we use systems with @xmath231 and 216 spins , randomly distributed on fcc lattices of linear size @xmath232 and 18 ( @xmath233 ) .
as expected , the computational time depends in a power - law manner on the size @xmath6 , with exponents 6.4 and 3.2 respectively for the standard and current approach .
clearly , our formulation can be successfully used for much larger sizes than those afforded by the standard approach .
eigenvalues @xmath154 for a system with @xmath234 random spins , evaluated for a grid @xmath157=0.05mev .
crossing and anti - crossings of neighboring eigenvalues is apparent .
the dotted line is @xmath214 .
since we are interested in the solutions of the equations @xmath235 , it is clear that too large a grid @xmath157 may lead to considerable errors . ]
one important aspect to keep in mind is that as the system size @xmath6 increases , neighboring eigenvalues become more closely spaced , and accidental crossings and anti - crossings appear .
some typical examples are shown in fig .
[ fig5 ] , for a disordered system with @xmath234 spins . tracking the crossing of the eigenvalues is essential , since such events lead to a change in the indexing of the eigenvalues from one grid point to the next .
we found that requiring continuity in first and second derivatives allows a unique identification of each continuous eigenvalue from one grid point to the next ( 3 steps are shown for each eigenvalue in fig .
[ fig4 ] ) .
in fact , the only relevant crossings are the ones that take place within the step @xmath157 where the eigenvalues intersects the @xmath214 line ( shown as a dotted line in fig . [ fig4 ] ) .
such an example is shown in the left panel of fig .
[ fig4 ] . for the largest system size investigated , of @xmath234 randomly distributed spins , we find that only 4.5% of the eigenvalues have such relevant crossings within @xmath157 of @xmath156
. the percentage decreases with decreasing @xmath6 , to 1.1% and 0.3% for @xmath236 216 and 125 , respectively .
this percentage also depends on the grid step @xmath157 ; for increasing @xmath157 the identification of crossings and anti - crossings becomes more difficult , leading to possibly large errors .
one can optimize the choice of the grid step @xmath157 by starting with a larger value .
the well separated eigenvalues ( such as the ones depicted in fig .
[ fig1 ] ) will provide unique identification and very accurate values for their corresponding spin - wave spectra values .
where considerable mixing and therefore non - linear variation of the eigenvalues @xmath154 is apparent , a finer mesh is necessary in order to correctly characterize their variation near @xmath214 . in all our simulations
, we use the grid @xmath213 mev , which allows for comfortable tracking of each eigenvalue and also is sufficiently small to allow us to approximate the variation of @xmath154 as being linear within each @xmath157 step .
these comparisons clearly demonstrate the accuracy and speed of our formulation of rpa , as compared to the standard rpa .
the biggest advantage , though , is that it can easily be applied to systems with large sizes , for which standard rpa is numerically cumbersome .
in this section , we use the rpa method described above to study the density of states ( dos ) , as well as nature ( extended or localized ) of the spin - wave spectrum of ( iii , mn)v dms described by the model introduced in section [ sec2 ] .
spin - wave dispersion @xmath237 for an ordered , simple cubic superlattice arrangement of mn spins inside the host semiconductor .
lattices of linear size @xmath238 and 150 , with a total of respectively @xmath239 and @xmath240 mn spins are considered .
this corresponds to @xmath241 and @xmath40 .
dispersion is plotted along three linear cuts in the brillouin zone : @xmath242 , @xmath243 and @xmath244 . ]
we study the spin - wave density of states for three types of arrangements of the mn spins inside the host semiconductor .
all samples studied correspond to @xmath241 and @xmath40 .
other parameters are as specified in section [ sec2 ] .    first ,
we consider fully ordered systems , in which the mn spins are arranged on a simple cubic superlattice . for @xmath241 ,
the superlattice constant is @xmath245 . in this case
, we can study the spin - wave dispersion and density of states for very large sizes , since we can use directly eq .
( [ 56 ] ) to compute the spin - wave frequency @xmath246 for each wave - vector @xmath247 inside the brillouin zone .
spin - wave dispersion obtained along three cuts in the brillouin zone is shown in fig .
the dispersion has quadratic behavior near the center of the brillouin zone , as expected from the discussion for an ordered case provided above .
the finite - size effects are reasonably small .
for the small sizes we used both eq .
( [ 56 ] ) and our method to compute the dispersion .
the results of the two agree with a relative error of less than @xmath248 .
one important aspect to notice is the small range of the spin - wave spectrum , as compared to the afm exchange @xmath37 mev .
this is a consequence of the fact that the mn spins do not interact directly with one another . instead
, their interaction is mediated by the rather small concentration of charge carriers present .
we compute the density of states ( dos ) @xmath249 associated with the spin - wave dispersion for the superlattice case employing the standard method of dividing the brillouin zone into tetrahedra and linearizing the dispersion ( ref . ) .
the dos obtained for the lattice with @xmath250 , @xmath251 is shown as a dotted line in fig .
we use a logarithmic scale for the energy , and normalize the density of states such that @xmath252 where @xmath253 .
this convention will be maintained throughout the rest of the paper .
upper panel : average density of states @xmath254 _ on a logarithmic scale _ for systems with @xmath255 and 512 mn spins in moderately disordered configurations ( full lines ) .
lower panel : average density of states for systems with @xmath255 and 512 mn spins in strongly disordered configurations ( full lines ) .
the dotted line is the spin - wave density of states of a dms with fully ordered ( superlattice ) configuration of mn spins .
all samples correspond to @xmath256 , @xmath40 . ]    in the upper panel of fig .
[ fig7 ] we also plot the densities of states obtained for moderately disordered configurations with @xmath257 and @xmath258 mn spins ( full lines ) .
these are configurations in which we place the mn spins randomly on the fcc ga sublattice of the host semiconductor , subject to the constraint that the distance between any two spins is _ larger than @xmath259_. [ in the ordered cubic superlattice , nearest neighbor spin separation is @xmath245 ] .
this moderate amount of disorder breaks translational invariance , and the wave - vectors are no longer good quantum numbers .
also , the large degeneracies of the superlattice spectrum are lifted .
we computed the spin - wave spectrum for @xmath261 realizations of the disorder with @xmath262 spins , 100 realizations with @xmath263 spins and 50 realizations with @xmath234 spins .
thus , we have a total of over 21,000 spin - wave energies for each size , and statistics can be comfortably carried out . in particular , from fig .
[ fig7 ] we see that the dos histograms corresponding to the three sizes are smooth functions , i.e. the average over disorder is well accounted for .
also , the curves are practically indistinguishable from one another , implying that finite size effects are negligible .    in the lower panel of fig .
[ fig7 ] we plot the dos corresponding to fully disordered configurations with @xmath257 and @xmath258 mn spins ( full lines ) . these are configurations in which we place the mn spins randomly on the fcc ga sublattice of the host semiconductor , with no restrictions ( except that each mn occupies a different site ) .
the number of samples analyzed is the same as in the previous case .
again , the curves show smooth behavior and no finite - size effects . the origin of the peak appearing near @xmath264mev will be discussed later .    before continuing the analysis
, we must also point out that the goldstone modes have been left out in the dos shown in fig .
the rpa system _ always _ has one solution of energy @xmath150 ( numerically , we find its magnitude to be less that @xmath265 ) , corresponding to @xmath266 .
this is the goldstone mode , describing the same overall rotation of all the spins .
these goldstone modes would appear in the dos as a @xmath267 function at @xmath150 ( @xmath268 ) .
the effect of disorder in the positions of the mn is reflected in the considerable widening of the dos ( on a logarithmic scale ) , and rounding of the van hove singularities of the superlattice dos , as the amount of disorder increases .
more significantly , however , is the _ substantial _ enhancement of the dos at _ low energies_. this behavior is in agreement with the general expectation of the effects of disorder on an energy - band dispersion .    a question of considerable interest concerns the _ nature _ of the spin - wave excitations .
we know that the ordered superlattice has extended , plane - wave type spin - waves .
the general expectation is that disorder will lead to localization , and thus one expects the appearance of localized spin - waves as the disorder increases .
one way to characterize the nature of the spin - waves is to compute their inverse participation ratio ( ipr ) , defined as : @xmath269 for an extended mode @xmath101 ( plane - wave , for instance ) , we expect that all @xmath270 coefficients are roughly of equal magnitude , since all spins are expected to participate equally in the spin - wave .
then , it follows that for an _ extended mode @xmath101 _ , ipr@xmath271 , i.e. it is inversely proportional to the size of the system . on the other hand , in a localized mode @xmath101 , only the spins within the localization volume
have non - vanishing values for @xmath270 .
thus , it follows that for a _
localized mode @xmath101 _ , ipr@xmath272 is independent of @xmath6 . instead
, its value is given by the inverse number of mn spins participating in the localized mode .
upper panel : average ipr(e ) for systems with @xmath255 and 512 mn spins in moderately disordered configurations ( circles , square and triangles , respectively ) .
lower panel : same for strongly disordered configurations .
all samples correspond to @xmath256 , @xmath40 . ]
we compute the ipr of the spin - waves in the following manner .
as we sweep the frequencies @xmath152 of interest and compute the eigenvalues @xmath154 of the @xmath151 matrix on the grid @xmath157 , we also compute the corresponding ipr of the eigenvectors @xmath273 ( which are normalized to unity ) .
the ipr is a continuous function of @xmath152 , and we use linear interpolation to find its value at the spin - wave frequencies @xmath134 of interest . whenever eigenvalues @xmath154 cross , the iprs are no longer well defined : one can choose any linear combination of the eigenvectors of the degenerate modes , which would lead to different values for the corresponding iprs .
we have checked that these accidental crossings are rather rare ( below a few percent ) and equally distributed over the entire energy scale . if we exclude all these degenerate cases , we get a dos which is virtually indistinguishable from the one obtained when these modes are kept .
more importantly , we find that modes which cross predominantly have the same nature ( either localized or extended ) . as a result , the ensemble averaged ipr is not sensitive to the precise treatment of these cases .    in fig .
[ fig8 ] we plot the ( geometric ) average ipr(e ) for the moderate ( upper panel ) and strong ( lower panel ) disorder configurations .
again , the goldstone modes are not shown .
we use the geometric mean ( i.e. , arithmetic mean of the @xmath274 values ) in order to insure proper weight for the extended modes , with low ipr .
for moderate disorder configurations , we see that the spin - wave modes at high energies @xmath275 mev are localized : the values corresponding to the three different sizes collapse on top of one another .
this is expected , given the fact that the upper band - edge of the superlattice spectra is just below @xmath276 mev ( see fig . [ fig7 ] ) .
thus , states of higher energy have been split off the band by disorder , and are expected to be localized .
the low - energy part of the band also appears to be localized : while the three curves do not quite meet , the value of the ipr is just below unity , showing that these spin - waves have considerable weight on only one spin ( i.e. , they correspond to a individual spin - flip ) .
however , for the central part of the spectrum , the spin - waves are extended , with the iprs for the three sizes clearly distinct and decreasing with increasing @xmath6 .
histogram of the ipr values for all spin - waves ( spin flips ) of energy @xmath277 , for @xmath278 and 512 .
results are shown only for strong disorder configurations .
except for the low - lying goldstone modes , the histograms are identical , suggesting only localized modes at these energies . ]    for the strong disorder configurations , this tendency is even more apparent .
the high- and low - energy regions , @xmath279 1 mev , respectively @xmath280 mev , contain localized spin - waves : in both limits , the ipr curves collapse on top of one another .
the low - energy localized modes are individual spin - flips .
the associated @xmath111 are vanishingly small at all sites except one , leading to ipr @xmath281 .
these sites are always situated far apart from all other mn spins , and the probability to be visited by charge carriers is exponentially small ( tailing from far occupied regions ) . as a result ,
the mn spins at these sites are virtually isolated , and their spin excitations are individual spin - flips .
the energy for such a spin - flip is equal to @xmath53 , if one neglects small corrections due to the extremely weak interactions [ see eq .
( [ a.9b ] ) and following discussion ] .
histograms of the ipr of all the modes with energies below @xmath282 are shown in fig .
[ fig9 ] . here
, we also show the goldstone modes , which have zero energy , and @xmath283 .
the histograms have been scaled by the total number of modes for all disorder realizations considered for each particular size . since the number of goldstone modes exactly equals the number of different realizations of disorder considered , their peaks are in a ratio of roughly 200/100/50 with respect to one another .
the main observation is that the histograms for the three sizes are very similar , with a huge peak just below ipr=1 .
the absence of a dependence on @xmath6 confirms that all these modes are localized ( except the goldstone modes , of course ) .
histogram of the ipr values for all spin - waves ( spin flips ) of energy @xmath284 , for @xmath278 and 512 .
results are shown only for strong disorder configurations .
the histograms are almost identical , suggesting only localized modes at these energies . ]
the high - energy localized states are of very different nature : their ipr is close to @xmath285 , suggesting that it has large contributions from only two sites .
analysis of the @xmath111 values for such modes shows that most of them are associated with nearest - neighbor ( n.n . )
mn spins on the fcc ga sublattice .
the characteristic energy for the spin - waves centered on such n.n .
sites is just below 5 mev , and this is the feature responsible for the peak appearing in the dos at these energies ( see fig .
[ fig7 ] , lower panel ) .
careful inspection of the dos reveals the appearance of smaller peaks at somewhat lower energies ( even the dos for moderate disorder has a small peak at around 2 mev , and its ipr for these modes is close to 0.5 ) .
these peaks are associated with excitations of spin - pairs ( or larger clusters ) with varying separations , but they are not as well defined as the one corresponding to nearest - neighbors spins .    the small clusters of mn spins giving rise to such modes are always found to be in regions densely populated with charge carriers . due to their closeness in space , these spins are much more strongly coupled to one another than they are to other spins with which they share charge carriers .
this leads to the `` resonance - like '' character of these modes : either of the spins can be flipped with equal probability . as a result of the strong - coupling , the cost of flipping either spin is high , since it frustrates their fm arrangement .
indeed , the characteristic energy of roughly @xmath286 reflects much stronger coupling than the average one present in the superlattice case .
the histogram of the ipr of all modes of energy @xmath284 shown in fig .
[ fig10 ] confirms all these conclusions .
we caution that these energies may be substantially modified due to direct antiferromagnetic mn - mn exchange ( left out in our model ) in real systems .
histogram of the ipr values for all spin - waves ( spin flips ) of energy @xmath287 , for @xmath288 ( upper ) , 216 ( middle ) and 512 ( lower panel ) .
results are shown only for strong disorder configurations . with increasing @xmath6 the histograms shift to lower values , suggesting that these spin - wave modes are extended . ]    finally , the spin - waves at intermediate energies @xmath289 correspond to extended modes .
it is rather difficult to establish exactly the `` mobility edge '' corresponding to the transitions from localized to extended modes at either end , given the rather small system sizes considered here .
thus , we use the values quoted above only as plausible estimates of these boundaries . clearly , for these energies the average ipr decreases with increasing @xmath6 ( see fig .
[ fig8 ] ) . in fig .
[ fig11 ] we plot the histograms of the ipr for all the modes with energy @xmath290 . in this case ,
size dependence is clearly visible , both in the position of the broad peak maximum as well as in the lower cut - off in ipr value .
the position of the maxima scale roughly like @xmath291 .
we have attempted both gaussian and lorentzian fits , but neither seems to capture the low - ipr tail properly . in the smaller samples , a second peak appears just below ipr=1 . at first sight , this might seem to indicate the existence of a finite density of localized states at these energies as well . however , we believe that this peak is a finite - size effect
. clearly , its amplitude decreases with increasing @xmath6 and would vanish in the thermodynamic limit .
this is very different from the behavior observed in either of the two localized energy ranges , where the distributions for all sizes are virtually identical ( see figs .
[ fig9 ] and [ fig10 ] ) .
the reason for the appearance of this second peak is simply the restriction ipr @xmath292 1 . as the broad peak moves towards higher values with decreasing @xmath6 ,
all the values in the upper tail `` bunch '' at ipr=1 .
one final important observation relates to the absolute values of the iprs in the extended spin - wave regime . although scaling with system size is present , the corresponding ipr values are much higher than the ones expected for modes extended over _ all _ mn sites ( these values are shown by the positions of the goldstone modes in fig .
[ fig9 ] , and are well below the cut - offs observed in fig .
[ fig11 ] ) .
this means that in the disordered system , the extended spin - waves are delocalized only over a fraction of the total number of sites . as the amount of disorder increases , this fraction decreases , since ipr averages for the moderate disorder are smaller than ipr averages for full disorder .
these results are consistent with the picture provided by the temperature dependent mean - field study of this model ( ref . ) .
there , we concluded that in the disordered system , the ( small fraction @xmath195 ) of charge carriers are concentrated in a small volume of the sample , where the local mn concentration is larger than average . as a result ,
the concentration of charge carriers is strongly enhanced in these regions , and the exchange with the mn spins inside these regions ( which is proportional to the probability of finding charge carriers nearby ) is greatly increased , preserving magnetization of these regions up to high temperatures . on the other hand ,
the mn spins in the regions devoid of charge carriers are very weakly coupled , and behave as free spins down to very low temperatures .
due to the very inhomogeneous assembly of weakly- and strongly- interacting spins , the magnetization curves have very unusual , concave shapes .
the present study of the spin - waves corroborates the same picture .
we find very low - energy , spin - flip like excitations , which are obviously due to the weakly interacting spins , and which are responsible for a very sharp decrease of the magnetization at exponentially small temperatures.@xcite this is is to be contrasted with the behavior of ordered conventional ferromagnets , where at low temperatures only low - energy long - wavelength spin - waves can be excited .
since their phase space is vanishing in the long - wavelength limit ( @xmath293 in @xmath294 dimensions ) , the magnetization of conventional ferromagnets decreases very slowly from its saturation value with increasing temperature , leading to convex upwards magnetization curves.@xcite    in our model , the extended spin - waves are concentrated around the high - density mn regions , where the charge carriers mediating the interactions are to be found .
this is consistent with the appearance of modes whose ipr , while scaling with the system size , shows modes extended only over a fraction of all system sites .
as the amount of disorder decreases all the way to a fully ordered superlattice , the charge carriers are more and more homogeneously spread throughout the sample and the ipr of the extended modes has lower and lower values , as observed from fig .
however , the more homogeneous a sample is , the less the average afm coupling between the mn spins and the charge carriers , since the charge carriers have now equal probability of being found anywhere in the sample , instead of being concentrated in a small fraction of the space .
this leads to a decrease of the critical temperature @xmath295 , as observed in both mean - field @xcite and monte - carlo@xcite analysis .
the enhancement of the mean - field @xmath295 with increasing disorder is suggested also by the existence of the high - energy localized modes , which show that ferromagnetic alignment will persist in high - density clusters up to very high temperatures .
it also suggests that local ferromagnetic fluctuations might be observed well above @xmath295 .
the aim of this paper is two - fold .
first , we demonstrate the accuracy and speed of a new scheme for computing rpa spectra , which allows tackling of systems of much larger sizes than the ones that can be analyzed with the standard rpa formulation .
investigation of large system sizes and averages over many disorder realizations facilitate a clear picture for the problem we are interested in , namely the spectrum and nature of spin - waves of a disordered dms .
we then demonstrate that disorder can significantly change the spectrum and the nature of the spin - waves .
this is likely to lead to important consequences not only as far as magnetic properties are concerned ( we have already commented on the fast de - magnetization with increasing temperature , due to low - energy spin - flip modes ) .
more importantly , this may have significant consequences for transport properties as well .
for instance , charge carrier spin scattering is likely to be very different in various regions of the sample .
large anomalous hall effect has been observed in ( ga , mn)as,@xcite but the theory used to interpret it is borrowed from phenomenology relevant to homogeneous ferromagnetic metals . in an inhomogeneous system
, some of the accepted ideas might have to be modified or at least verified to still hold true .
this research was supported by nsf dmr-9809483 .
m.b . was supported in part by a postdoctoral fellowship from the natural sciences and engineering research council of canada .
t. dietl , a. haury and y. m. daubign , phys .
b * 55 * , r3347 ( 1997 ) ; m. takahashi , phys . rev .
b * 56 * , 7389 ( 1997 ) ; t. jungwirth , w. a. atkinson , b. h. lee and a. h. macdonald , phys . rev .
b * 59 * , 9818 ( 1999 ) ; t. dietl , h. ohno , f. matsukura , j. cibert and d. ferrand , science * 287 * , 1019 ( 2000 ) .
uniform magnets as well as metallic alloys are characterized by weak deviations from a universal curve on a @xmath296 versus @xmath297 and @xmath298 plot .
[ see t. kaneyoshi , _ introduction to amorphous magnets _ , p 56 - 61 ( world scientific , 1992 ) ] .