Functional Mn–Mg k cation complexes in GaN featured by Raman spectroscopy T. Devillers, D. M. G. Leite, J. H. Dias da Silva, and A. Bonanni Citation: Applied Physics Letters 103, 211909 (2013); doi: 10.1063/1.4833024 View online: http://dx.doi.org/10.1063/1.4833024 View Table of Contents: http://scitation.aip.org/content/aip/journal/apl/103/21?ver=pdfcov Published by the AIP Publishing Articles you may be interested in Manganese-induced long-range lattice disorder and vacancy formation in metal-organic chemical vapor deposition grown and ion-implanted Ga 1 − x Mn x N J. Vac. Sci. Technol. A 24, 1640 (2006); 10.1116/1.2201052 Local vibrational modes of H complexes in Mg-doped GaN grown by molecular beam epitaxy Appl. Phys. Lett. 84, 897 (2004); 10.1063/1.1645668 Raman study of Mg, Si, O, and N implanted GaN J. Appl. Phys. 94, 4389 (2003); 10.1063/1.1606521 Improvement of magnetic property of GaMnN by codoping of Mg J. Appl. Phys. 93, 6793 (2003); 10.1063/1.1556114 Mg doping of GaN layers grown by plasma-assisted molecular-beam epitaxy Appl. Phys. Lett. 76, 718 (2000); 10.1063/1.125872 This article is copyrighted as indicated in the article. Reuse of AIP content is subject to the terms at: http://scitation.aip.org/termsconditions. Downloaded to IP: 186.217.234.103 On: Fri, 10 Oct 2014 16:27:05 http://scitation.aip.org/content/aip/journal/apl?ver=pdfcov http://oasc12039.247realmedia.com/RealMedia/ads/click_lx.ads/www.aip.org/pt/adcenter/pdfcover_test/L-37/1077920811/x01/AIP-PT/APL_ArticleDL_1014/AIP-2293_Chaos_Call_for_EIC_1640x440.jpg/47344656396c504a5a37344142416b75?x http://scitation.aip.org/search?value1=T.+Devillers&option1=author http://scitation.aip.org/search?value1=D.+M.+G.+Leite&option1=author http://scitation.aip.org/search?value1=J.+H.+Dias+da+Silva&option1=author http://scitation.aip.org/search?value1=A.+Bonanni&option1=author http://scitation.aip.org/content/aip/journal/apl?ver=pdfcov http://dx.doi.org/10.1063/1.4833024 http://scitation.aip.org/content/aip/journal/apl/103/21?ver=pdfcov http://scitation.aip.org/content/aip?ver=pdfcov http://scitation.aip.org/content/avs/journal/jvsta/24/4/10.1116/1.2201052?ver=pdfcov http://scitation.aip.org/content/avs/journal/jvsta/24/4/10.1116/1.2201052?ver=pdfcov http://scitation.aip.org/content/aip/journal/apl/84/6/10.1063/1.1645668?ver=pdfcov http://scitation.aip.org/content/aip/journal/jap/94/7/10.1063/1.1606521?ver=pdfcov http://scitation.aip.org/content/aip/journal/jap/93/10/10.1063/1.1556114?ver=pdfcov http://scitation.aip.org/content/aip/journal/apl/76/6/10.1063/1.125872?ver=pdfcov Functional Mn–Mgk cation complexes in GaN featured by Raman spectroscopy T. Devillers,1,a) D. M. G. Leite,2,3 J. H. Dias da Silva,3 and A. Bonanni1,b) 1Institut f€ur Halbleiter-und-Festk€orperphysik, Johannes Kepler University, Altenbergerstr. 69, A-4040 Linz, Austria 2Instituto de F�ısica e Qu�ımica, Universidade Federal de Itajub�a, 37500-903, Itajub�a–MG, Brazil 3Department of Physics, S~ao Paulo State University, Bauru–SP, Brazil (Received 23 October 2013; accepted 8 November 2013; published online 22 November 2013) The evolution of the optical branch in the Raman spectra of (Ga,Mn)N:Mg epitaxial layers as a function of the Mn and Mg concentrations, reveals the interplay between the two dopants. We demonstrate that the various Mn-Mg-induced vibrational modes can be understood in the picture of functional Mn–Mgk complexes formed when substitutional Mn cations are bound to k substitutional Mg through nitrogen atoms, the number of ligands k being driven by the ratio between the Mg and the Mn concentrations. VC 2013 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution 3.0 Unported License. [http://dx.doi.org/10.1063/1.4833024] Nitride semiconductors have increasingly attracted attention over the last decades, particularly due to their remarkably broad spectrum of technologically relevant applications.1 While nitride-based compounds undoubtedly prevail as building-bocks for visible and ultraviolet solid- state optoelectronic and high power electronic devices, the fields of spintronics,2 thermoelectricity,3 and terahertz elec- tronics4 are among the ones where nitrides are currently giv- ing rise to a growing interest for prospective functionalities. The realization of the next generation of nitride-based devi- ces requires a profound understanding and control of the doping mechanisms in these materials. In this context, it is becoming increasingly clear that the traditional picture of isolated impurities is not sufficient to properly illustrate the properties of these materials systems, particularly when two or more dopant species are involved and/or in regime of high concentration. We have recently demonstrated how the for- mation of Mn–Mgk complexes in (Ga,Mn)N:Mg substan- tially affects the behaviour of each dopant, resulting in remarkable and unforeseen material properties.5 The inter- play between the dopants, but also with the host lattice must be understood in order to ensure a controlled and effective functionality of these nano-objects in innovative devices. Raman spectroscopy is particularly suitable to investigate impurities and impurity complexes in semiconductors6 due to its sensitivity to variations in the structural and electronic properties of the material under study. In the case of doping of nitrides with Mg, this technique was decisive to identify the Mg–H complexes7,8 responsible for the low efficiency of the Mg activation, and for the consequent poor p-type con- ductivity especially in GaN:Mg.9–11 While Mg is the most relevant p-type dopant in GaN, and is employed in the ma- jority of optoelectronic devices based on GaN, Mn generated interest lately for its potential in spintronics, and particularly in the perspective of turning p-type Mn-doped GaN (i.e., (Ga,Mn)N:Mg) into a ferromagnetic dilute magnetic semiconductor with a high Curie transition temperature.12 Moreover, the recent demonstration of unexpected (infrared) optical and magnetic responses in (Ga,Mn)N:Mg,5 has opened appealing fundamental and technological perspec- tives for this system. In this Letter, we apply Raman spectroscopy at room temperature in order to investigate the local environment of Mn and Mg in (Ga,Mn)N:Mg. The samples consist of a 1 lm thick GaN buffer layer grown on c-sapphire (0001), on which a 600 nm layer of GaN doped with Mn and/or Mg is depos- ited. The samples are fabricated by metalorganic vapor phase epitaxy (MOVPE) according to a procedure detailed else- where.5,13 Prior to Raman spectroscopy, the structural, mag- netic, and chemical properties of the samples are extensively investigated5 and in particular the crystalline structure is ana- lysed by synchrotron high-resolution x-ray diffraction and high-resolution transmission electron microscopy, showing that all samples are single-crystalline and do not include sec- ondary phases. Furthermore, synchrotron x-ray absorption fine structure (EXAFS) studies demonstrate that at least 95% of the Mn ions substitutes Ga cations,14 and we have reported that in the case of Mn-Mg codoping, both Mn and Mg occupy Ga substitutional sites, and Mg tends to stay close to Mn, sep- arated by one N atom, to form Mn–Mgk cation complexes, where k is the number of Mg ions in the first cation shell of a single Mn.5 The Mn and Mg concentrations xMn and xMg are extracted from secondary ion mass spectroscopy (SIMS) and lie between 0% and 1%. The Raman spectra are acquired with a Jobin-Yvon LabRAM HTS confocal micro-Raman setup equipped with 532 nm, 633 nm, and 785 nm lasers with a spectral resolution of 1 cm�1. All spectra presented here are acquired using the 532 nm beam, in ZðX�Þ �Z configuration, i.e., at normal incidence and emergence, with a linearly polar- ized incident beam, and unpolarized scattered beam. After the subtraction of a linear background, all the spectra are renormalized with respect to the intensity of the EH 2 peak. Particular attention is given to the optical branch of GaN, especially to the energy range between 500 cm�1 and 800 cm�1. a)thibaut.devillers@jku.at b)alberta.bonanni@jku.at 0003-6951/2013/103(21)/211909/4 VC Author(s) 2013103, 211909-1 APPLIED PHYSICS LETTERS 103, 211909 (2013) This article is copyrighted as indicated in the article. Reuse of AIP content is subject to the terms at: http://scitation.aip.org/termsconditions. Downloaded to IP: 186.217.234.103 On: Fri, 10 Oct 2014 16:27:05 http://dx.doi.org/10.1063/1.4833024 http://dx.doi.org/10.1063/1.4833024 http://dx.doi.org/10.1063/1.4833024 http://dx.doi.org/10.1063/1.4833024 mailto:thibaut.devillers@jku.at mailto:alberta.bonanni@jku.at http://crossmark.crossref.org/dialog/?doi=10.1063/1.4833024&domain=pdf&date_stamp=2013-11-22 In Fig. 1, the Raman spectra for four representative sam- ples are reported, namely GaN, GaN:Mn (xMn¼ 0.75%), GaN:Mg (xMg¼ 0.25%), and GaN:(Mn,Mg) (xMn¼ 0.25%, xMg¼ 0.45%). In the present geometry, the optical branch of the GaN spectra is dominated by two modes: the transverse EH 2 mode at 568 cm�1 and the longitudinal A1(LO) at 733 cm�1. The A1(TO) is not visible, and only a small residual contribution of the E1(TO) can be observed, in the form of a component at the low energy side of EH 2 . The very low intensity of this E1(TO) is a confirmation of the good crystalline quality of the layers. The influence of Mg doping on the Raman spectra of GaN:Mg has already been widely discussed in literature.9–11,15 In the studied energy range, only the local vibrational mode (LVM) at 657 cm�1 for substitutional Mg in GaN is expected. However, as-grown samples generally do not exhibit this mode, since H is bound to Mg in Mg–N–H complexes, which are Raman-silent in this range.16 An annealing step is necessary in order to break the complex, to allow H to diffuse towards the sample surface, and to finally reveal the free-Mg LVM.17–19 As expected and as seen in Fig. 1, the as-grown GaN:Mg considered here does not exhibit the Mg-LVM, confirming that H is bound to Mg. The introduction of Mn into GaN induces a much richer structure in the spectra and the resulting Raman signature was previously discussed, though the origin of the various contributions is still a matter of controversy.20–22 To analyse the spectra of (Ga,Mn)N, we employ an approach similar to the one of Gebicki et al.,23 but based on the deconvolution of the Mn-related signal into six Lorentzian contributions, as shown in Fig. 2. In addition to the GaN EH 2 , E1(TO), and A1(LO), we identify a contribution labeled L which corresponds to the Mn LVM associated to the GaN op- tical branch. According to the light impurity model,23–25 the frequency xTO Mn�LVM of the TO Mn-LVM can be estimated via the formula xTO Mn-LVM ¼ xTO GaN ffiffiffiffiffiffiffiffiffiffi lGaN lMnN r ; where lGa�N and lMn�N correspond to the reduced mass of the Ga–N and M–N pairs, respectively.25 Considering the frequency of the EH 2 peak (569 cm�1)—which is the only transverse optical mode allowed in our geometry—we obtain for the TO Mn-LVM peak a theoretical value of 582 cm�1, which is in agreement with the experimental value of 579 cm�1. In addition to the Mn-LVM, we can also identify five lorentzian contributions (A,B,C,D,E). The integrated in- tensity of the C peak depends linearly on the Mn concentra- tion, and the whole set (ABCDE) vanishes when increasing the laser wavelength from 633 nm to 785 nm (not shown), confirming the resonant character of this broad band hinted at in Refs. 21 and 23. In the samples containing both Mn and Mg, the presence of Mg heavily affects the Mn-related signature in the Raman spectrum. In Fig. 3, the Raman spectra for (Ga,Mn)N:Mg samples with different Mn and Mg concentrations ranging from 0.1% up to 0.5%, are reported. One can identify four re- markable features occurring in the presence of both Mn and Mg, namely: (1) a significant quenching of the Mn-induced FIG. 1. Raman spectra of: nominally undoped GaN, GaN:Mg, (Ga,Mn)N, and (Ga,Mn)N:Mg. To clarity, an y-offset of 0.05 has been added between the curves. FIG. 2. Deconvolution of the Mn-induced Raman broad band in multiple Lorentzians. FIG. 3. Raman spectra of (Ga,Mn)N:Mg as a function of the Mn and Mg concentrations. The curves are grouped by three: in each group the Mn concentration is kept constant while the Mg concentration is varied. The groups are separated by a y-offset of 0.05. The characteristic features of (Ga,Mn)N:Mg are indicated by (1), (2), (3), and (4). 211909-2 Devillers et al. Appl. Phys. Lett. 103, 211909 (2013) This article is copyrighted as indicated in the article. Reuse of AIP content is subject to the terms at: http://scitation.aip.org/termsconditions. Downloaded to IP: 186.217.234.103 On: Fri, 10 Oct 2014 16:27:05 broad band (ABDE and to a lesser extent of C); (2) a shift in the position of the C peak; (3) the appearance of an addi- tional sharp peak around 650 cm�1; (4) the presence of a spe- cific broad peak around 688 cm�1. In the following, we discuss these features in the frame of the recently demon- strated evidence of Mn–Mgk complexes in (Ga,Mn)N:Mg.5 In order to provide a quantitative description, we deconvo- lute the spectra in a similar way as we did in the case of (Ga,Mn)N (without Mg). To take the influence of Mg into account, three Lorentzian contributions are added, one around 650 cm�1 to describe feature (3) and two between 680 cm�1 and 695 cm�1 to describe feature (4). With these two components, it is possible to reproduce the experimental data, and to extrapolate the relevant parameters to describe the evolution of the spectra as a function of Mn and Mg, respectively. The Mn-induced broad band covering the range between EH 2 and A1(LO) (peaks A,B,D,E) is highly depend- ent on the Mg concentration. The intensity of this broad band is at its maximum in the absence of Mg and is rapidly damped as soon as Mg is added, until it goes under the detection limit for xMg � 2xMn. The fact that the attenuation of the C peak is in comparison much less evident—also upon renormalized with the Mn concentration—is a hint of the fact that, despite the symmetry common to all these peaks as demonstrated by Gebicki et al.,23 the physical mechanism behind the origin of the C peak is different from the one giving rise to A,B,D,E. The C peak, i.e., the main feature correlated with Mn, and which scales with the Mn concentration, is affected by the Mg incorporation both in energy position and in inten- sity. We can observe that the intensity variation is mostly due to the variations in the Mn concentration, but the shift in energy is clearly correlated to the introduction of Mg, and does not directly depend on the Mn concentration. The gen- eral trend is that the increment of the ratio xMg/xMn between the Mg and Mn concentrations induces a blue shift of the C peak. In Fig. 4(a) we report the shift of the C peak as a func- tion of the ratio between the concentrations of Mg and Mn. Remarkably, for xMg � xMn, we observe that this shift and the evolution of the Mn–Mg coordination5 follow the same trend. The shift in the C peak is correlated to the formation of Mn–Mgk complexes as the Mg concentration increases. Previous ab initio calculations5 have evidenced that the incorporation of Mg in the first cation coordination shell of Mn induces a displacement of the N atom binding Mg to Mn towards the Mn atom. The increment in the stretching fre- quency with the Mg/Mn ratio is completely coherent with the shortening of the Mn–N bonds induced by the presence of Mg. The saturation of these shifts reveals that for a high Mg/Mn ratio, the Mn environment is not perturbed anymore, as shown by synchrotron EXAFS measurements.5 However, this simplified model does not account for the red shift of the C peak at low xMg /xMn. To understand this phenomenon, it is necessary to take into account the spin density calculations reported in Ref. 5. The simple Mn–Mg1 complex tends to delocalize the spin density on the nitrogen atoms, while the complexes of higher order (Mn–Mg2,…) reduce this delocal- ization of the spin density and enhance the electron-phonon coupling. This effect is likely to be responsible for the non-monotonous behaviour of the C peak position as a func- tion of xMg/xMn. The peak at 650 cm�1 is clearly visible in all samples containing both Mn and Mg. From its energy position, this feature can reasonably be assigned to the Mg local vibration mode (Mg-LVM) which is generally reported at 657 cm�1 in activated Mg-doped GaN. However, this LVM is usually not visible in as-grown GaN:Mg, due to the fact that the for- mation of a Mg–N–H complex significantly distorts the MgGa–N4 tetrahedron.10,26 In our case, the presence of this peak in as-grown samples can be explained taking into account the atomic arrangement of Mg and Mn in GaN. As shown by SIMS measurements,5 the presence of Mn in (Ga,Mn)N:Mg induces the H concentration to be between one and two orders of magnitude lower than in GaN:Mg. The incorporation of H is actually hindered by the formation of the Mn–Mgk complexes, since Mg is electrically passivated by the proximity of a Mn ion. On the other hand, the presence of Mn in the vicinity of Mg shifts N towards Mn, increasing in this way the Mg–N bond length from 2.03 Å to 2.11 Å. This distortion of the Mg–N4 tetrahedron is much smaller FIG. 4. (a) Position of the C peak as a function of the xMg/xMn ratio; (b,c,d,e) Structure of Mn–Mgk complexes for k¼ (0,1,2,3), corresponding to the ratio xMg/xMn¼ (0,1,2,3); (f) Integrated intensity of the Mn–Mgk-related double peak (680 cm�1–695 cm�1) as a function of the xMg /xMn ratio. 211909-3 Devillers et al. Appl. Phys. Lett. 103, 211909 (2013) This article is copyrighted as indicated in the article. Reuse of AIP content is subject to the terms at: http://scitation.aip.org/termsconditions. Downloaded to IP: 186.217.234.103 On: Fri, 10 Oct 2014 16:27:05 than the one induced by the formation of a Mg–N–H com- plex26 and does not impair the Mg-LVM. However, the red-shift observed (peak maximum at 650 cm�1 instead than at 657 cm�1) is likely to originate from this distortion. This red-shift does not depend on the Mg concentration, in accord- ance with the fact that the environment of the Mg atom does not depend massively on the Mn and Mg concentrations (in- dependently of the Mn and Mg concentrations, each Mg is bound to one and only one Mn). The last feature induced by the presence of Mg in (Ga,Mn)N:Mg is a relatively broad band around 688 cm�1, which—particularly at high xMg /xMn—is actually composed of at least two peaks, between 680 cm�1 and 695 cm�1. To describe phenomenologically this mode, we report in Fig. 4(f) the integrated intensity of this band (the sum of the inte- grated intensity of the two peaks), normalized by the Mn concentration, as a function of the ratio between the Mg and Mn concentrations. If we compare the trend in Fig. 4(f) to the populations of the different complexes calculated in Ref. 5, it appears that all three complexes give a contribution, and particularly the Mn–Mg2 complex. Surprisingly, according to previous ab initio calculations, this complex is the one with the lowest symmetry, the Mn–N4 tetrahedron being more distorted than for other complexes. Since a purely structural argument cannot account for the Raman activity of the Mn–Mg2 complex, one should also consider the influence of the electronic structure. In addition to differences in the localization of carriers for different complexes that we al- ready mentioned, the presence of free Mg in GaN for xMg=xMn � 3 and the consequent introduction of holes is also expected to affect the Raman activity of this Mn-Mgk-induced Raman band, and to be responsible for the quenching at high xMg/xMn. In conclusion, Raman spectroscopy measurements at room temperature of (Ga,Mn)N:Mg have shown that the codoping of GaN with Mn and Mg generates Mn–Mgk cation complexes, in agreement with previously reported synchro- tron EXAFS measurements. Moreover, the evolution of the different features relative to Mg and Mn in (Ga,Mn)N:Mg supports the ab initio calculations of the local arrangements around the complexes, in particular the displacement of the N atoms surrounding Mn. This information, contributes to the understanding and control of the functional complexes Mn–Mgk that have been shown to introduce into the system two quite spectacular features, namely: (i) a broadband infra- red luminescence persisting at room temperature and (ii) the possibility to tune the spin state of Mn in a controlled way as a function of the degree of coordination k, paving the way to the realization of optically active and tunable spintronic devices. This work was supported by the FunDMS Advanced Grant of the European Research Council (ERC Grant No. 227690) within the Ideas 7th Framework Programme of the European Community, by the Austrian Science Fundation–FWF (P22471, P20065, and P22477) and by the Fapesp (Grant Nos. 2006/05627-8 and 2012/21147-7) 1H. Morkoç, Handbook of Nitride Semiconductors (Wiley-VCH Verlag GmbH, 2008), Vol. 3. 2T. Dietl, Nature Mater. 9, 965 (2010). 3N. Lu and I. Ferguson, Semicond. Sci. Technol. 28, 074023 (2013). 4N. Nepal, V. D. Wheeler, T. J. Anderson, F. J. Kub, M. A. Mastro, R. L. Myers-Ward, S. B. Qadri, J. A. Freitas, S. C. Hernandez, L. O. Nyakiti et al., Appl. Phys. Express 6, 061003 (2013). 5T. Devillers, M. Rovezzi, N. G. Szwacki, S. Dobkowska, W. Stefanowicz, D. Sztenkiel, A. Grois, J. Suffczy�nski, A. Navarro-Quezada, B. Faina et al., Sci. Rep. 2, 722 (2012). 6J. Wagner, Appl. Surf. Sci. 50, 79 (1991). 7W. G€otz, N. M. Johnson, D. P. Bour, M. D. McCluskey, and E. E. Haller, Appl. Phys. Lett. 69, 3725 (1996). 8R. Cusco, N. Domenech-Amador, L. Artus, K. Wank, T. Yamaguchi, and Y. Nanishi, J. Appl. Phys. 112, 053528 (2012). 9M. S. Brandt, W. Gotzt, and N. M. Johnson, Phys. Rev. B 49, 14758 (1994). 10H. Harima, T. Inoue, S. Nakashima, M. Ishida, and M. Taneya, Appl. Phys. Lett. 75, 1383 (1999). 11A. Kaschner, H. Siegle, G. Kaczmarczyk, M. Straßburg, A. Hoffmann, and C. Thomsen, Appl. Phys. Lett. 74, 3281 (1999). 12T. Dietl, H. Ohno, and F. Matsukura, Phys. Rev. B 63, 195205 (2001). 13A. Bonanni, M. Sawicki, T. Devillers, W. Stefanowicz, B. Faina, T. Li, T. E. Winkler, D. Sztenkiel, A. Navarro-Quezada, M. Rovezzi et al., Phys. Rev. B 84, 035206 (2011). 14W. Stefanowicz, D. Sztenkiel, B. Faina, A. Grois, M. Rovezzi, T. Devillers, F. D’Acapito, A. Navarro-Quezada, T. Li, R. Jakieła et al., Phys. Rev. B 81, 235210 (2010). 15R. Kirste, M. P. Hoffmann, J. Tweedie, Z. Bryan, G. Callsen, T. Kure, C. Nenstiel, M. R. Wagner, R. Collazo, A. Hoffmann et al., J. Appl. Phys. 113, 103504 (2013). 16F. J. Manjon, B. Mari, J. Serrano, and A. H. Romero, J. Appl. Phys. 97, 053516 (2005). 17S. Nakamura, N. Iwasa, M. Senoh, and T. Tokushima, Jpn. J. Appl. Phys., Part 1 31, 1258 (1992). 18W. G€otz, N. M. Johnson, D. P. Bour, and R. A. Street, Appl. Phys. Lett. 68, 667 (1996). 19G. Pozina, C. Hemmingsson, P. P. Paskov, J. P. Bergman, B. Monemar, T. Kawashima, H. Amano, I. Akasaki, and A. Usui, Appl. Phys. Lett. 92, 151904 (2008). 20L. L. Guo, Y. H. Zhang, and W. Z. Shen, Appl. Phys. Lett. 89, 161920 (2006). 21J. H. Dias da Silva, D. M. G. Leite, and A. R. Zanatta, J. Phys. Condens. Matter 20, 252201 (2008). 22E. Alarcon-Llado, J. Iba~nez, R. Cusco, L. Artus, S. V. Novikov, and T. C. Foxon, Semicond. Sci. Technol. 24, 115019 (2009). 23W. Gebicki, P. Dominik, and S. Podsiadło, Phys. Rev. B 77, 245213 (2008). 24W. Gebicki, J. Strzeszewski, G. Kamler, T. Szyszko, and S. Podsiado, Appl. Phys. Lett. 76, 3870 (2000). 25A. S. Barker and A. J. Sievers, Rev. Mod. Phys. 47, S1 (1975). 26S. Limpijumnong, J. E. Northrup, and C. G. Van de Walle, Phys. Rev. B 68, 075206 (2003). 211909-4 Devillers et al. Appl. Phys. Lett. 103, 211909 (2013) This article is copyrighted as indicated in the article. Reuse of AIP content is subject to the terms at: http://scitation.aip.org/termsconditions. Downloaded to IP: 186.217.234.103 On: Fri, 10 Oct 2014 16:27:05 http://dx.doi.org/10.1038/nmat2898 http://dx.doi.org/10.1088/0268-1242/28/7/074023 http://dx.doi.org/10.7567/APEX.6.061003 http://dx.doi.org/10.1038/srep00722 http://dx.doi.org/10.1016/0169-4332(91)90142-7 http://dx.doi.org/10.1063/1.117202 http://dx.doi.org/10.1063/1.4749266 http://dx.doi.org/10.1103/PhysRevB.49.14758 http://dx.doi.org/10.1063/1.124701 http://dx.doi.org/10.1063/1.124701 http://dx.doi.org/10.1063/1.123320 http://dx.doi.org/10.1103/PhysRevB.63.195205 http://dx.doi.org/10.1103/PhysRevB.84.035206 http://dx.doi.org/10.1103/PhysRevB.84.035206 http://dx.doi.org/10.1103/PhysRevB.81.235210 http://dx.doi.org/10.1063/1.4794094 http://dx.doi.org/10.1063/1.1856222 http://dx.doi.org/10.1143/JJAP.31.1258 http://dx.doi.org/10.1143/JJAP.31.1258 http://dx.doi.org/10.1063/1.116503 http://dx.doi.org/10.1063/1.2909541 http://dx.doi.org/10.1063/1.2364472 http://dx.doi.org/10.1088/0953-8984/20/25/252201 http://dx.doi.org/10.1088/0953-8984/20/25/252201 http://dx.doi.org/10.1088/0268-1242/24/11/115019 http://dx.doi.org/10.1103/PhysRevB.77.245213 http://dx.doi.org/10.1063/1.126804 http://dx.doi.org/10.1103/RevModPhys.47.S1.2 http://dx.doi.org/10.1103/PhysRevB.68.075206