This is an open access article published under an ACS AuthorChoice License, which permits copying and redistribution of the article or any adaptations for non-commercial purposes.
Article Cite This: ACS Omega 2018, 3, 5301−5307
Theory of Ferromagnetism in Reduced ZrO2−x Nanoparticles Elisa Albanese, Antonio Ruiz Puigdollers, and Gianfranco Pacchioni* Dipartimento di Scienza dei Materiali, Università di Milano-Bicocca, via R. Cozzi 55, 20125 Milano, Italy S Supporting Information *
ABSTRACT: Bulk ZrO2 is both nonreducible and nonmagnetic. Recent experimental results show that dopant-free, oxygendeficient ZrO2−x nanostructures exhibit a ferromagnetic behavior at room temperature (RT). Here, we provide a comprehensive theoretical foundation for the observed RT ferromagnetism of zirconia nanostructures. ZrO2 nanoparticles containing up to 700 atoms (3 nm) have been studied with the help of density functional theory. Oxygen vacancies in ZrO2 nanoparticles form more easily than in bulk zirconia and result in electrons trapped in 4d levels of low-coordinated Zr ions. Provided the number of these sites exceeds that of excess electrons, the resulting ground state is high spin and the ordering is ferromagnetic. The work provides a general basis to explain magnetism in intrinsically nonmagnetic oxides without the help of dopants.
1. INTRODUCTION Ferromagnetism in semiconductor metal oxides finds important applications in spintronics and optoelectronics.1,2 Many efforts have been dedicated to the study of ferromagnetism (FM) induced by transition metal (TM) doping of semiconductor oxides such as ZnO,3,4 TiO2,5−75−7 and ZrO2,8,98,9 also called diluted magnetic semiconductors (DMSs). However, it remains unclear whether FM is an intrinsic property of the system or is due to the aggregation of magnetic metal dopants or other magnetic impurities.10,11 Oxides doped by light 2p elements (e.g., C- and N-ZnO,12,13 C- and N-TiO2,14,15 N-ZrO2,16 and N-BaTiO317) have been also considered and predicted as ferromagnets at room temperature (RT), but the origin of the magnetic behavior is a matter of debate.18−20 In other cases, staying with the material of interest for this study, tetragonal zirconia, no ferromagnetism has been observed by doping the oxide with Mn.21 Encouraging results have been recently obtained for undoped thin films of TiO2, ZnO, In2O3, and HfO2,22−26 which are FM at RT without any type of dopant element. Also, thin films of ZrO2 exhibit RT FM. It has been suggested that this is related to the presence of intrinsic defects. The crystallographic phase also appears to play a role, with tetragonal ZrO2 that seems to be essential.27−29 In thin films, however, the necessity to precisely control the nature of the interface with the support is a critical issue. For this reason, metal oxide nanoparticles (NPs) and nanostructures represent ideal systems to obtain FM at high Curie temperatures (TC > RT).30 Recently, a FM behavior at RT has been reported for undoped ZrO2 nanostructures.2 The phenomenon has been attributed to the presence of surface oxygen vacancies (VO). However, there is a conceptual problem with this hypothesis. A neutral oxygen vacancy in the bulk and © 2018 American Chemical Society
on the surface of zirconia consists of two electrons trapped in a cavity, with a diamagnetic ground state.31,32 This means that removing oxygen from the bulk or the surface of zirconia does not necessarily lead to the appearance of magnetic states. Another problem is the low reducibility of the material. The formation of VO centers in the bulk or on the surface of ZrO2 is thermodynamically very unfavorable, with vacancy formation energies of about 6 eV (computed with respect to 1/2 O2).33 In fact, ZrO2 belongs to the class of nonreducible oxides. However, density functional theory (DFT) calculations have shown that the cost to remove oxygen in ZrO2 NPs is substantially lower and that the excess electrons associated with a VO defect can lead in some cases to the formation of Zr3+ at the position of low-coordinated (LC) Zr ions.33 This suggests a special role of nanostructuring in turning the nonreducible ZrO2 into a reducible material. But under which conditions the Zr3+ centers can lead to a FM ordering? Is this ordering stable at RT? What size of NP is required to reach these conditions? What is the level of reduction requested? To answer these questions, we have performed an extensive theoretical investigation of nonstoichiometric ZrO2−x NPs characterized by increasing the size and level of oxygen deficiency. For the calculations, we have used both DFT + U and hybrid (PBE0) functionals. Hybrid functionals, in particular, are needed to provide a correct description of the electronic structure34 and to accurately reproduce the strength of the magnetic coupling and of the resulting Curie temperature (TC) of the finite systems.35−37 Received: April 6, 2018 Accepted: May 4, 2018 Published: May 17, 2018 5301
DOI: 10.1021/acsomega.8b00667 ACS Omega 2018, 3, 5301−5307
Article
ACS Omega
2.1. Nonmagnetic versus Magnetic Ground State. The systems considered in this study are finite in size, and as such, their electronic structure should be discussed in terms of discrete levels and molecular orbitals. However, at the sizes considered, bands start to form; furthermore, the study aims at comparing NPs with the extended solid (bulk). For these reasons, we will make use of the language of the band structure to address the electronic structure of the zirconia NPs. The relative stability of diamagnetic and high-spin (HS) solutions of ZrO2 NPs has been evaluated by means of PBE + U exploratory calculations. For the smallest NP considered, Zr44O80, the ground state is a singlet closed shell. The lowest HS solution, characterized by 16 unpaired electrons, lies 1.69 eV higher in energy. In the ground state (singlet), the excess electrons are localized on the dz2 orbital of the six Zr4c corner ions, which are reduced from Zr4+(4d0) to Zr2+(4d2) (see Figure 2). This is because the 4d states of the low-coordinated Zr ions (ZrLC) are lower in energy than those of the fully coordinated Zr ones and appear as mid-gap states. In Zr44O80, they lie about 1 eV below the conduction band (CB) (see Figure S2).38 The remaining excess of charge is distributed inside the NP. In the HS model, the excess electrons are localized on the ZrLC ions at corners and edges, which are reduced from Zr4+(4d0) to Zr3+(4d1). The electron localization on the Zr ions leads to a remarkable deformation of the Zr first coordination sphere, corresponding to the formation of a small polaron. The small number of Zr6c ions along each edge (only two in Zr44O80) implies a dense distribution of the unpaired electrons, with a higher cost for the polaronic distortion. This is the origin of the low stability of the HS state. On increasing the NP size, the nature of the ground state changes. In the medium-size NP considered, Zr85O160, there are 20 excess electrons (10 missing O atoms) and the ground state is HS and FM, with the singlet closed-shell state lying 0.86 eV higher in energy (see Table 1). The order of diamagnetic and magnetic solutions is thus reversed with respect to the smaller NP. The distribution of the 20 unpaired electrons is similar to that in Zr44O80, but now the ZrLC3+ ions are more uniformly and homogeneously distributed over the entire surface due to the higher number of Zr6c sites along the edges (see Table 1). No evidence of clustering of Zr3+ sites is found. This is an important conclusion. The switch from the diamagnetic to the magnetic ground state depends on the ratio between the number of excess electrons and the number of ZrLC sites available (see Table 1). If ZrLC ≫ n(e−), the magnetic ground state is preferred. In the absence of a sufficient number of lowlying acceptor levels, in fact, electrons are forced to doubly occupy these states, leading to a nonmagnetic ground state. This interpretation is further corroborated by the results for the larger Zr146O280 NPs. Here, only the FM solution is obtained; various attempts to converge on a diamagnetic state failed because of the intrinsic instability of the singlet closedshell state compared to that of the HS one (see Table 1). To summarize, a larger size of the ZrO2−x NPs implies that an increase of the number of ZrLC ions is able to stably trap excess electrons, favoring a single occupancy of the Zr 4d levels (Zr3+). The results clearly show that an O-deficient ZrO2−x NP has a magnetic ground state, at variance with the extended surface or the bulk of reduced zirconia. The spontaneous magnetization is thus an intrinsic property of nanostructured, O-deficient zirconia. 2.2. Magnetic Ordering and the Origin of Magnetic Behavior in Zirconia NPs. To study the existence of RT FM
We show that the high-spin (HS) FM solution is, beyond a given size of the NP, the ground state of the system. This work provides a solid theoretical basis for the observed magnetic behavior of nonstoichiometric ZrO2 NPs and nanostructures.
2. RESULTS AND DISCUSSION We model ZrO2−x NPs following the relative stabilities of the orientation of the lattice planes in bulk tetragonal zirconia (tZrO2). This is analogous to a Wulff construction at the nanoscale.33,38 In particular, bulk t-ZrO2 was cleaved along the O-terminated (101) surface, the most stable crystallographic face (see Supporting Information). This structure has also been observed experimentally for CeO2 NPs, which have similar characteristics.39 For a discussion on the stability of oxide nanostructures, see also ref 40. We considered tetragonal ZrO2 because FM has been observed for this specific polymorph.27−29 We built four nonstoichiometric (oxygen-deficient) octahedral nanostructures: Zr44O80, Zr85O160, Zr146O280, and Zr231O448 (see Figure 1). In the smallest NP, Zr44O80, eight O atoms are
Figure 1. Optimized structures of ZrO2 nanoparticles. Red spheres represent the O atoms, and light blue spheres, Zr atoms.
missing with respect to the stoichiometric composition. In the larger particles, the numbers of missing O atoms are 10, 12, and 14, respectively. The size of the NPs goes from 1.4 nm to almost 3.0 nm, which is the same range as that of synthesized octahedral-based zirconia NPs.41 The NPs are characterized by the presence of six 4-coordinated Zr corners (Zr4c) and by a size-dependent number of 6-coordinated Zr edges (Zr6c): from two in Zr44O80 up to five in Zr231O448. 5302
DOI: 10.1021/acsomega.8b00667 ACS Omega 2018, 3, 5301−5307
Article
ACS Omega
Figure 2. (a) Spin densities of three ZrO2−x NPs in the high-spin solution; (b) charge density plots of the excess electrons of the singlet NPs. ρ = 0.01 e−/Å3.
Table 1. Number of Excess Electrons (e−) and of Zr Edge and Zr Corner Atoms and Relative Stabilities (eV) of the Singlet Closed Shell (S) vs Those of High-Spin (HS) Solutions of ZrO2−x NPsa Zr44O80 Zr85O160 Zr146O280 a
e−
Zr4c
Zr6c
S
HS
16 20 24
6 6 6
24 36 48
0.0 +0.86 not conv
+1.69 0.0 0.0
PBE + U results.
in ZrO2−x NPs, we have compared the stabilities of FM and antiferromagnetic (AFM) solutions using a broken-symmetry solution (the calculations described in what follows are based on the PBE0 hybrid functional and CRYSTAL14 code; all structures have been geometrically optimized). The first-principles AFM solution has been computed for Zr146O280 starting from the optimal FM geometry and then fully relaxing the structures (see Figure 3). However, the two structures are virtually identical. In Zr146O280, the cost to flip 12 spins, going from 24 up ↑ electrons (FM) to 12 up ↑ and 12 down ↓ electrons (AFM), is 13.3 meV, confirming that the FM ordering is the ground state (see Table 2). We considered an even larger zirconia NP, Zr231O448 (see Figure 3). Notice that this particle has nearly 700 atoms and a total of 28 unpaired electrons. For this large NP, the FM state is more stable than the AFM one by 30.4 meV (see Table 2). The trend is then confirmed: a higher number of missing O atoms in the NPs leads to a more stable FM state. The FM coupling can be explained by bound magnetic polaron (BMP) theory.11,42 The model works for transition metal (TM) dopants and implies an exchange interaction
Figure 3. Spin density plots of (a, c) high-spin ferromagnetic (FM) and (b, d) antiferromagnetic (AFM) solutions of Zr146O280 and Zr231O448 NPs, respectively. Yellow contours represents spin up ↑ and blue contour represents spin down ↓ electrons. ρ = 0.01 e−/Å3.
5303
DOI: 10.1021/acsomega.8b00667 ACS Omega 2018, 3, 5301−5307
Article
ACS Omega
represented by the Zr3+ ions. The stability of the FM state depends on the hybridization degree between the impurity band and the CB. Because the nature of the CB in turn depends on the polymorph and the crystallographic facets involved, this model can explain the effect of the structure on the FM alignment. As reported in Table 2, the FM ordering in Zr146O280 is more stable than the AFM one by about 13 meV. This means that the AFM state is not expected to survive at RT because the value of 13.3 meV is well below the thermal energy (26 meV). In contrast, the higher magnetic coupling obtained for Zr231O448 (30 meV) suggests that the ferromagnetism could exist at higher temperatures and could be stable against thermal fluctuations.45 In fact, increasing the number of BMPs, the total magnetization and then the TC are also expected to increase. 2.3. Magnetic Coupling in Zirconia NPs. Finally, we have computed the exchange coupling (Jij) of a pair of isolated Zr ions. This requires us to compare the stabilities of FM and antiferromagnetic (AFM) solutions using a broken-symmetry solution (PBE0 results, CRYSTAL14 code). An approximated approach has been adopted. We considered the stoichiometric Zr40O80 NP, where a LC oxygen has been removed (Zr40O80VO) from one of the most favorable sites (low O vacancy formation energy) to study with first principles the exchange interaction between just two magnetic moments. This leaves two electrons that can give a FM (triplet), an AFM (singlet open shell), or a diamagnetic (singlet closed shell, S) solution. In the open-shell solutions, one of the two excess electrons is localized on a Zr4c corner ion and the other in the vacancy (Figure 5). The singlet closed-shell configuration with electron
Table 2. Energy Difference between AFM and FM Solutions, ΔE (AFM − FM), for Zr146O280 and Zr231O448 NPsa
a
NP
ΔE (AFM − FM), meV
Zr146O280 Zr231O448
+13.3 +30.4
PBE0 results.
between shallow electrons associated with defects and localized d electrons of the TM dopants. This approach has been reconsidered by taking into account only the oxygen vacancy (VO) as the defect responsible for the observed FM in undoped-TiO2 nanoribbons,43 and the polaron considered is associated with oxygen vacancies. Once the VO concentration reaches a certain limit, an overlap among the BMPs is established, thus enhancing the FM behavior. However, it should be mentioned that the model cannot explain the absence of RT FM in some oxide films, such as monoclinic HfO2 and ZrO2, characterized by a large number of oxygen vacancies.27 Recently, the charge transfer ferromagnetism (CTF) model has been proposed.44 According to this model, the presence of defects introduces an impurity band below the CB. The CTF model involves the presence of a charge reservoir that facilitates the hopping of electrons between the impurity band and the CB, leading to the splitting of spin states (see the inset of Figure 4).
Figure 4. Projected density of states of Zr146O280 NP. The contributions of O (red), Zr4+ 4d0 (blue), Zr3+ 4d1 corner (black), and Zr3+ 4d1 edge (purple) atoms are reported. The Fermi energy is set to the highest occupied level. The schematically proposed mechanism of the FM interaction is due to the hybridization of the Zr3+ states with the CB (composed by Zr4+ 4d and 5s orbitals), which leads to the splitting of up and down states, as shown in the inset. A simplified Zr146O280 model with excess electrons in the FM configuration is also shown.
Figure 5. Spin density plot of high-spin (a) ferromagnetic and (b) antiferromagnetic solutions for a Zr40O80 NP with VO at the surface. Yellow represents up ↑ and blue represents down ↓ electrons. ρ = 0.008 e/Å3.
localization in the vacancy is favored by the presence of the void formed in the position of the missing oxygen atom (as for zirconia surface and bulk). This further shows that the magnetic behavior appears only in the presence of several vacancies and the related excess of charge. Then, the exchange coupling (Jij) of two electrons localized on two low-coordinated Zr atoms and considered as an isolated pair (ij), Jij = ΔE (AFM − FM), is computed. The FM solution is 1.0 meV more stable than the AFM one (J = 1.0 meV). Of course, this value may slightly change if the vacancy is formed in another position or if it is created in a particle of different
Here, we propose a model that is a hybrid of these two theories to explain the FM behavior of undoped ZrO2 and possibly other oxides.2,43 The model implies the formation of a BMP among electrons localized on two nearby Zr centers. As shown in Figure 4, at high vacancy concentrations, these defects lead to the formation of a band in the gap just below the CB. The defect states, composed by the Zr 4d1 states, can merge with the CB, leading to a splitting of up and down spin states that favors the FM ordering. The charge reservoir is 5304
DOI: 10.1021/acsomega.8b00667 ACS Omega 2018, 3, 5301−5307
Article
ACS Omega size. However, the final goal of this model is to compute a realistic J value for two magnetic Zr sites with similar environment and separated by more or less the same average distance observed among magnetic Zr ions in Zr146O280 and Zr231O448 nanoparticles. From this point of view, we believe that the J value computed is representative of the magnetic interaction on larger zirconia reduced nanoparticles. Various FM and AFM solutions with different geometries have been considered (see Figure S3 and Table S3). It turns out that the FM state is always lower in energy. The preference for a FM coupling seems to be an intrinsic property of the system. To estimate the Curie temperature, one can use the meanfield approximation (MFA) and the Heisenberg model, according to the following equation46 kBTC =
2 J 3c
4. COMPUTATIONAL DETAILS For the study of ZrO2−x NPs, we used two different approaches. Both are based on the use of codes that use periodic boundary conditions. A NP is contained in a supercell and separated by sufficiently large distances from their replicas in adjacent supercells. The structures of ZrO2−x have been geometrically optimized by means of spin-polarized DFT calculations with the Vienna Ab-initio Simulation Package (VASP 5.3)51−53 (plane wave basis set with a kinetic energy cutoff of 400 eV; core electrons described by the projector-augmented wave method).54,55 The generalized gradient approximation (GGA) for the exchange functional was used within the Perdew, Burke, and Ernzerhof (PBE) formulation.56 The self-interaction error of GGA functionals was partly corrected using the PBE + U approach57,58 with the onsite Coulomb correction (Ueff = U − J) for the Zr(4d) states set to 4 eV.31 Geometry optimizations of the NP were carried out at the Γ-point, with all atoms free to relax until the ionic forces become smaller than |0.01| eV/Å (self-consistency of the electron density set to 10−5 eV). The magnetic nature of the finite nanoparticles has been evaluated with the hybrid DFT functional PBE059,60 using the CRYSTAL14 code.61 The structures have been fully optimized at the PBE0 level. Crystalline orbitals are represented as linear combinations of Bloch functions (BFs). Each BF is built from linear combinations of Gaussian-type orbitals. An all-electron basis set has been used for O, 8-411(d1), whereas for Zr, the Hay and Wadt small-core effective core potential62 associated with a 311(d2) basis set was selected. For the numerical integration of the exchange-correlation term, 75 radial points and 974 angular points (XLGRID) in a Lebedev scheme in the region of chemical interest were adopted. The accuracy of the integral calculations was increased with respect to the default value by setting the tolerances to 10, 10, 10, 10, and 20. The self-consistent field iterative procedure converged to a tolerance in total energy of ΔE = 1 × 10−8 au. The threshold for the maximum and the root-mean-square forces were set to 0.00045 and 0.0003 au, respectively.
(1)
The MFA method does not deal with the statistical problem of the local magnetic moments and tends to overestimate TC.47−49 Nevertheless, it has given acceptable results.49 The J value is that obtained for a pair of Zr ions as discussed above, 1.0 meV, and c is the Zr3+ concentration. Although this is well defined in the bulk of a magnetic solid, it is much more complex to define on a nanoparticle where the magnetic moments localize on lowcoordinated ions at the surface of the nanoparticle. The value of TC depends critically on this quantity. If we define c as the ratio between the number of Zr3+ ions in the NP and the total number of Zr ions, for Zr231O448 (c = 28/231 = 0.12), we estimate a Curie temperature TC = 64 K. If we assume a hypothetical concentration of 2% of Zr3+ ions (c = 0.02), we obtain TC = 387 K. Clearly, it is impossible to provide a quantitative estimate of TC without knowing the exact concentration of Zr3+ ions.
3. CONCLUSIONS Bulk ZrO2 is both nonreducible and nonmagnetic. Recent experimental results show that dopant-free, oxygen-deficient ZrO2−x nanostructures exhibit FM behavior at RT. We have recently shown that ZrO2 NPs are easy to reduce via O2 or H2O desorption from low-coordinated sites.50 In this work, we provide a theoretical foundation for the observed magnetic character of zirconia nanostructures. In ZrO2−x NPs, the presence of low-coordinated Zr ions favors the localization of the excess electrons originated from the oxygen deficiency in the low-lying Zr 4d levels (Zr3+). The magnetic Zr3+(4d1) defects can form only if a sufficient number of ZrLC ions exist. When this is not the case, the excess of charge is accommodated in doubly occupied 4d levels, Zr2+(4d2), leading to an overall diamagnetic ground state. Thus, there is a clear relationship among the maximum possible level of reduction (number of excess electrons), the surface area of the NP, and the nature of the ground state (magnetic or diamagnetic). If ZrLC sites exceed the number of excess electrons, the ground state is high spin and the preferred ordering is ferromagnetic. Increasing the number of magnetic defects also leads to a stronger exchange coupling and higher Curie temperatures. This work provides a solid theoretical ground to explain the completely different magnetic behavior of oxygen-deficient nanostructured zirconia compared to that of the bulk material.
■
ASSOCIATED CONTENT
S Supporting Information *
The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acsomega.8b00667. Details about the construction of zirconia NPs, the density of states of all of the nonstoichiometric NPs, and various solutions for the exchange coupling in the Zr40O80-VO nanoparticle (PDF)
■
AUTHOR INFORMATION
Corresponding Author
*E-mail:
[email protected]. ORCID
Elisa Albanese: 0000-0002-8255-3658 Gianfranco Pacchioni: 0000-0002-4749-0751 Author Contributions
G.P. proposed the theoretical studies. E.A. and A.R.P. performed the calculations and the analysis of the results. All authors discussed and interpreted the results, developed the theoretical models, and contributed to the preparation and writing of the manuscript. Notes
The authors declare no competing financial interest. 5305
DOI: 10.1021/acsomega.8b00667 ACS Omega 2018, 3, 5301−5307
Article
ACS Omega
■
(19) Albanese, E.; Pacchioni, G. Ferromagnetism in Nitrogen-doped BaO: A Self-interaction Corrected DFT Study. Phys. Chem. Chem. Phys. 2017, 19, 3279−3286. (20) Albanese, E.; Leccese, M.; Di Valentin, C.; Pacchioni, G. Magnetic Properties of Nitrogen-doped ZrO2: Theoretical Evidence of Absence of Room Temperature Ferromagnetism. Sci. Rep. 2016, 6, No. 31435. (21) Srivastava, S. K.; Lejay, P.; Barbara, B.; Boisron, O.; Pailhès, S.; Bouzerar. Absence of ferromagnetism in Mn-doped tetragonal zirconia. J. Appl. Phys. 2011, 110, No. 043929. (22) Venkatesan, M.; Fitzgerald, C. B.; Coey, J. M. Thin Films: Unexpected Magnetism in a Dielectric Oxide. Nature 2004, 430, 630. (23) Hong, N. H.; Sakai, J.; Poirot, N.; Brizé, V. Room-Temperature Ferromagnetism Observed in Undoped Semiconducting and Insulating Oxide Thin Films. Phys. Rev. B 2006, 73, No. 132404. (24) Ning, S.; Zhan, P.; Xie, Q.; Li, Z.; Zhang, Z. Room-Temperature Ferromagnetism in Undoped ZrO2 Thin Films. J. Phys. D: Appl. Phys. 2013, 46, No. 445004. (25) Sun, S.; Wu, P.; Xing, X. d0 Ferromagnetism in Undoped n and p-type In2O3 Films. Appl. Phys. Lett. 2012, 101, No. 132417. (26) Gupta, A.; Cao, H.; Parekh, K.; Rao, K. V.; Raju, A. R.; Waghmare, U. V. Room Temperature Ferromagnetism in Transition Metal (V, Cr, Ti) Doped In2O3. J. Appl. Phys. 2007, 101, No. 09N513. (27) Ning, S.; Zhan, P.; Xie, Q.; Wang, W.; Zhang, Z. Defect Driven Ferromagnetism in Undoped Dilute Magnetic Oxides: A Review. J. Mater. Sci. Technol. 2015, 31, 969−978. (28) Kumar, S.; Bhunia, S.; Ojha, A. K. Experimental and Theoretical Evidence for the Presence of Room Temperature Ferromagnetism in Undoped and Mn Doped Tetragonal ZrO2 Nanostructures. Chem. Phys Lett. 2016, 644, 271−275. (29) Ning, S.; Zhan, P. Phase-Dependent and Defect-Driven d0 Ferromagnetism in Undoped ZrO2 Thin Films. RSC Adv. 2015, 5, 3636−3641. (30) Long, L.; Xiong, S.; Meng, M.; Liu, L.; Zhang, J.; Wu, X.; Chu, P. K. Enhancement of Ferromagnetism in Nonmagnetic Metal Oxide Nanoparticles by Facet Engineering. Small 2017, 13, No. 1602951. (31) Syzgantseva, O. A.; Calatayud, M.; Minot, C. Revealing the Surface Reactivity of Zirconia by Periodic DFT Calculations. J. Phys. Chem. C 2012, 116, 6636−6644. (32) Hur, J.-H.; Park, S.; Chung, U.-I. First Principles Study of Oxygen Vacancy States in Monoclinic ZrO2: Interpretation of Conduction Characteristics. J. Appl. Phys. 2012, 112, No. 113719. (33) Puigdollers, A. R.; Illas, F.; Pacchioni, G. Structure and Properties of Zirconia Nanoparticles from Density Functional Theory Calculations. J. Phys. Chem. C 2016, 120, 4392−4402. (34) Gerosa, M.; Bottani, C. E.; Di Valentin, C.; Onida, G.; Pacchioni, G. Accuracy of Dielectric-Dependent Hybrid Functionals in the Prediction of Optoelectronic Properties of Metal Oxide Semiconductors: a Comprehensive Comparison with Many-Body GW and Experiments. J. Phys.: Condens. Matter 2018, 30, No. 044003. (35) Dinh, V. A.; Toyoda, M.; Sato, K.; Katayama-Yoshida, H. Exchange Interaction and Tc in Alkaline-Earth-Metal-Oxide-Based DMS without Magnetic Impurities: First principle Pseudo-SIC and Monte Carlo Calculation. J. Phys. Soc. Jpn. 2006, 75, No. 093705. (36) Illas, F.; Moreira, I. P. R.; de Graaf, C.; Barone, V. Magnetic Coupling in Biradicals, Binuclear Complexes and Wide-Gap Insulators: a Survey of Ab Initio Wave Function and Density Functional Theory Approaches. Theor. Chem. Acc. 2000, 104, 265−272. (37) Moreira, I. d. P. R.; Illas, F. Unified View of the Theoretical Description of Magnetic Coupling in Molecular Chemistry and Solid State Physics. Phys. Chem. Chem. Phys. 2006, 8, 1645−1659. (38) Ruiz-Puigdollers, A.; Illas, F.; Pacchioni, G. ZrO2 Nanoparticles: A Density Functional Theory Study of Structure, Properties and Reactivity. Rend. Lincei 2017, 28, 19−27. (39) Zhang, F.; Jin, Q.; Chon, S. W. Ceria Nanoparticles: Size, Size Distribution, and Shape. J. Appl. Phys. 2004, 95, 4319−4326. (40) Bromley, S. T.; Moreira, I.; Neyman, K. M.; Illas, F. Approaching Nanoscale Oxides: Models and Theoretical Methods. Chem. Soc. Rev. 2009, 38, 2657−2670.
ACKNOWLEDGMENTS We thank Prof. F. Illas for useful discussions. This work has been supported by Italian MIUR through the PRIN Project “2015K7FZLH SMARTENESS” and by the FP7/2007-2013 Marie Curie Network CATSENSE under Grant Agreement number 607417. CINECA-LISA Awards are also acknowledged for the availability of high-performance computing resources.
■
REFERENCES
(1) Fert, A. Nobel Lecture: Origin, Development, and Future of Spintronic. Rev. Mod. Phys. 2008, 80, 1517−1529. (2) Rahman, M. A.; Rout, S.; Thomas, J. P.; McGillivray, D.; Leung, K. T. Defect-Rich Dopant-free ZrO2 Nanostructures with Superior Dilute Ferromagnetic Semiconductor Properties. J. Am. Chem. Soc. 2016, 138, 11896−11906. (3) Pereira, L. M. C.; Araújo, J. P.; Wahl, U.; Decoster, S.; Van Bael, M. J.; Temst, K.; Vantomme, A. Searching for Room Temperature Ferromagnetism in Transition Metal Implanted ZnO and GaN. J. Appl. Phys. 2013, 113, No. 023903. (4) Straumal, B. B.; Mazilkin, A. A.; Protasova, S. G.; Myatiev, A. A.; Straumal, P. B.; Schütz, G.; van Aken, P. A.; Goering, E.; Baretzky, B. Magnetization Study of Nanograined Pure and Mn-doped ZnO Films: Formation of a Ferromagnetic Grain-Boundary Foam. Phys. Rev. B. 2009, 79, No. 205206. (5) Matsumoto, Y.; et al. Room-Temperature Ferromagnetism in Transparent Transition Metal-doped Titanium Dioxide. Science 2001, 291, 854−856. (6) Zhang, S.; et al. Electronic Manifestation of Cation-VacancyInduced Magnetic Moments in a Transparent Oxide Semiconductor: Anatase Nb:TiO2. Adv. Mater. 2009, 21, 2282−2287. (7) Wang, S.; Pan, L.; Song, J. J.; Mi, W.; Zou, J. J.; Wang, L.; Zhang, X. Titanium-Defected Undoped Anatase TiO2 with p-type Conductivity, Room-Temperature Ferromagnetism, and Remarkable Photocatalytic Performance. J. Am. Chem. Soc. 2015, 137, 2975−2983. (8) Ostanin, S.; et al. Mn-Stabilized Zirconia: from Imitation Diamonds to a New Potential High-Tc Ferromagnetic Spintronics Material. Phys. Rev. Lett. 2007, 98, No. 016101. (9) Zippel, J.; et al. Defect-Induced Ferromagnetism in Undoped and Mn-Doped Zirconia Thin Films. Phys. Rev. B 2010, 82, No. 125209. (10) Dietl, T. A Ten-year Perspective on Dilute Magnetic Semiconductors and Oxides. Nat. Mater. 2010, 9, 965−974. (11) Ogale, S. B. Dilute Doping, Defect, and Ferromagnetism in Metal Oxide Systems. Adv. Mater. 2010, 22, 3125−3155. (12) Pan, H.; Yi, J. B.; Shen, L.; Wu, R. Q.; Yang, J. H.; Lin, J. Y.; Feng, Y. P.; Ding, J.; Van, L. H.; Yin, J. H. Room-Temperature Ferromagnetism in Carbon-Doped ZnO. Phys. Rev. Lett. 2007, 99, No. 127201. (13) Shen, L.; Wu, R. Q.; Pan, H.; Peng, G. W.; Yang, M.; Sha, Z. D.; Feng, Y. P. Mechanism of Ferromagnetism in Nitrogen-Doped ZnO: First-principle Calculations. Phys. Rev. B: Condens. Matter Mater. Phys. 2008, 78, No. 073306. (14) Tao, J. G.; Guan, L. X.; Pan, J. S.; Huan, C. H. A.; Wang, L.; Kuo, J. L.; Zhang, Z.; Chai, J. W.; Wang, S. Density Functional Study on Ferromagnetism in Nitrogen-Doped Anatase TiO2. Appl. Phys. Lett. 2009, 95, No. 062505. (15) Bao, N. N.; Fan, H. M.; Ding, J.; Yi, J. B. Room Temperature Ferromagnetism in N-doped Rutile TiO2 Films. J. Appl. Phys. 2011, 109, No. 07C302. (16) Zhu, H.; Li, J.; Chen, K.; Yi, X.; Cheng, S.; Gan, F. Nature of Charge Transport and p- electron Ferromagnetism in Nitrogen-doped ZrO2: An Ab-initio Perspective. Sci. Rep. 2015, 5, No. 8586. (17) Tan, X.; Chen, C.; Jin, K.; Luo, B. Room-Temperature Ferromagnetism in Nitrogen-doped BaTiO3. J. Alloys Compd. 2011, 509, L311−L313. (18) Seike, M.; Dinh, V. A.; Sato, K.; Katayama-Yoshida, H. Firstprinciples Study of the Magnetic Properties of Nitrogen-doped Alkaline Earth Metal Oxides. Phys. B 2012, 407, 2875−2878. 5306
DOI: 10.1021/acsomega.8b00667 ACS Omega 2018, 3, 5301−5307
Article
ACS Omega (41) Grena, R.; Masson, O.; Portal, L.; Reḿ ondier̀e, F.; Berghout, A.; Jouin, J.; Thomas, P. Stabilization Effect of Surface Impurities on the Structure of Ultrasmall ZrO2 Nanoparticles: An ab-initio Study. J. Phys. Chem. C 2015, 119, 15618−15626. (42) Kaminski, A.; Das Sarma, S. Polaron Percolation in Diluted Magnetic Semiconductors. Phys. Rev. Lett. 2002, 88, No. 247202. (43) Santara, B.; Giri, P. K.; Imakita, K.; Mujii, M. Evidence of Oxygen Vacancy Induced Room Temperature Ferromagnetism in Solvothermally Synthesized Undoped TiO2 Nanoribbons. Nanoscale 2013, 5, 5476−5488. (44) Coey, J. M. D.; Wongsaprom, K.; Alaria, J.; Venkatesan, M. Charge-Transfer Ferromagnetism in Oxide Nanoparticles. J. Phys. D: Appl. Phys. 2008, 41, No. 134012. (45) Wang, Q.; Sun, Q.; Chen, G.; Kawazoe, Y.; Jena, P. VacancyInduced Magnetism in ZnO Thin Films and Nanowires. Phys. Rev. B 2008, 77, No. 205411. (46) Sato, K.; Dederics, P. H.; Katayama-Yoshida, H. Curie Temperatures of Diluted Magnetic Semiconductors Calculated from First Principles. Europhys. Lett. 2003, 61, 403−408. (47) Bouzerar, G.; Ziman, T.; Kudrnovský, J. Calculating the Curie Temperature Reliably in Diluted III-V Ferromagnetic Semiconductors. Europhys. Lett. 2005, 69, 812−818. (48) Chakraborty, A.; Bouzerar, G. Reconsidering the Possibility of Room Temperature Ferromagnetism in Mn-doped Zirconium Oxide. Europhys. Lett. 2013, 104, 57010−57015. (49) Sato, K.; Bergqvist, L.; Kudrnovský, J.; Dederichs, P. H.; Eriksson, O.; Turek, I.; Sanyal, B.; Bouzerar, G.; Katayama-Yoshida, H.; Dinh, V. A.; Fukushima, T.; Kizaki, H.; Zeller, R. First-Principles Theory of Dilute Magnetic Semiconductors. Rev. Mod. Phys. 2010, 82, 1633. (50) Ruiz-Puigdollers, A.; Illas, F.; Pacchioni, G. Reduction of Hydrogenated ZrO2 Nanoparticles by Water Desorption. ACS Omega 2017, 2, 3878−3885. (51) Kresse, G.; Hafner, J. Ab-initio Molecular Dynamics for Liquid Metals. Phys. Rev. B: Condens. Matter Mater. Phys. 1993, 47, 558−561. (52) Kresse, G.; Hafner, J. Ab-initio Molecular-Dynamics Simulation of the Liquid-Metal Amorphous-Semiconductor Transition in Germanium. Phys. Rev. B: Condens. Matter Mater. Phys. 1994, 49, 14251−14271. (53) Kresse, G.; Furthmuller, J. Efficiency of Ab-initio Total Energy Calculations for Metals and Semiconductors Using a Plane-Wave Basis Set. Comput. Mater. Sci. 1996, 6, 15−50. (54) Blöchl, P. E. Projector Augmented-Wave Method. Phys. Rev. B: Condens. Matter Mater. Phys. 1994, 50, 17953−17979. (55) Kresse, G.; Joubert, J. From Ultrasoft Pseudopotentials to the Projector Augmented-Wave Method. Phys. Rev. B: Condens. Matter Mater. Phys. 1999, 59, 1758−1775. (56) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865−3868. (57) Dudarev, S. L.; Botton, G. A.; Savrasov, S. Y.; Humphreys, C. J.; Sutton, A. P. Electro-Energy-Loss Spectra and the Structural Stability of Nickel Oxide: An LSDA+U Study. Phys. Rev. B: Condens. Matter Mater. Phys. 1998, 57, 1505−1509. (58) Anisimov, V. I.; Aryasetiawan, F.; Lichtenstein, A. I. FirstPrinciples Calculations of the Electronic Structure and Spectra of Strongly Correlated Systems: the LDA+U Method. J. Phys.: Condens. Matter 1997, 9, 767−808. (59) Perdew, J. P.; Burke, K.; Ernzerhof, M. Rationale for Mixing Exact Exchange with Density Functional Approximations. J. Chem. Phys. 1996, 105, 9982−9985. (60) Adamo, C.; Barone, V. Toward Reliable Density Functional Methods without Adjustable Parameters: the PBE0 Model. J. Chem. Phys. 1999, 110, 6158−6170. (61) Dovesi, R.; et al. CRYSTAL14: A Program for the Ab Initio Investigation of Crystalline Solids. Int. J. Quantum Chem. 2014, 114, 1287. (62) Hay, P. J.; Wadt, W. R. Ab Initio Effective Core Potentials for Molecular Calculations. Potentials for K to Au Including the Outermost Core Orbitals. J. Chem. Phys. 1985, 82, 299−310. 5307
DOI: 10.1021/acsomega.8b00667 ACS Omega 2018, 3, 5301−5307
