[ ]{}

[*Introduction.*]{} The recent discovery of carrier-induced ferromagnetism[@Story86; @Ohno92/96] in diluted magnetic semiconductors (DMS) has generated intense interest, in part because it opens the prospect of developing devices which combine information processing and storage functionalities in one material [@Ohno99; @Prinz98; @Hayashi98; @Pekarek98; @VanEsch97; @Oiwa99; @Beschoten99]. Critical temperatures $T_c$ exceeding $100 {\rm K}$ have been realized[@Ohno99] by using low-temperature molecular beam epitaxy growth to introduce a high concentration $c$ of randomly distributed ${\rm Mn}^{2+}$ ions in GaAs systems with a high hole density $c^*$. The tendency toward ferromagnetism and trends in the observed $T_c$’s have been explained using a picture[@Dietl97; @Takahashi97; @Jungwirth99; @Dietl00; @Lee00] in which uniform itinerant-carrier spin polarization mediates a long-range ferromagnetic interaction between the ${\rm Mn}^{2+}$ ions with spin $S=5/2$.

We present here the first theory which accounts for dynamic correlations in the ordered state and is able to describe its fundamental properties. We use a path-integral formulation and, in the RKKY spirit of the previous theory, integrate out the itinerant carriers, and expand the effective action for the impurity spins up to quadratic order. In addition to the usual Goldstone spin waves, we find itinerant-carrier dominated collective excitations analogous to the optical spin-waves in a ferrimagnet. We also find that the spin stiffness is inversely proportional to the itinerant-carrier mass and independent of the strength of the exchange coupling between magnetic ions and itinerant carriers, in contrast to the mean-field critical temperature which is proportional to the itinerant-carrier mass and to the square of the exchange coupling. We extend our calculations to the critical temperature using an approximate self-consistent spin-wave scheme and address the temperature-dependent magnetization and specific heat.

