Magnetism mandates that the electrons in a material collectively align their spins. In a ferromagnet there is an imbalance in the occupancy of spin-up and spin-down states; such a collective spin alignment leads to a macroscopic spontaneous magnetization that persists up to a Curie temperature . Many common ferromagnets (such as iron or nickel) are metals, exhibiting a set amount of free electrons that make them conduct electrically. In semiconductors, on the other hand, the concentration and type of electronic carriers is controllable externally, for instance, by doping with select atom types or directly by injecting carriers. Most semiconductors are not magnetic. A material that exhibits both ferromagnetic and semiconductor properties offers the exciting prospect of combining nonvolatile magnetic storage and conventional semiconductor electronics in a single device. Magnetic semiconductors offer a number of interesting possibilities in the pursuit of “spintronics,”a branch of science and technology that exploits the spin dimension of the electron in addition to its charge, for novel electronic devices. These materials combine the properties of a semiconductor and a magnetic material, providing, for instance, a way to create spin-polarized currents, and by the same token, the promise of electrical control of magnetic effects. While in some magnetic semiconductors, for example, magnetite, all of the material’s constituent ions are intrinsically magnetic (“concentrated magnets”), the most recent focus has been on nonmagnetic semiconductor host materials that can be doped by a small amount of magnetic transition-metal ions or by defects that promote magnetism (“dilute magnets”).
The first experiments on dilute magnetic semiconductors focused on host semiconductors with narrow band gaps, such as and , which are infrared-absorbing black substances. When doped with manganese, substitute for ions, thus releasing “holes” (a positive charge resulting from the absence of an electron) that are said to collectively align the spins ferromagnetically when the sample is cooled below . Ferromagnetism in cannot, however, be enhanced indefinitely by adding more ions, as some of the additional manganese will occupy interstitial sites or form some clustered phases that inhibit a global ferromagnetic order. This has limited the ferromagnetism in and other doped III-V semiconductors to experimentally relevant, but not technologically useful, temperatures. For this reason, the possibility of doping common insulators like , and instead was seized upon with great enthusiasm. These materials have the additional bonus feature of a wide band gap, promising the control of magnetism through charge carriers in an optically transparent medium. Reports on ferromagnetism above room temperature in such wide-gap materials quickly appeared—even in samples doped with elements that were not in themselves magnetic, or not doped at all! The excitement brought back fond memories of the early days of room-temperature superconductivity, when unusual phenomena were discovered in a group of unsuspected oxides and nitrides. What was special about these developments in pushing the boundaries in dilute magnets was the role of theory, which in many cases preceded experiments and influenced the type of systems that were studied. The nature of the delicate tango danced by theory and experiment in this field is the subject of this article.
The tools of solid-state theory of magnetic systems have played a central role in this rapidly developing field, offering not only ex post facto explanations of measured phenomena, but very often leading the field by predicting specific combinations of host insulators with impurity atoms that would lead to room-temperature ferromagnetism. In this pursuit, two proactive cultures of predictive theories have emerged. The initial basic philosophy was to use model Hamiltonians, in which certain specific magnetic interactions between ions are postulated (and others are excluded). For this technique, the “user” needs to guess at the outset the type of magnetic interactions that will dominate. This style of model-Hamiltonian theory has had a glorious past in numerous areas of condensed matter physics. The most prominent supposition for the interaction type between spin-polarized impurity ions (such as , , and other common dopants from the transition-metal series) in insulators assumed it was similar to the interaction of nuclear spins mediated by the conduction electrons in metals, and the model also assumes that the host crystal is largely unperturbed by the presence of the impurity ion. This type of Ruderman-Kittel-Kasuya-Yosida (RKKY)  model-Hamiltonian approach  went all the way to predict ferromagnetism in various impurity-doped compound semiconductors. In particular, the prediction that ferromagnetism persists to a high Curie temperature with ’s well above room temperature in various wide-gap oxides and nitrides —has resulted in thousands of papers in the past decade, trying to capitalize on these exciting predictions offering high in common insulators.
Alongside the culture of model-Hamiltonian theories of dilute magnets, first-principles calculations have rapidly penetrated this field. The philosophy here is to describe each host and impurity system by considering its fundamental electronic structure, thereby obtaining the magnetic interaction type as the answer, not the question. Instead of prejudging the preferred type of magnetic interactions between ions, one articulates fundamental electron-ion and electron-electron interactions and lets the magnetic interactions emerge as solutions to the fundamental many-electron Hamiltonian. The availability of density-functional approximations (exchange and correlation functionals for the inter-electronic interactions), accurate pseudopotentials (simplifying the calculation greatly), and the ubiquitous computer packages encoding such developments into friendly interfaces freed the theorist from having to guess at the outset how a particular impurity would affect a particular host crystal. At the same time, this apparent freedom sometimes created the illusion that such modeling can be done on “automatic pilot,” disregarding potential pitfalls. Regrettably, such pitfalls sometimes created “false positive” results (suggesting high ’s when no magnetic ordering should occur), spurring enthusiasm and optimism about high-temperature ferromagnetism even when not warranted. Indeed, it has often been the case that materials synthesized from computational recipes have not behaved in the way that was expected.
In the following we will (i) discuss some of the issues and potential pitfalls that can come with uncritical application of theoretical methodology, (ii) analyze the physics of the artifact, and (iii) suggest specific ways out of these dilemmas. The good news is that the insights gained from such critical evaluation can now be embodied (in different ways!) into practical first-principles calculations, offering a truly predictive design option. Furthermore, such calculations clarify what types of interactions are indeed critical and offer model-Hamiltonian approaches a safer, albeit ex-post-facto, method to construct more realistic Hamiltonians.
Transition-metal impurities in semiconductors: Interesting physics, no ferromagnetism
The electronic properties of isolated transition-metal impurities in semiconductors were studied in great detail in the 1980s [3, 4], when the growth methods produced ultradilute samples (defined as containing – impurities/) with virtually no impurity-impurity interaction. The high crystalline perfection of such samples enabled a clear experimental and theoretical understanding of the properties and chemical trends in the periodic table of isolated impurities, including their spin configuration, donor and acceptor transitions, and optical excitations. A number of intriguing physical phenomena were revealed due to coexisting localized impurity states. One of these is the “self-regulating response” [5, 6], recently revisited by Raebiger et al. , which shows that adding electrons to the derived impurity levels does not change the integrated electron density around the impurity atom because the orbitals of the surrounding ligand atoms rehybridize to minimize the impact of impurity charging. A second phenomenon is the “vacuum pinning rule,” whereby transition levels of the same impurity atom in different host semiconductors line up in approximately the same way, regardless of the host chosen [8, 9], thus permitting the deduction of band offsets between different pure materials from the knowledge of impurity level positions . Third, another consequence is called “exchange-correlation negative ,”  whereby the strongly varying exchange splitting of differently charged impurities can lead to effective electron-electron attraction. These systems all show interesting physics, however, such systems revealed no ferromagnetism, as the generally low solubility of dopants precluded collective interactions.
Collective magnetic ordering requires the individual magnetic impurity ions to be close enough to one another so as to interact. In practice, this implies large concentrations, often well above typical thermodynamic solubility limits. The low solubility limit of most ions in semiconductors can, however, be overcome by using nonequilibrium growth techniques (such as very low-temperature molecular beam epitaxy). These techniques increased the impurity concentration by a few orders of magnitude to dopants/. This, however, comes at the price of severely lowering the structural quality of the samples. Such high concentrations of magnetic ions in semiconductors, such as [12, 13], leads to ferromagnetism with relatively high Curie temperatures of around —a temperature range that has opened the era of spintronics [14, 15].
The basic driving force of this form of ferromagnetism of impurities in III-V semiconductors is now understood  as the energy stabilization ensuing from the interaction between partially occupied hybrid orbitals located on different impurity sites, and , in the host crystal. Figure 1 (top) illustrates that if two such orbitals interact, there will be a gain in energy, i.e., a stabilization of the FM state (relative to the antiferromagnetic state ).
One might ask, how would such individual partially occupied hybrid orbitals come to exist in the first place? The answer to this is illustrated in Fig. 1 (bottom). The formation of a substitutional impurity can be thought of as a two step process, involving first the removal of a host atom, and second, the placement of a impurity atom in its place. The combined impurity/host orbitals simply result from coupling between the host-crystal cation vacancy state (“dangling bonds” of t(p) symmetry) and the impurity orbital of the atom filling this vacancy (having t(d) symmetry as well as nonbonding e(d) symmetry). The interaction between these states leads to a bonding state, and an antibonding state. The partial occupancy which ultimately drives ferromagnetism resides in the higher of these states. For in the antibonding state, dominated by t(p), is called a “dangling bond hybrid” (DBH) state, and the bonding/nonbonding states, dominated by t(d) and e(d), are called “crystal field resonances” (CFR). As we move to the left from in the periodic table, keeping the host, we approach the light impurity elements such as . For these, the atomic t(d) and e(d) are higher in energy and the roles of DBH and CFR are interchanged. The same interchange occurs in the case where is placed in a host crystal with light anion elements, such as nitrogen in , where the t(p) level is deeper and t(d) dominates the hole. It is the interaction between the partially occupied hybrid states, centered on different impurity sites in the crystal, which energetically stabilizes the ferromagnetic configuration.
This simple mechanism implies simple rules  for impurity-host material combinations that might produce ferromagnetism: one obtains the electronic configuration by constructing an Aufbau ground state for - electrons ( is the number of TM valence electrons, and is the valence of the cation site) occupying the relevant symmetry orbitals. For example, a ferromagnetism-promoting partial occupancy of DBH can be expected when occupies a trivalent site, but not the divalent site. This mechanism was directly verified by first-principles calculations of the energy of impurities in a solid as a function of the occupancy of the DBH, and noting the stabilization (destabilization) of the ferromagnetism for partial (full) occupancy.
Even higher promised in wider-gap semiconductors
First-principles electronic structure theory of doped into III-V nitrides, phosphides, and arsenides [18, 19, 20] have shown that the ferromagnetism-promoting hole (i.e., the orbital with partial occupancy, in Fig. 1) resides inside the host band gap. When the host anion changes from to and then to , the band gap increases primarily by shifting the valence band maximum (VBM) to lower energies. Thus, as one moves from arsenides to phosphides and then to nitride semiconductors, the hole-carrying impurity level, being roughly pinned, appears farther and farther from the VBM, deeper in the band gap. This makes the hole-carrying impurity wave function progressively more localized, weakening its communication with other magnetic ions; thus we lose ferromagnetism. Thus the conclusion of such microscopic calculations [4, 16, 17] was that will have lower as we go from a host arsenide to a host phospide and then to a nitride. While initial estimates of derived from directly associating the Weiss field with first-principles calculated magnetic coupling energies predicted optimistically higher in the wider gap (i.e., nitride) host crystals , it was quickly realized (by the same authors) that a proper statistical treatment of these exchange energies (via Monte Carlo, see Ref. ) predicted a much reduced as one moves to wider gap host crystals.
In parallel with such first-principles calculations, more traditional model-Hamiltonian approaches were developed [2, 23] that assume a scenario for the underlying interactions, and then solve the scenario mathematically. The classic approach in semiconductor physics, as described in Ref. , expands the states of a system in terms of a small number of preselected bands of the host crystal at the Brillouin-zone center. In the field of dilute magnetic semiconductors, the use of such concepts assumes that the hole that promotes magnetism can be described by just a few host crystal states, becoming essentially a “hostlike hole” [2, 23]). Such simple concepts differ substantially from the atomistic picture of Fig. 1 (bottom). The model Hamiltonian approach is then able to produce, analytically, the dependence of the ferromagnetism Curie temperature on basic material parameters. This lead to the highly influential bar-diagram published by Dietl et al.  (cited more than times) predicting the highest for binary semiconductors having the shortest bond length and thus the largest band gap. This design principle pointed towards wide-gap binary compounds containing a first-row element such as , , or as the best candidates for high- ferromagnets. This prediction started a world-wide quest for high- ferromagnets in wide-gap metal oxides (), nitrides (), and carbides based on dilute transition metal impurities. As exciting as such predictions were, they did not have the additional virtue of becoming true. Indeed, the assumed scenario that the ferromagnetism-promoting hole has the personality of the host valence band did not seem to be true for wide-gap insulators doped with transition metals such as . The simple lesson from microscopic models that the acceptor (hole-carrying) level becomes progressively more localized (i.e., contributing to lower ) as one moves up the group elements from arsenic to nitrogen was lost in the global enthusiasm to seek high in the widest gap host systems.
The role of theory as a guide: Pitfalls and fixes
Magnetic ions were not the only dilute species in wide-gap insulators that were implicated with high- ferromagnetism. Indeed, simple structural defects such as cation vacancies [25, 26, 27, 28] or carbon or nitrogen impurities [29, 30, 31] also became prime candidates for point sources of ferromagnetism behavior in wide-gap insulators. The basic thought here is that a proper combination of defect/impurity and a host crystal will lead to a hole-carrying (partially occupied) level, centered about this defect/impurity site, and that a ferromagnetic spin arrangement will be stabilized when different such centers interact [Fig. 1 (top)]. The search for the particular combination of defect/impurity and host crystal was, in many instances, guided by theory. Indeed, an impressive sign of the maturity of electronic structure theory based on first-principles was the central role that it had occupied in guiding the international quest for ferromagnetism in dilute wide-gap insulators [32, 33, 34, 35].
Unfortunately, a number of overlooked factors contributed to confusion. We next discuss a number of ways in which such calculations can produce unfounded optimism, such as a too high- ferromagnetism. Fortunately, recognition of such factors has lead to proposing a fix. Efforts worldwide to improve these fixes promise to create a robust guide for new design principles for such materials.
If it’s not one thing, it’s another
Guiding experiments towards promising defects or impurities that can promote carrier-induced ferromagnetism in certain wide-gap materials has proven to be difficult. Even if you correctly formulate the underlying electronic structure problem so that false-occupancy is avoided, the LDA failure to localize holes can get you. And if these are fixed, you still need to establish the defect concentration that satisfies percolation, yet is tolerated mechanically by the lattice. And you still need to consider the microstructure, i.e., the microscopic distribution of defects, since ideal single oxide crystals with ideal impurities do not lead to ferromagnetism. Indeed, if it’s not one thing that leads to false FM prediction, it’s another.
Finally, one should ask whether it is appropriate to describe the collective magnetism of a sample in terms of interacting single isolated point defects? Statistics tells us that the probability of two or more defects on one of the sublattices in a zinc blende, wurtzite, or rock salt crystal occupying neighboring sites is for defect concentrations ) . Once a defect pair has another defect next to it (e.g., with a probability), it has a lowered point group symmetry and consequently a different electronic spectrum, and is likely to exhibit completely different magnetic interactions compared to its isolated relatives [48, 70]. This has the detrimental effect that magnetic interactions no longer can be mapped on Dirac-van Vleck-Heisenberg type models ( . As expected, this problem can be circumvented by introducing additional correction terms or explicit configuration dependences to the model, as was done in Refs. [49, 51], which, however, mandates the explicit calculation of various additional clustered/microstructured configurations and their energetics. Moreover, the flaws in the underlying spin models seem to be of a fundamental nature, since, as it turns out, magnetic coupling energies, even of the simplest hydrogen molecules, have no mapping thereupon . This failure of Heisenberg’s theory  can further be traced back onto the Heitler-London theory , of which it is a generalization, and which has been abandoned on several occasions [75, 76, 77].
Bad is Good? No, this is not an attempt at moral revisionism. But consider this: The functionality of optoelectronic semiconductor devices such as lasers, light-emitting diodes, or photovoltaic solar cells is predicated on the structural quality of the active material. “Good materials,” i.e., those leading to efficient functionality, invariably need to have good crystalinity, purity, and very few defects. Yet if one were to use such near-perfect bulk samples of oxides, no ferromagnetism of the type discussed here would be present. Considered from this viewpoint, structurally imperfect and compositionally nonstoichiometric oxide materials (which would be considered “bad materials” in the semiconductor device culture) are “good materials” in this field of dilute magnetic oxides, as they are the only ones that yield, so far, the desired functionality.
The future of theory guide
Although much of the article has pointed out the various pitfalls encountered in running LDA on automatic pilot, the conclusions are more optimistic: Electronic structure calculations are a powerful tool to describe magnetic ground states, even in complicated clustered or nanostructured configurations. The ensuing results clarify which types of interactions are indeed critical, offering model-Hamiltonian approaches a safer, albeit ex-post-facto method to construct safe, simpler Hamiltonians. Furthermore, the results may be used as input for statistical models constructed to avoid the obvious pitfalls in earlier Heisenberg-like attempts.
This work was funded by the U.S. Department of Energy, Office of Science, Basic Energy Science, Materials Science and Energy Division, under Contract No. DE-AC36-08GO28308 to NREL. HR is supported by a Grant-in-Aid for Young Scientists (A) grant (No. 21686003) from the Japan Society for the Promotion of Science.
- M. A. Ruderman and C. Kittel, Phys. Rev. 96, 99 (1954); T. Kasuya, Prog. Theor. Phys. 16, 45 (1956); K. Yosida, Phys. Rev. 106, 893 (1957).
- T. Dietl et al., Science 287, 1019 (2000).
- V. F. Masterov, Fiz. Tekh. Poluprovodn. 18, 3 (1984).
- A. Zunger in Solid State Physics, edited by F. Seitz, D. Turnbull, and H. Ehrenreich (Academic Press, Orlando, 1986), Vol. 39, pp. 275-464[Amazon][WorldCat].
- F. Haldane and P.W. Anderson, Phys. Rev. B 13, 2553 (1976).
- A. Zunger and U. Lindefelt, Solid State Commun. 45, 343 (1983).
- H. Raebiger, S. Lany, and A. Zunger, Nature 453, 763 (2008).
- L. A. Ledebo and B. K. Ridley, J. Phys. 15, L961 (1982).
- M. J. Caldas, A. Fazzio, and A. Zunger, Appl. Phys. Lett. 45, 671 (1984).
- A. Zunger, Phys. Rev. Lett. 54, 849 (1985).
- H. Katayama-Yoshida and A. Zunger, Phys. Rev. Lett. 55, 1618 (1985).
- H. Munekata et al., Phys. Rev. Lett. 63, 1849 (1989).
- H. Ohno et al., Appl. Phys. Lett. 69, 363 (1996).
- A. Bonanni and T. Dietl, Chem. Soc. Rev. 39, 528 (2010).
- K. Kikoin, J. Magn. Magn. Mater. 321, 702 (2009).
- P. Mahadevan et al., Phys. Rev. Lett. 93, 177201 (2004); Phys. Rev. B 69, 115211 (2004).
- Y. J. Zhao et al., J. Appl. Phys. 98, 113901 (2005).
- K. Sato and H. Katayama-Yoshida, Jpn. J. Appl. Phys. 40, 485 (2001).
- L. Sandratskii and P. Bruno, Phys. Rev. B 67, 214402 (2003).
- T. Schulthess et al., Nature Mat. 4, 838 (2005).
- K. Sato et al., Europhys. Lett. 61, 403 (2003).
- K. Sato and H. Katayama-Yoshida, Phys. Rev. B 70, 201202 (2004).
- T. Jungwirth et al., Rev. Mod. Phys. 78, 809 (2006).
- P. Yu and M. Cardona, Fundamentals of Semiconductors: Physics and Materials Properties (Springer, New York, 2005)[Amazon][WorldCat].
- I. S. Elfimov, S. Yunoki, and G.A. Sawatzky, Phys. Rev. Lett. 89, 216403 (2002).
- M. Venkatesan, C. B. Fitzgerald, and J. M. D. Coey, Nature 430, 630 (2004).
- N. H. Hong et al., Phys. Rev. B 73, 132404 (2006).
- P. Dev, Y. Xue, and P. H. Zhang, Phys. Rev. Lett. 100, 117204 (2008).
- K. Kenmochi et al., Jpn. J. Appl. Phys. 43, L934 (2004).
- B. Gu et al., Phys. Rev. B 79, 024407 (2009).
- H. Peng et al., Phys. Rev. Lett. 102, 017201 (2009).
- K. Sato and H. Katayama-Yoshida, Jpn. J. Appl. Phys. 39, L555 (2000).
- O. Eriksson et al., J. Appl. Phys. 101, 09H114 (2007).
- H. Pan et al., Phys. Rev. B 77, 125211 (2008).
- K. Griffin Roberts et al., Phys. Rev. B 78, 014409 (2008).
- A. Filippetti et al., Chem. Phys. 309, 59 (2005).
- D. Iusan et al., Phys. Status Solidi A 204 53 (2007).
- S. Lany, H. Raebiger, and A. Zunger, Phys. Rev. B 77, 241201(R) (2008).
- N. E. Christensen, Phys. Rev. B 30, 5753 (1984).
- C. H. Patterson, Phys. Rev. B 74, 144432 (2006).
- C. D. Pemmaraju, T. Archer, D. Sánchez-Portal, and S. Sanvito, Phys. Rev. B 75, 045101 (2007).
- H. Kizaki et al., Appl. Phys. Express 2, 053004 (2009).
- H. Raebiger, S. Lany, and A. Zunger, Phys. Rev. Lett. 101, 027203 (2008).
- H. Raebiger, S. Lany, and A. Zunger, Phys. Rev. B 79, 165202 (2009).
- A. Walsh, J.L.F. Da Silva, and S.H. Wei, Phys. Rev. Lett. 100, 256401 (2008).
- S. Sanvito and C.Das Pemmaraju, Phys. Rev. Lett. 102, 159701 (2009).
- J. Osorio-Guillén, S. Lany, S. V. Barabash, and A. Zunger, Phys. Rev. Lett. 96, 107203 (2006).
- H. Raebiger et al., Phys. Rev. B 72, 014465 (2005).
- T. Hynninen et al., J Phys. Condens. Matter 18, 1561 (2006).
- A. I. Liechtenstein et al., J. Magn. Magn. Mater. 67, 65 (1987).
- A. Franceschetti et al., Phys. Rev. Lett. 97, 047202 (2006).
- J. Osorio-Guillén, S. Lany, S. V. Barabash, and A. Zunger, Phys. Rev. B 75, 184421 (2007).
- C. G. van de Walle and J. Neugebauer, J. Appl. Phys. 95, 3851 (2004).
- C. Persson, Y. J. Zhao, S. Lany, and A. Zunger, Phys. Rev. B 72, 035211 (2005).
- S. Lany and A. Zunger, Phys. Rev. B 78, 235104 (2008).
- A. M. Stoneham, J. Phys. Condens. Matter 22, 074211 (2010).
- J. P. Perdew et al., Phys. Rev. Lett. 49, 1691 (1982).
- J. P. Perdew et al., Phys. Rev. A 76, 040501 (2007).
- P. Mori-Sánchez, A. J. Cohen, and W. Yang, Phys. Rev. Lett. 100, 146401 (2008).
- S. Lany and A. Zunger, Phys. Rev. B 80, 085202 (2009).
- A. L. Shluger et al., Modelling Simul. Mater. Sci. Eng. 17, 084004 (2009).
- M. d’Avezac, M. Calandra, and F. Mauri, Phys. Rev. B 71, 205210 (2005).
- A. Droghetti, C. D. Pemmaraju, and S. Sanvito, Phys. Rev. B 78, 140404 (2008).
- S. Lany and A. Zunger (unpublished).
- S. Lany and A. Zunger, Phys. Rev. B 81, 205209 (2010).
- S. Lany and A. Zunger, Appl. Phys. Lett. 96, 142114 (2010).
- Q. Wang et al., Phys. Rev. B 77, 205411 (2008).
- J. A. Chan, S. Lany, and A. Zunger, Phys. Rev. Lett. 103, 016404 (2009).
- R. Behringer, J. Chem. Phys. 29, 537 (1958).
- H. Raebiger et al., J Phys. Condens. Matter 16, L457 (2004).
- M. van Schilfgaarde and O. Mryasov, Phys. Rev. B 63, 233205 (2001).
- I. Ciofini et al., Chem. Phys. 309, 133 (2005).
- W. Heisenberg, Z. Phys. 49, 619 (1928).
- W. Heitler and F. London, Z. Phys. 44, 455 (1927).
- L. Pauling, The Nature of the Chemical Bond (Cornell University Press, Ithaca, 1960)[Amazon][WorldCat].
- J. C. Slater, Rev. Mod. Phys. 25, 199 (1953).
- K. Hongo et al., Mat. Trans. 48, 662 (2007).