[*Hamiltonian and effective action.*]{} We study a model which provides an accurate[@modelrefs] description of Mn based zincblende DMS’s. Magnetic ions with $S=5/2$ at positions $\vec R_I$ are antiferromagnetically coupled to valence-band carriers described by envelope functions, $$H = H_0 + J_{pd} \int \hspace{-1mm}d^3 r \;
  \vec S(\vec r) \cdot \vec s(\vec r) ,$$ where $\vec S(\vec r) = \sum_I \vec S_I \delta( \vec r - \vec R_I)$ is the impurity-spin density. The itinerant-carrier spin density is expressed in terms of carrier field operators by $\vec s(\vec r)= {1\over 2} \sum_{\sigma \sigma'} \Psi^\dagger_\sigma (\vec r)
 \vec \tau_{\sigma \sigma'} \Psi_{\sigma'} (\vec r)$. ($\vec \tau$ is the vector of Pauli spin matrices.) $H_0$ includes the valence-band envelope-function Hamiltonian[@band] and, if an external magnetic field $\vec B$ is present, the Zeeman energy, $$\begin{aligned}
  H_0 &=& \int \hspace{-1mm}d^3 r\; \bigg\{ \sum_\sigma
      \hat\Psi^\dagger_\sigma (\vec r)
      \left( -{\hbar^2 \vec\nabla^2\over 2m^*}-\mu^{*} \right)
      \hat \Psi_\sigma (\vec r)
\nonumber \\
  &&- g\mu_B \vec B \cdot \vec S(\vec r)
 - g^* \mu_B \vec B \cdot \vec s (\vec r)\bigg\} .\end{aligned}$$ The effective mass, chemical potential, and $g$-factor of the itinerant carriers are labeled by $m^*$, $\mu^*$, and $g^*$. The model we use here is, thus, related to colossol magnetoresistance (CMR) materials[@Millis95] and identical to those for dense Kondo systems, which simplify when the itinerant-carrier density $c^*$ is much smaller than the magnetic ion density $c$[@Sigrist91]. The fact that $c^*/c \ll 1$ in ferromagnetic-semiconductor materials is essential to their ferromagnetism. Similar models have been used for ferromagnetism induced by magnetic ions in nearly ferromagnetic metals such as palladium[@Doniach67].

We represent the impurity spins in terms of Holstein-Primakoff (HP) bosons[@Auerbach94]. By coarse graining, the spin density $\vec S(\vec r)$ can be replaced by a smooth function, $$\begin{aligned}
   S^+(\vec r) &=& \left(\sqrt{2cS - b^{\dag}(\vec r) b(\vec r)} \, \right)
   b(\vec r)
\\
   S^-(\vec r) &=& b^\dag(\vec r) \sqrt{2cS - b^{\dag}(\vec r) b(\vec r)}
\\
   S^z(\vec r) &=& cS - b^{\dag}(\vec r) b(\vec r)\end{aligned}$$ with bosonic fields $b^{\dag}(\vec r), b(\vec r)$. The partition function as a coherent-state path-integral in imaginary times reads, $$\label{partition function}
 Z = \int \hspace{-1mm}
 {\cal D} [\bar z z] {\cal D} [\bar\Psi \Psi]
 e^{-\int_0^\beta d\tau L( \bar z z, \bar\Psi \Psi)}$$ with $L = \int d^3 r \left[ \bar z \partial_\tau z +
\sum_\sigma \bar \Psi_\sigma \partial_\tau \Psi_\sigma \right] +
H (\bar z z, \bar\Psi \Psi)$. The bosonic (impurity spins) and fermionic (itinerant carriers) degrees of freedom are labeled by the complex variables $\bar z,z$ and the Grassmann numbers $\bar \Psi, \Psi$, respectively.

Since the Hamiltonian is bilinear in fermionic fields, we can integrate out the itinerant carriers and arrive at an effective description in terms of the localized spin density only, $Z = \int {\cal D} [\bar z z] 
\exp (- S_{\rm eff} [\bar z z])$ with the action $$\begin{aligned}
  S_{\rm eff} [\bar z z] = \int_0^\beta \hspace{-1mm}d\tau
  \int \hspace{-1mm}d^3 r \left[
  \bar z \partial_\tau z - g\mu_B B (cS - \bar zz)
  \right]
\nonumber \\
  - \ln \det \left[ (G^{MF})^{-1} + \delta G^{-1}(\bar zz) \right] \, .
\label{effective action}\end{aligned}$$ Here, we have already split the total kernel $G^{-1}$ into a mean-field part $(G^{MF})^{-1}$ and a fluctuating part $\delta G^{-1}$, $$\begin{aligned}
   (G^{MF})^{-1} &=&
 \left( \partial_\tau - {\hbar^2 \vec \nabla^2 \over 2m^*} -\mu^* \right)
 {\bf 1} + {\Delta \over 2}\tau^z
\\
   \delta G^{-1} &=& {J_{pd}\over 2} \left[
    (z \tau^- + \bar z \tau^+) \sqrt{2cS - \bar z z}  - \bar z z \tau^z
\right]\end{aligned}$$ where $\Delta = cJ_{pd}S - g^*\mu_B B$ is the zero-temperature spin-splitting gap for the itinerant carriers. The physics of the itinerant carriers is embedded in the effective action of the magnetic ions. It is responsible for the retarded and non-local character of the interactions between magnetic ions, which is described here for the first time.

[*Independent spin-wave theory.*]{} The independent spin-wave theory, which is a good approximation at low temperatures, is obtained by expanding Eq. (\[effective action\]) up to quadratic order in $z$. We find (in Fourier representation) $$S_{\rm eff} [\bar z z] = {1\over \beta V} \hspace{-1mm}
  \sum_{|\vec p| < p_c, m} \! \!
  \bar z(\vec p, \nu_m) D^{-1}(\vec p, \nu_m) z(\vec p, \nu_m),$$ in addition to the temperature-dependent mean-field contribution. A Debye cutoff ($p_c^3 = 6 \pi^2 c$) ensures that we include the correct number of magnetic ion degrees of freedom. The kernel of the quadratic action is the inverse of the spin-wave propagator, $$\begin{aligned}
   D^{-1}(\vec p, \nu_m) = 
 -i \nu_m + g\mu_B B + J_{pd} n^{*}
 \hspace*{2.cm}
\nonumber \\
 + { c J_{pd}^2 S \over 2\beta V} \sum_{n, \vec k}
 G^{MF}_\uparrow (\vec k, \omega_n)
 G^{MF}_\downarrow (\vec k+ \vec p, \omega_n+\nu_m)
\label{inverse propagator}\end{aligned}$$ with the mean-field itinerant carrier Green’s function $G^{MF}_\sigma (\vec k, \omega_n) = - \left[ i\omega_n - \left(
\epsilon_{\vec k} +\sigma \Delta /2 -\mu^* \right) \right]^{-1}$. The mean-field spin density is denoted by $n^{*}=(n_\downarrow - n_\uparrow )/2$, and $\epsilon_{\vec k} = \hbar^2 k^2/(2 m^*)$.

[*Excitations.*]{} We obtain the spectral density of the spin-fluctuation propagator by analytical continuation, $i\nu_m \rightarrow \Omega+i0^+$ and $A(\vec p,\Omega)={\rm Im} \, D(\vec p,\Omega)/\pi$. We find three different types of spin excitations. In all figures we take $B=0$ and use as typical parameters[@Ohno99] $m^*=0.5m_e$, $J_{pd}=0.15{\rm eVnm}^3$, and $c=1{\rm nm}^{-3}$, where $m_e$ is the free-electron mass. For these parameters the mean-field itinerant-carrier system is fully polarized at $T=0$.

i\) Our model has a gapless Goldstone-mode branch (see Fig. \[fig1\]) reflecting the spontaneous breaking of rotational symmetry[@gap]. Expansion of the $T=0$ propagator at $\Delta > \epsilon_F$, where $\epsilon_F$ is the Fermi energy of the majority-spin band, yields for small and large momenta the dispersion of the collective modes $$\Omega^{(1)}_p = { x\over 1- x} \epsilon_p \left( 1 -
 {4 \epsilon_F\over 5\Delta} \right) + {\cal O}(p^4)$$ and $\Omega^{(1)}_p = x \Delta ( 1 - \Delta / \epsilon_p ) + {\cal O} (1/p^4)$. At short wavelengths we obtain the mean-field result $x\Delta$, the spin-splitting of a magnetic ion in the effective field produced by fully spin-polarized itinerant carriers. Note that the itinerant-carrier and magnetic-ion mean-field spin splittings differ by a factor of $x=c^*/(2cS) \ll 1$. At long wavelengths the magnon dispersion in an isotropic ferromagnet is proportional the spin stiffness $\rho$ divided by the magnetization $M$. In the adiabatic limit[@Millis95; @Nagaev98], $\epsilon_F\ll\Delta$, our long-wavelength result reflects a spin-stiffness due entirely to the increase in kinetic energy of a fully spin-polarized band when the orientation has a spatial dependence, $\rho = c^* \hbar^2/(4m^*)$, and a magnetization $M$ which has opposing contributions from magnetic ions and itinerant carriers, $M = c S -c^*/2 = c S (1-x)$. In this limit, the mean-field critical temperature and the spin stiffness have opposite dependences on the itinerant-carrier mass.

The low-energy mode describes spin waves in the local-impurity system. This contrasts with the case of ferromagnetism induced by local moments[@Doniach67] where the low-energy mode is primarily in the free-carrier part.

![Spin-wave dispersion for $c^*= 0.1 {\rm nm}^{-3}$. The short wavelength limit is the mean-field result $x\Delta$. For comparison, we show also the result obtained from an RKKY picture.[]{data-label="fig1"}](fig1.eps){width="8cm"}

ii\) We find a continuum of Stoner spin-flip particle-hole excitations. They correspond to flipping a single spin in the itinerant-carrier system and, therefore, occur at much larger energies $\sim \Delta$ (see Fig. \[fig2\]). For $\Delta > \epsilon_F$ and zero temperature, all these excitations carry spin $S^z=+1$, i.e., increase the spin polarization, and therefore turn up at negative frequencies in the boson propagator we study. (When $\Delta < \epsilon_F$, excitations with both $S^z=+1$ and $S^z=-1$ contribute to the spectral function.) This continuum lies between the curves $-\Delta-\epsilon_p\pm 
2\sqrt{\epsilon_p\epsilon_{F}}$ and for $\Delta < \epsilon_F$ also between $-\Delta+\epsilon_p\pm 2\sqrt{\epsilon_p(\epsilon_F-\Delta)}$.

![Stoner excitations and collective modes in the free-carrier system for $c^*=0.1 {\rm nm}^{-3}$.[]{data-label="fig2"}](fig2.eps){width="8cm"}

iii\) We find additional collective modes associated primarily with the itinerant-carrier system at energies below the Stoner continuum (see Fig. \[fig2\]). At $T=0$ we obtain $$-\Omega^{(2)}_p = \Delta (1-x) - { 1\over 1- x}
    \epsilon_p \left( {4 \epsilon_F\over 5x\Delta} -1 \right)
    + {\cal O}(p^4)$$ for $\Delta>\epsilon_F$. The spectral weight of these modes is $-x/(1-x)$ at zero momentum.

The finite spectral weight at negative energies indicates that the ground state is not fully spin polarized because of quantum fluctuations.

[*Comparison to RKKY picture.*]{} For comparison we evaluate the $T=0$ magnon dispersion assuming an RKKY interaction between magnetic ions. This approximation results from our theory if we neglect the spin polarization in the itinerant carriers and evaluate the static limit of the resulting spin-wave propagator. The Stoner excitations and magnons shown in Fig. \[fig2\] do not emerge and the spin-wave dispersion (Fig. \[fig1\]) is incorrect except for the limit $\Delta\ll\epsilon_F$. While the RKKY picture does, in some circumstances, provide a realistic estimate of $T_c$, it completely fails as a theory of the ferromagnetic state.

[*Comparison to a ferrimagnet.*]{} Some features of the excitation spectrum are, not coincidentally, like those of a localized spin ferrimagnet with antiferromagnetically coupled large spin $S$ (corresponding to the magnetic ions) and small spin $s$ (corresponding to the free carriers) subsystems on a bipartite cubic lattice. Representing the spins by HP bosons and expanding up to quadratic order, we either diagonalize the resulting Hamiltonian directly or, as in our ferromagnetic-semiconductor calculations, integrate out the smaller spins using a path-integral formulation to obtain equivalent results. We find that $$\Omega^{(1)/(2)}_{\vec p} = {\Delta\over 2} \left[ -(1-x) \pm
 \sqrt{(1-x)^2+4x\gamma_{\vec p}} \, \right]$$ with $x=s/S$ and $\Delta=6JS$, where $J$ is the exchange coupling, and $\gamma_{\vec p}=(1/3)\sum_{i}[1-\cos(p_{i}a)]$ with lattice constant $a$. We recover two collective modes, the coupled spin waves of the two subsystems. One is gapless, the other one gapped with $\Delta(1-x)$. The bandwidth is $x\Delta$, and the spectral weights at zero momentum are $1/(1-x)$ and $-x/(1-x)$, respectively, as in our model.

[*Self-consistent scheme.*]{} Near the transition temperature $T_c$, the spin-wave density is of the order of $cS$, and expanding around a fully-polarized state is not a good starting point. Instead, we adopt a scheme which accounts self-consistently for the reduction of spin density in the localized impurity and itinerant carrier subsystems[@foot]. Our approach is motivated by the Weiss mean-field theory, which is expected to be accurate for a model with static long-range interactions between the magnetic ions. This model can be obtained in our approach by taking the Ising limit (i.e., replacing $\vec S \cdot \vec s$ by $S^z s^z$) and letting the mean-field itinerant-carrier spin splitting, $\Delta(T)= J_{pd} \langle S^z \rangle$, reflect the thermal suppression of the impurity-spin density $\langle S^z \rangle$. In this simple model the spin-wave spectrum is independent of momentum, $\Omega_{p}(T) = J_{pd} n^{*}(T)$, allowing a unitary transformation to spin-wave eigenstates with a site label. The constraint on the number of spin bosons ($ \leq 2S$) on a site is then easily applied and we obtain $$\begin{aligned}
   &&\langle S^z \rangle = {1\over V}\sum_{|\vec p| < p_c}
 S B_S (\beta S \Omega_p)
\nonumber\\
   && = {1\over V}\sum_{|\vec p| < p_c} \left\{ S - n(\Omega_p)
 + (2S+1) n[(2S+1)\Omega_p] \right\},
\label{SelfConsistent}\end{aligned}$$ where $B_S(x)$ is the Brillouin function, and $n(\omega)$ is the Bose function. The first two terms in the second form of Eq. (\[SelfConsistent\]) give the spin density obtained by treating spin waves as bosons, while the last term is the correction from spin kinematics which rules out unphysical states.

Our self-consistent spin-wave approximation consists of using Eq. (\[SelfConsistent\]) when $\Omega_p$ is $p-$dependent (see Fig. \[fig3\]). Although the characteristic $T^{3/2}$ law for the localized-ions magnetization is recovered at low temperatures, the prefactor is reduced by $1-(2S+1)^{-1/2}$, compared to the correct[@Dyson56] value of linearized spin-wave theory. That is, constraining boson populations in momentum space is too restrictive at low temperatures, although the error is not serious for $S=5/2$. The transition temperature of the self-consistent spin-wave theory is lower than the estimate from Weiss mean-field theory.

The presence of spin waves also shows up in the specific heat $C_V$, which we calculate within the same scheme (see Fig. \[fig4\]). Calculating the entropy $S$ we obtain $$C_V = {T\over V}{dS\over dT} =
         \frac{1}{V} \sum_{|\vec p|<p_{c}} \Omega_{p}(T)
    \left(\frac{dN_{p}}{dT} \right),$$ where $N_{p} = S - S B_{S}(\beta S \Omega_{p})$ is the average number of spin bosons at each momentum. The specific heat from the kinetic energy of itinerant carriers turns out to be negligibly small. The specific heat of magnetic ions $C_{V}$ is proportional to $T^{3/2}$ at low temperature, shows a Schottky-like anomaly, and has a jump at the transition temperature. We estimate that the jump is 5% of the lattice specific heat when $T_c \sim 100{\rm K}$, suggesting that this should be observable.

![Average impurity and itinerant-carrier spin as a function of temperature for $c^*=0.1{\rm nm}^{-3}$ in Weiss mean-field (Weiss) and self-consistent spin-wave theory (SSW). Inset: $T_c$ as a function of $c^*$.[]{data-label="fig3"}](fig3.eps){width="7.9cm"}

![Specific heat due to spin waves.[]{data-label="fig4"}](fig4.eps){width="7.8cm"}

We acknowledge useful discussions with M. Abolfath, B. Beschoten, A. Burkov, J. Furdyna, S. Girvin, and T. Jungwirth. This work was supported by the Deutsche Forschungsgemeinschaft under grant KO 1987-1/1 and by the National Science Foundation, DMR-9714055.

[10]{}

T. Story, R.R. Gałazka, R.B. Frankel, and P.A. Wolff, Phys. Rev. Lett. [**56**]{}, 777 (1986).

H. Ohno [*et al.*]{}, Phys. Rev. Lett. [**68**]{}, 2664 (1992); H. Ohno [*et al.*]{}, Appl. Phys. Lett. [**69**]{}, 363 (1996).

H. Ohno, J. of Magnetism and Magnetic Materials [**200**]{}, 110 (1999); Science [**281**]{}, 951 (1998).

G.A. Prinz, Science [**282**]{}, 1660 (1998).

T. Hayashi [*et al.*]{}, J. Appl. Phys. [**83**]{}, 6551 (1998).

T.M. Pekarek, B.C. Crooker, I. Miotkowski, and A.K. Ramdas, J. Appl. Phys. [**83**]{}, 6557 (1998).

A. Van Esch [*et al.*]{}, J. Phys. Condensed Matter [**9**]{}, L361 (1997); A. Van Esch [*et al.*]{}, Phys. Rev. B [**56**]{}, 13103 (1997).

A. Oiwa [*et al.*]{}, Phys. Rev. B [**59**]{}, 5826 (1999); F. Matsukura, H. Ohno, A. Shen, and Y. Sugawara, Phys. Rev. B [**57**]{}, R2037 (1998).

B. Beschoten [*et al.*]{}, Phys. Rev. Lett. [**83**]{}, 3073 (1999).

T. Dietl, A. Haury, and Y.M. d’Aubign[é]{}, Phys. Rev. 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 [*et al.*]{}, Science [**287**]{}, 1019 (2000).

B.H. Lee, T. Jungwirth, and A.H. MacDonald, cond-mat/0001391, Phys. Rev. B [**61**]{}, June 15 issue (2000).

The central cell corrections required for an accurate envelope-function description of an isolated ${\rm Mn}^{2+}$ ion in GaAs \[A.K. Bhattacharjee and C. Benoit à la Guillaume, Solid State Commun. [**113**]{}, 17 (2000)\] will be reduced by free-carrier screening.

To simplify the present discussion we use a generic single-band model and neglect interactions between the free carriers as well as the antiferromagnetic Mn-Mn interaction. In quantitative calculations for particular materials it is necessary to account for enhanced ferromagnetism due to interactions [@Jungwirth99], realistic heavy and light hole bands, and at high carrier densities, also the spin-orbit split-off band [@Dietl00]. These important effects will be discussed separately \[J. K" onig, H. Lin, and A.H. MacDonald (to be published)\].

A.J. Millis, P.B. Littlewood, and B.I. Shraiman, Phys. Rev. Lett. [**74**]{}, 5144 (1995).

M. Sigrist, H.Tsunetsugu, and K. Ueda, Phys. Rev. Lett. [**67**]{}, 2211 (1991).

S. Doniach and E.P. Wohlfarth, Proc. R. Soc. London Ser. A [**296**]{}, 442 (1967).

A. Auerbach, [*Interacting Electrons and Quantum Magnetism*]{} (Springer, New York, 1994).

Because of spin-orbit coupling in the valence bands, the spin-wave spectrum in realistic models will have a strain- and disorder-dependent gap \[M. Abolfath, J. Brum, T. Jungwirth, and A.H. MacDonald (to be published)\].

E.L. Nagaev, Phys. Rev. B [**58**]{}, 827 (1998).

In the following, we neglect Stoner excitations and magnons in the itinerant carriers, i.e., we assume sharp spin-wave excitations with spectral weight one.

F.J. Dyson, Phys. Rev. [**102**]{}, 1230 (1956).
