DOE Final Report on Grant DE-FG02-03ER46063 “The Controlled Synthesis of Metastable Oxides Utilizing Epitaxy and Epitaxial Stabilization” August 2003 - February 2007
Principal Investigator: Darrell G. Schlom Period Covered by this Final Report: August 2003 - February 2007 DOE Award Number: DE-FG02-03ER46063 Amount of Unexpended Funds: $0 The research enabled by this DOE grant led to 13 publications in leading refereed journals including Physical Review Letters and Applied Physics Letters as well as feature articles in the Journal of the American Ceramic Society and Physik Journal (the German equivalent of Physics Today distributed to all members of the German Physical Society) on the controlled synthesis of metastable oxides utilizing epitaxy and epitaxial stabilization. In total our results fill over 100 pages of archived journals and are attached. Much of our work has been performed in collaboration with other DOE researchers. Specifically, we have active collaborations with (please see publications on the next page) with the following DOE groups.
DOE COLLABORATORS ON THE WORK SUPPORTED BY THIS DOE GRANT
DOE Laboratory Collaborators Argonne National Laboratory Stephen K. Streiffer, Mark A. Zurbuchen, and Dillon D. Fong Pacific Northwest National Laboratory Scott A. Chambers Los Alamos National Laboratory Marilyn E. Hawley
PERSONNEL SUPPORTED BY THIS DOE GRANT
Name Position Darrell G. Schlom Principal Investigator Wei Tian Postdoctoral Researcher
PUBLICATIONS AND PRESENTATIONS The support from this DOE program has led to the publications listed below and attached to the end of this document. Publications 1. M.A. Zurbuchen, Y. Jia, S. Knapp, A.H. Carim, D.G. Schlom, and X.Q. Pan, “Defect Generation
by Preferred Nucleation in Epitaxial Sr2RuO4 / LaAlO3,” Applied Physics Letters 83 (2003) 3891-3893.
1 2. A.V. Pogrebnyakov, J.M. Redwing, S. Raghavan, V. Vaithyanathan, D.G. Schlom, S.Y. Xu, Q. Li, D.A. Tenne, A. Soukiassian, X.X. Xi, M.D. Johannes, D. Kasinathan, W.E. Pickett, J.S. Wu,
and J.C.H. Spence, “Enhancement of Superconducting Transition Temperature of MgB2 by a Strain-Induced Bond-Stretching Mode Softening,” Physical Review Letters 93 (2004) 147006-1 - 147006-4. 3. A.V. Pogrebnyakov, X.X. Xi, J.M. Redwing, V. Vaithyanathan, D.G. Schlom, A. Soukiassian, S.B. Mi, C.L. Jia, J.E. Giencke, C.B. Eom, J. Chen, Y.F. Hu, Y. Cui, and Q. Li, “Properties of
MgB2 Thin Films with Carbon Doping,” Applied Physics Letters 85 (2004) 2017-2019. 4. D.A. Tenne, X.X. Xi, Y.L. Li, L.Q. Chen, A. Soukiassian, M.H. Zhu, A.R. James, J. Lettieri, D.G. Schlom, W. Tian, and X.Q. Pan, “Absence of Low-Temperature Phase Transitions in
Epitaxial BaTiO3 Thin Films,” Physical Review B 69 (2004) 174101-1 - 174101-5. 5. X.X. Xi, A.V. Pogrebnyakov, X.H. Zeng, J.M. Redwing, S.Y. Xu, Q. Li, Z.K. Liu, J. Lettieri, V. Vaithyanathan, D.G. Schlom, H.M. Christen, H.Y. Zhai, and A. Goyal, “Progress in the
Deposition of MgB2 Thin Films,” Superconductor Science and Technology 17 (2004) S196-S201 and in: Applied Superconductivity 2003: Proceedings of the Sixth European Conference on Applied Superconductivity, IOP Vol. 181 (Institute of Physics, Bristol, 2004), pp. 37-44. 6. J. Mannhart and D.G. Schlom, “Oxide—Tausendsassa für die Elektronik,” Physik Journal 4 (2005) 45-51. 7. A. Venimadhav, A. Soukiassian, D.A. Tenne, Q. Li, X.X. Xi, D.G. Schlom, R. Arroyave, Z.K. Liu, H.P. Sun, X.Q. Pan, M.Y. Lee, and N.P. Ong, “Structural and Transport Properties of
Epitaxial NaxCoO2 Thin Films,” Applied Physics Letters 87 (2005) 172104-1 – 172104-3. 8. P. Orgiani, Y. Cui, A.V. Pogrebnyakov, J.M. Redwing, V. Vaithyanathan, D.G. Schlom, and
X.X. Xi, “Investigations of MgB2 / MgO and MgB2 / AlN Heterostructures for Josephson Devices,” IEEE Transactions on Applied Superconductivity 15 (2005) 228-231.
9. B.T. Liu, X.X. Xi, V. Vaithyanathan, and D.G. Schlom, “MgB2 Thin Films Grown at Different Temperatures by Hybrid Physical-Chemical Vapor Deposition,” IEEE Transactions on Applied Superconductivity 15 (2005) 3249-3252. 10. M.A. Zurbuchen, J. Lettieri, Y. Jia, A.H. Carim, S.K. Streiffer, and D.G. Schlom, “Out-of-phase Boundary (OPB) Nucleation in Layered Oxides,” in: Ferroelectric Thin Films XIII edited by R. Ramesh, J.-P. Maria, M. Alexe, and V. Joshi, Vol. 902E (Materials Research Society, Warrendale, 2006), pp. T10-55.1 – T10-55.6. 11. W. Tian, J.H. Haeni, D.G. Schlom, E. Hutchinson, B.L. Sheu, M.M. Rosario, P. Schiffer, Y. Liu, M.A. Zurbuchen, and X.Q. Pan, “Epitaxial Growth and Magnetic Properties of the First Five
Members of the Layered Srn+1RunO3n+1 Oxide Series,” Applied Physics Letters 90 (2007) 022507- 1 – 022507-3. 12. M.A. Zurbuchen, W. Tian, X.Q. Pan, D. Fong, S.K. Streiffer, M.E. Hawley, J. Lettieri, Y. Jia, G. Asayama, S.J. Fulk, D.J. Comstock, S. Knapp, A.H. Carim, and D.G. Schlom, “Morphology, Structure, and Nucleation of Out-of-Phase Boundaries (OPBs) in Epitaxial Films of Layered Oxides,” Journal of Materials Research 22 (2007) 1439-1471. 13. M.A. Zurbuchen, R.S. Freitas, M.J. Wilson, P. Schiffer, M. Roeckerath, J. Schubert, M.D. Biegalski, G.H. Mehta, D.J. Comstock, J.H. Lee, Y. Jia, and D.G. Schlom, “Synthesis and Characterization of an n = 6 Aurivillius Phase Incorporating Magnetically Active Manganese,
Bi7(Mn,Ti)6O21,” Applied Physics Letters 91 (2007) 033113-1 – 033113-3. 14. D.G. Schlom, L.Q. Chen, X.Q. Pan, A. Schmehl, and M.A. Zurbuchen, “A Thin Film Approach to Engineering Functionality into Oxides,” Journal of the American Ceramic Society 91 (2008) 2429-2454. (Feature Article with Journal Cover)
2 APPLIED PHYSICS LETTERS VOLUME 83, NUMBER 19 10 NOVEMBER 2003 Õ Defect generation by preferred nucleation in epitaxial Sr2RuO4 LaAlO3 Mark A. Zurbuchen,a) Yunfa Jia, Stacy Knapp, Altaf H. Carim,b) and Darrell G. Schlom Department of Materials Science and Engineering, The Pennsylvania State University, University Park, Pennsylvania 16803-6602 X. Q. Pan Department of Materials Science and Engineering, University of Michigan, Ann Arbor, Michigan 48109-2136 ͑Received 27 May 2003; accepted 12 September 2003͒ ͑ ͒ The atomic structure of the film–substrate interface of a 001 Sr2RuO4/(100)c LaAlO3 film, determined by high-resolution transmission electron microscopy and simulation, is reported. The structure of superconductivity-quenching ⌬cϷ0.25 nm out-of-phase boundaries ͑OPBs͒ in the film is also reported. Growth in one region on the La-terminated surface is observed to nucleate with a SrO layer. Because two structurally equivalent SrO layers exist within the unit cell, two neighboring nuclei with differing growth order (SrO-RuO2-SrO or RuO2-SrO-SrO) will nucleate an OPB where their misaligned growth fronts meet. Strategies to avoid OPB generation by this mechanism are suggested, which it is hoped may ultimately lead to superconducting Sr2RuO4 films. © 2003 American Institute of Physics. ͓DOI: 10.1063/1.1624631͔
Sr2RuO4 is the only known copper-free layered perov- In the heteroepitaxial growth of layered structures, sub- skite superconductor,1 and is believed to have an unconven- strate surface quality has a great effect on final film quality.14 tional spin-triplet (p-wave͒ symmetry of the superconduct- Advances in substrate engineering have enabled the prepara- 2 14,15 ing order parameter. Electrical transport studies of Sr2RuO4 tion of substrates with known atomic terminations. Some are expected to further the understanding of the super- disagreement currently exists regarding the atomic termina- 3 ͑ ͒ 16–19 conducting state. Sr2RuO4 is well lattice-matched to tion of 001 LaAlO3 . YBa2Cu3O7Ϫ␦ and other cuprate superconductors. Epitaxial It would be useful experimentally to determine which heterostructures would be useful for probing the interaction atomic layers are energetically preferred for the nucleation of between d-wave and p-wave superconductors. Although the Sr2RuO4 on a given substrate, so that high-quality films can 4,5 ͓ epitaxial growth of Sr2RuO4 films and epitaxial be grown by sequential deposition techniques e.g., 6 ͑ ͔͒ YBa2Cu3O7Ϫ␦ /Sr2RuO4 heterostructures have been dem- molecular-beam epitaxy MBE . Few studies of the pre- onstrated, no superconducting Sr2RuO4 films have been ferred nucleation layer of epitaxial thin films of complex reported.4,5 As we show in this letter, greater attention to oxides on substrates with known termination have been nucleation and growth is important to reducing defects and to reported.20–22 Understanding defect nucleation mechanisms possibly achieving superconductivity in Sr2RuO4 films. and microstructural development in the growth of complex Both impurities and crystallographic defects can quench oxide films will allow greater control over their growth. Sub- superconductivity in Sr2RuO4 . As little as 300 ppm of alu- strate surface features have been demonstrated to result in minum is sufficient to destroy superconductivity.7 Crystallo- micron-scale crystallographic defects in some similarly graphic defects that quench superconductivity include out- structured layered perovskites (nϭ2 Aurivillius phases͒.23,24 of-phase boundaries ͑OPBs͒.8 These are translation In this letter, the structure of the film–substrate interface of ͑ ͒ 16 boundaries consisting of a fractional misalignment in the an epitaxial 001 Sr2RuO4 /(001)c LaAlO3 film is re- long axis (c-axis͒ direction between two neighboring regions ported, and the mechanism whereby the preferred nucleation of the same crystal. Because any interruption of the structure layer can lead to the generation of superconductivity- can act as a pair-breaker, a linear density of OPBs on the quenching OPBs is described. ͑ ͒ order of 1/ ab(0), where ab(0) is the superconducting co- (001)c Sr2RuO4 epitaxial films were grown on 001 ͑ ͒ Ϸ 9 16 ͑ ͒ herence length in the 001 plane, ab(0) 66 nm, renders a LaAlO3 substrates by pulsed-laser deposition PLD at high-purity film nonsuperconducting.9 A correlation between 1000 °C in a radiatively heated sample chamber. Specifics of 8 Tc and the residual resistivity ( 0), a measure of disorder in the film growth are described in detail elsewhere. The film equivalent-purity crystals, has also been reported.10 Although described here was a high-quality film, judged by the full OPBs are common defects in complex oxides, little is known width at half-maximum of the 006 Sr2RuO4 peak of approxi- about their nucleation mechanisms or impact upon mately 0.23° in 2, and 0.33° in ͑instrumental resolution properties.11–13 Ϸ0.20°͒ in x-ray diffraction measurements. Samples for cross-sectional transmission electron microscopy ͑TEM͒ ex- amination were prepared by standard sandwiching, slicing, a͒Presently at: Argonne National Laboratory, Argonne, IL; electronic mail: dimpling, and argon ion milling on a liquid nitrogen-cooled mark–[email protected] b͒Presently at: the U.S. Department of Energy, Germantown, MD. stage at 4 kV and 8° to 11°.
0003-6951/2003/83(19)/3891/3/$20.003891 © 2003 American Institute of Physics Downloaded 04 Nov 2003 to 146.186.179.145. Redistribution subject to AIP license or copyright, see http://ojps.aip.org/aplo/aplcr.jsp 3892 Appl. Phys. Lett., Vol. 83, No. 19, 10 November 2003 Zurbuchen et al.
FIG. 1. Schematic of the planar defect morphology of OPBs in SrRu2O4 /LaAlO3 films; diagonal lines represent OPBs. Two OPBs with opposing inclination, on the far right, have annihilated where they met. On the far left, two OPBs with parallel inclination penetrate through the full thickness of the film.
High-resolution TEM ͑HRTEM͒ imaging was performed using JEOL 4000-EX and Hitachi HF-2000 microscopes. Re- gions of appropriate and consistent thickness of the film– ͑ ͒ substrate phase boundary and an individual OPB were lo- FIG. 2. HRTEM image of a single OPB in a 001 SrRu2O4 /(001)c LaAlO3 film, with a model and simulation of the defect inset. Large circles cated, and through-focus image series were acquired. represent strontium, small circles ruthenium. An extra SrO layer is incorpo- Appropriate multislice simulation parameters were deter- rated along the OPB. The c-direction misregistry across the OPB is 0.25 nm, mined by matching of images to simulated images of the equal to the c-direction dimension of a single layer of SrO. individual phases, and supercells based on possible interfa- cial and defect structures were constructed. HRTEM image found. Close examination of the images reveals that only the simulations used sufficiently large supercells to remove cell- model in Fig. 3͑d͒, with a charge-neutral SrO layer nucleat- edge distortions, caused by their nonrepeating nature, from ing on a LaO-terminated surface, matches the HRTEM im- the region of interest. Bulk values for a and b of Sr2RuO4 age. This result is consistent with the previous determination Ϸ and LaAlO3 were used in constructing the interfacial models, of a LaO surface termination of LaAlO3 above 300 °C in with cells aligned in the center. The 2.1% lattice mismatch vacuum or near vacuum.17,19 Growth proceeds with a 16 ͑ ͒ between (001)c LaAlO3 and 001 Sr2RuO4 at room tem- SrO-RuO2-SrO layer ordering. The Sr2RuO4 unit cell con- perature resulted in negligible variation in alignment ͑0.12 tains two structurally equivalent SrO layers, and only one ͒ Å over the length scale of the supercell. As a check, simu- RuO2 layer, and the OPB offset is equivalent to one SrO lations of a coherent interface were also performed, with layer. There are, then, two possibilities for the nucleation identical results. Only the 2.1% lattice-mismatched supercell mechanism. ͑i͒ If SrO is strongly energetically preferred as simulations are presented here. the starting layers, then a LaO-SrO-SrO-RuO2 growth order, Figure 1 shows a schematic of the generalized defect as shown in Fig. 3͑c͒, might also be possible. However, this morphology we have observed in (001) Sr2RuO4/ layer order contains three AO ͑rocksalt͒-type layers in a row, 16 8,23 (001)c LaAlO3 epitaxial films. Films contain a number which has never been observed in these RP nor in substo- of ⌬cϷ0.25 nm (ϳc/5) OPBs, inclined 17° from the c di- ichiometric perovskite type phases,26 so that it is considered rection. These features either penetrate the full film thick- to be unlikely. ͑ii͒ If the energy difference of a SrO versus a ness, or else annihilate where pairs of opposite inclination RuO2 nucleation layer is small, then a LaO-RuO2-SrO-SrO meet. Similar features in other nϭ1 Ruddlesden-Popper ͑RP͒25 phases have been reported,11,12 although no detailed information on their structure or nucleation is available. A ͓100͔ HRTEM image of a single OPB, showing a ͑011͒ habit, is shown in Fig. 2. A model and HRTEM simu- lation are inset, showing a very good match with the image. Both chemical and physical interruptions of the structure ex- ist at the OPB. A distorted rocksalt SrO layer is incorporated along the OPB. This lattice interruption not only will scatter phonons,8 but also interrupts the continuity of the a – b plane RuO2 layers, believed to be responsible for superconductiv- ity in the material. Figure 3͑a͒ shows a ͓100͔ HRTEM image of the film– substrate interface, with a Fourier-filtered region of the im- age on the right. Bright spots correspond to atomic columns. ͑ ͒ ͑ ͒ ͑ FIG. 3. a HRTEM image of the 001 Sr2RuO4 /(001)c LaAlO3 see Ref. Models and HRTEM simulations for each of the possible 16͒ film–substrate interface, taken at Scherzer defocus. A Fourier-filtered termination-nucleation layer pairs are shown in Figs. 3͑b͒– region of the image is overlaid on the right. ͑b͒–͑f͒ Models of the possible 3͑f͒, with the exception of growth beginning with RuO on termination-nucleation layer pairs for the film, shown with HRTEM simu- 2 lations (⌬ f ϭϪ58 nm,tϭ4.5 nm). The model in ͑d͒ shows the only reason- an AlO -terminated surface, which was ruled out because no 2 able match with the image. The film nucleated with a SrO-RuO2-SrO layer known phases with the corresponding structure could be sequence on a LaO-terminated substrate. Downloaded 04 Nov 2003 to 146.186.179.145. Redistribution subject to AIP license or copyright, see http://ojps.aip.org/aplo/aplcr.jsp Appl. Phys. Lett., Vol. 83, No. 19, 10 November 2003 Zurbuchen et al. 3893
LaO-terminated substrate, and continuing with either se- quence. Perhaps the best strategy would be to use a ͑001͒ 16 LaSrAlO4 substrate instead of (001)c LaAlO3 , as LaSrAlO4 is isostructural to Sr2RuO4. Lastly, extra kinetic activation ͑e.g., higher temperature, higher laser fluence͒ during growth of the first few atomic layers could increase surface diffusivities, increasing the average nuclei 28 FIG. 4. Nucleation-layer growth-order OPB nucleation mechanism in the diameter, thereby reducing the density of OPBs. ͑ ͒ ͑ ͒ 001 Sr2RuO4 /(001)c LaAlO3 see Ref. 16 film. Because two structurally equivalent SrO layers exist in the Sr2RuO4 unit cell, adjacent nuclei growing The authors gratefully acknowledge the financial support with either SrO-RuO2-SrO or RuO2-SrO-SrO ordering, will form an OPB of the Department of Energy through contract DE-FG-02- where they meet. The coordination octahedra of the ruthenium and alumi- num are shown, with circles representing strontium and lanthanum, in 03ER46063. Microscopy was performed in the University of shades of gray equivalent to those in Fig. 3. A unit cell of Sr2RuO4 is shown Michigan’s EMAL, and Penn State’s MCL. to the right for reference. growth order, as shown in Fig. 3͑b͒, might also be possible. 1 Y. Maeno, H. Hashimoto, K. Yoshida, S. Nishizaki, T. Fujita, J. J. Bed- A schematic of this OPB nucleation mechanism, considered norz, and F. Lichtenberg, Nature ͑London͒ 372, 532 ͑1994͒. 2 T. M. Rice and M. Sigrist, J. Phys. C 7, L643 ͑1995͒. most likely to be correct, is shown in Fig. 4. 3 Y. Maeno, T. M. Rice, and M. Sigrist, Phys. Today 54,42͑2001͒. Growth proceeding in either of these scenarios would 4 D. G. Schlom, Y. Jia, L.-N. Zou, J. H. Haeni, S. Briczinski, M. A. Zur- result in the generation of OPBs through the preferred nucle- buchen, C. W. Leitz, S. Madhavan, S. Wozniak, Y. Liu, M. E. Hawley, G. ation mechanism. The energetics of continuing growth after W. Brown, A. Dabkowski, H. A. Dabkowska, R. Uecker, and P. Reiche, Proc. SPIE 3481, 226 ͑1998͒. the first layer would be expected to be quite similar, consid- 5 S. Ohashi, M. Lippmaa, N. Nakagawa, H. Nagasawa, H. Koinuma, and M. ering the close lattice match and the fact that each layer is Kawasaki, Rev. Sci. Instrum. 70,178͑1999͒. charge neutral, so electrostatic effects do not affect growth 6 F. Lichtenberg, A. Catana, J. Mannhart, and D. G. Schlom, Appl. Phys. ͑ ͒ ordering, as they undoubtedly do in other layered perovskites Lett. 60, 1138 1992 . 7 A. P. Mackenzie, R. K. W. Haselwimmer, A. W. Tyler, G. G. Lonzarich, Y. such as the Aurivillius phases. The c-direction misregistry of Mori, S. Nishizaki, and Y. Maeno, Phys. Rev. Lett. 80, 161 ͑1998͒. the OPBs is equal to the c dimension of a single SrO layer 8 M. A. Zurbuchen, Y. Jia, S. K. Knapp, A. H. Carim, D. G. Schlom, L.-N. ͑i.e., an extra one or a missing one͒ within the unit cell. This Zou, and Y. Liu, Appl. Phys. Lett. 78, 2351 ͑2001͒. 9 T. Akima, S. Nishizaki, and Y. Maeno, J. Phys. Soc. Jpn. 68,694͑1999͒. implies that OPBs nucleated at the film–substrate interface 10 ͑ ͒ ͑ ͒ Z. Q. Mao, Y. Mori, and Y. Maeno, Phys. Rev. B 60,610 1999 . almost all of the OPBs observed by the meeting of two 11 J. G. Wen, C. Traeholt, and H. W. Zandbergen, Physica C 205,354͑1993͒. nuclei with differing growth order. 12 G. Kong, M. O. Jones, J. S. Abell, P. P. Edwards, S. T. Lees, K. E. Bearing this model in mind, we reconsider the general Gibbons, I. Gameson, and M. Aindow, J. Mater. Res. 16, 3309 ͑2001͒. 13 M. A. Zurbuchen, J. Lettieri, Y. Jia, D. G. Schlom, S. K. Streiffer, and M. OPB morphology depicted in Fig. 1. Because neighboring E. Hawley, J. Mater. Res. 16, 489 ͑2001͒. out-of-phase nuclei have only one possible offset, OPB in- 14 M. Kawasaki, K. Takahashi, T. Maeda, R. Tsuchiya, M. Shinohara, O. teraction is uncomplicated. Two OPBs with opposing incli- Ishiyama, T. Yonezawa, M. Yoshimoto, and H. Koinuma, Science 266, nations, like those shown on the right in Fig. 1, will annihi- 1540 ͑1994͒. 15 T. Ohnishi, K. Takahashi, M. Nakamura, M. Kawasaki, M. Yoshimoto, late when the film growth front proceeds past the point and H. Koinuma, Appl. Phys. Lett. 74, 2531 ͑1999͒. 16 where they meet. Two nuclei with parallel inclinations, like Cubic indices; LaAlO3 is cubic at the film growth temperature, and trans- those shown on the left, will not interact and will penetrate forms to rhombohedral on cooling. the full film thickness. 17 J.-P. Jacobs, M. A. San Miguel, and L. J. Alvarez, J. Mol. Struct.: THEOCHEM 390, 193 ͑1997͒. Additional mechanisms that may also lead to the genera- 18 D.-W. Kim, D.-H. Kim, B.-S. Kang, T. W. Noh, D. R. Lee, and K.-B. Lee, tion of OPBs in nϭ1 RP phases on perovskite substrates are Appl. Phys. Lett. 74,2176͑1999͒. 23 19 surface steps, and perhaps in the case of films on LaAlO3 P. A. W. van der Heide and J. W. Rabalais, Chem. Phys. Lett. 297,350 ͑1998͒. substrates, twinning of the substrate upon transformation 20 16 S. Bals, G. Rijnders, D. H. A. Blank, and G. Van Tendeloo, Physica C 355, during cooling from the growth temperature. The slight tilt 225 ͑2001͒. between domains could conceivably be accommodated by a 21 J. C. Jiang and X. Q. Pan, J. Appl. Phys. 89, 6365 ͑2001͒. defect having a similar appearance to the OPBs described 22 M. Salluzzo, C. Aruta, I. Maggio-Aprile, Ø. Fischer, S. Gals, and J. Ze- genhagen, Phys. Status Solidi A 186, 339 ͑2001͒. here, although the linear OPB density in the samples we have 23 ͑у ͒8 M. A. Zurbuchen, J. Lettieri, G. Asayama, Y. Jia, A. H. Carim, D. G. examined 1 of every 17 perovskite cells is much greater Schlom, and S. K. Streiffer ͑unpublished͒. 24 than the linear density of domain boundaries in LaAlO3 is B. Aurivillius, Ark. Kemi 1, 463 ͑1950͒; 1,499͑1950͒; 2,519͑1951͒; 5, likely to be. 39 ͑1953͒. 25 ͑ ͒ Several strategies that might suppress the nucleation- S. N. Ruddlesden and P. Popper, Acta Crystallogr. 10, 538 1957 ; 11,54 ͑1958͒. layer OPB generation mechanism in Sr2RuO4 films are pos- 26 C. N. R. Rao and B. Raveau, Transition Metal Oxides: Structure, Proper- sible, and could lead to lower OPB densities and possibly to ties, and Synthesis of Ceramic Oxides, 2nd ed. ͑Wiley, New York, 1998͒, the first superconducting Sr RuO films. An atomic layer-by- p. 61. 2 4 27 layer growth technique such as MBE27 could be used to en- D. G. Schlom, J. H. Haeni, J. Lettieri, C. D. Theis, W. Tian, J. C. Jiang, and X. Q. Pan, Mater. Sci. Eng., B 87, 282 ͑2001͒. gineer the growth order, for example, beginning growth with 28 B. Dam and B. Sta¨uble-Pu¨mpin, J. Mater. Sci.: Mater. Electron. 9,217 ͑ ͒ either the preferred SrO layer, or perhaps a RuO2 layer, on a 1998 .
Downloaded 04 Nov 2003 to 146.186.179.145. Redistribution subject to AIP license or copyright, see http://ojps.aip.org/aplo/aplcr.jsp PHYSICAL REVIEW LETTERS week ending VOLUME 93, NUMBER 14 1 OCTOBER 2004
Enhancement of the Superconducting Transition Temperature of MgB2 by a Strain-Induced Bond-Stretching Mode Softening
A.V. Pogrebnyakov,1,2,3 J. M. Redwing,2,3 S. Raghavan,2,3 V.Vaithyanathan,2,3 D. G. Schlom,2,3 S.Y. Xu,1,3 Qi Li,1,3 D. A. Tenne,1,3 A. Soukiassian,2,3 X. X. Xi,1,2,3 M. D. Johannes,4 D. Kasinathan,4 W. E. Pickett,4 J. S. Wu,5 and J. C. H. Spence5 1Department of Physics, The Pennsylvania State University, University Park, Pennsylvania 16802, USA 2Department of Materials Science and Engineering, The Pennsylvania State University, University Park, Pennsylvania 16802, USA 3Materials Research Institute, The Pennsylvania State University, University Park, Pennsylvania 16802, USA 4Department of Physics, University of California, Davis, California 95616, USA 5Department of Physics and Astronomy, Arizona State University, Tempe, Arizona 85287, USA (Received 5 June 2004; published 30 September 2004)
We report a systematic increase of the superconducting transition temperature Tc with a biaxial tensile strain in MgB2 films to well beyond the bulk value. The tensile strain increases with the MgB2 film thickness, caused primarily by the coalescence of initially nucleated discrete islands (the Volmer- Weber growth mode.) The Tc increase was observed in epitaxial films on SiC and sapphire substrates, although the Tc values were different for the two substrates due to different lattice parameters and thermal expansion coefficients. We identified, by first-principles calculations, the underlying mechanism for the Tc increase to be the softening of the bond-stretching E2g phonon mode, and we confirmed this conclusion by Raman scattering measurements. The result suggests that the E2g phonon softening is a possible avenue to achieve even higher Tc in MgB2-related material systems.
DOI: 10.1103/PhysRevLett.93.147006 PACS numbers: 74.62.–c, 74.70.–b, 74.78.–w
Since the discovery of superconductivity in MgB2 by [13,14]. Therefore, it suffices to consider only the Nagamatsu et al. [1], many techniques, including substi- contribution to Tc. According to the McMillan-Allen- tution, disorder, and pressure, have been used in an at- Dynes analysis, tempt to further increase the transition temperature Tc. ÿf ;^ Tc / !e ; (1) However, the highest Tc in MgB2 has remained at about 40 K. Doping and chemical substitution in MgB2 are where ! is the phonon frequency, f ; ^ 1 difficult [2]. The successful cases such as substitution of = ÿ ^ ,and^ is similar to the Coulomb repulsion Mg by Al [3] and B with C [4] have so far always sup- . The electron-phonon coupling due to the band and pressed Tc. Although atomic disorder induced by proton the bond-stretching mode becomes [14] irradiation enhances the pinning of vortices, it reduces Tc mjDj2 [5]. Subjecting MgB2 to pressure also causes Tc to de- / ; (2) E2g 2 crease [6]. The only higher Tc values reported are from M! 10 the B isotope effect (Tc 40:2K) [7] and by Hur et al. where the band effective mass m is proportional to the when they exposed B crystals to Mg vapor (zero resist- density of states of holes in the band at the Fermi level, ance Tc0 39:8K) [8]. Recently, we have shown that Tc D is the -band deformation potential, and M is the B of MgB2 films on (0001) SiC increases with film thickness mass. A change in Tc can arise from any combination of [9]. There have been speculations that a tensile strain may changes in !, D, m,or^ . be the cause of enhanced Tc in the films [8,9]. In this The epitaxial MgB2 films in this work were deposited Letter, we show unambiguously a systematic increase of by hybrid physical-chemical vapor deposition (HPCVD) Tc with epitaxial tensile strain in MgB2 films on both SiC [15]. Pure magnesium chips were heated simultaneously and sapphire substrates to well beyond the bulk value (the with the substrate to 720 C to generate a high Mg pres- highest Tc0 41:8K). Moreover, we identify the under- sure, and 1000 ppm diborane (B2H6)inH2 was used as the lying mechanism, among several materials’ characteris- B precursor. The carrier gas was a H2 flow of 450 sccm at tics that determine Tc, to be the softening of the E2g a pressure of 100 Torr. The films were deposited on both phonon. (0001) 4H-SiC and (0001) sapphire substrates. Films on MgB2 is a clear and rather extreme example of a ‘‘two (0001) SiC are c-axis oriented and epitaxial with an in- gap’’ superconductor [10–12], arising due to two qualita- plane alignment of the a axis of MgB2 with that of SiC tively different Fermi surfaces (called and ) and their [16]. Films on (0001) sapphire are also c-axis oriented and different pairing strengths (extremely strong and weak, epitaxial, but the a axis of MgB2 is rotated by 30 in respectively). The contribution is dominant, however, plane from the a axis of sapphire [15]. and specifically the contribution from the B-B stretch We have reported previously that the properties of modes (of E2g symmetry) in the B2 graphene layer MgB2 films on (0001) SiC depend on film thickness [9].
147006-1 0031-9007=04=93(14)=147006(4)$22.50 2004 The American Physical Society 147006-1 PHYSICAL REVIEW LETTERS week ending VOLUME 93, NUMBER 14 1 OCTOBER 2004 3.53 (a)
Å) 3.52 3.51 c axis 3.11 3.10 ce Parameter ( Parameter ce 3.09 Latti 3.08 a axis
3.07 1.15 (b)
Ratio 1.14 c/a
1.13 29.4 ) 3
Å 29.2
FIG. 1. Superconducting transition in MgB2 films of different thicknesses on (a) (0001) sapphire and (b) (0001) SiC sub- 29.0 strates. (c) Zero-resistance Tc0 as a function of film thickness Volume ( 28.8 on both sapphire and SiC substrates. (c) 40 41 42 We found a similar trend of thickness dependence on T (K) (0001) sapphire substrates. Figure 1 shows resistivity ver- c0 sus temperature curves for MgB films with different 2 FIG. 2. (a) Tc0 versus a-axis and c-axis lattice constants. thicknesses grown on sapphire [Fig. 1(a)] and SiC (b) Tc0 versus the c=a ratio. (c) Tc0 versus the MgB2 unit cell [Fig. 1(b)] substrates. Both figures clearly show that Tc volume. The open symbols are for sapphire, and the solid becomes higher as the film thickness increases. The thick- symbols are for the SiC substrate. The error bars indicate the ness dependence of the zero-resistance Tc0 for both sub- accuracy of the measurement. The dashed lines are the bulk strates is plotted in Fig. 1(c). A clear trend of increasing values. Tc0 with film thickness is seen for both substrates with Tc0 consistently 1–1.5 K higher in films on SiC than on beyond the bulk value, and the thicker the films the larger sapphire. The highest Tc0 observed (41.8 K) is 2 K higher the tensile strain. than the bulk value. An increase of tensile strain with film thickness has X-ray diffraction analysis, from both ÿ 2 and been observed in films grown in the Volmer-Weber mode: scans, shows a direct correlation between the film strain film growth by initial nucleation of discrete islands which and Tc0. Figure 2(a) shows the a-axis and c-axis lattice later coalesce [17,18]. When the islands coalesce, they constants of several MgB2 films, of different thicknesses ‘‘zip up’’ because the surface energy of the islands is and on sapphire (open symbols) and SiC (solid symbols) larger than the free energy of the grain boundaries, thus substrates, versus their Tc0 values. Clearly, a higher Tc0 creating a tensile strain which increases with film thick- corresponds to a larger a-axis and smaller c-axis lattice ness [19]. In the so-called ‘‘low-mobility’’ films, the constants. Consequently, as Tc0 increases, the c=a ratio height of the zipped boundary is less than the film thick- decreases [see Fig. 2(b)] and the unit cell volume in- ness when the islands coalesce, and the ‘‘zipping up’’ creases [see Fig. 2(c)]. Most strikingly, although the re- continues at the surface terrace level long afterwards sults on SiC and sapphire follow two different Tc0 versus [20]. Therefore, the tensile strain continues to increase thickness curves in Fig. 1(c), they fall on the same curves with film thickness toward the upper bound value pre- in Fig. 2. Compared to the bulk values indicated by the dicted by the Nix-Clemens model [19]. dashed lines in Fig. 2 (a 3:086 A , c 3:524 A [1], Atomic force microscopy (AFM) images confirm that 3 c=a 1:142, and cell volume of 29:06 A ), it is evident MgB2 films grow in the Volmer-Weber mode. AFM im- that the tensile strain in the films causes Tc0 to increase ages of two MgB2 films on SiC are shown in Figs. 3(a)
147006-2 147006-2 PHYSICAL REVIEW LETTERS week ending VOLUME 93, NUMBER 14 1 OCTOBER 2004
(thickness 75 A ) and 3(b) (thickness 900 A ). Hexagonal- combined effect of all the sources of strain results in shaped MgB2 crystallites are seen in the thinner film, larger tensile strain in MgB2 films on SiC than on which coalesce into a continuous film at larger film thick- sapphire. It is difficult to accurately predict the ness. The continued increase of the tensile strain beyond maximum tensile strain attainable in MgB2 films. A coalescence and the surface terraces that are readily ob- rough estimate from Eq. (1) in Ref. [19] [taking the servable in thicker films indicate that MgB2 behaves like a following parameters: grain boundary energy low-mobility system at 720 C. This is consistent with the 0J=cm2, solid/vapor surface energy 1J=m2, biaxial high melting temperature of B (2075 C)orMgB2 (it modulus of MgB2 421 GPa [26], and grain size at melts congruently at 2430 C with pressure higher than coalescence 1300 A from Fig. 3(a)] indicates an upper 6.5 MPa [21]). The grain growth as the film becomes bound value of 0:8%, which is larger than the maxi- thicker may also contribute to the tensile strain [22]. mum tensile strain reported here ( 0:55%). However, This is supported by the cross-section transmission elec- this stress estimation may be oversimplified and the films tron microscopy (TEM) images of a 2100 A thick MgB2 may partially relax as the thickness increases. film on (0001) SiC in Fig. 3(c) (low magnification) and The systematic increase of Tc with tensile strain shown Fig. 3(d) (high resolution), which show that the defect in Fig. 2 has led us to search for the mechanism respon- density is high near the film/substrate interface and de- sible for the rise in Tc of MgB2 by first-principles calcu- creases as the film thickness increases. lations. For the highest Tc sample in Fig. 2, The lattice mismatch could result in 11% compressive a=a 0:55% and c=c ÿ0:25%, which leads to strain (with a 30 in-plane rotation) in MgB2 films on changes m =m ÿ1:2% and 2jDj=jDj ÿ2:6%. sapphire (a 3:07 A ) and 0.1% compressive strain on The calculated change in the E2g phonon frequency, SiC (a 4:765 A ). The coefficient of thermal expansion !=! ÿ5:5%, is in agreement with the experimen- in the (0001) plane at room temperature is 5:5 tally extrapolated value [27]. From these values and using ÿ6 ÿ1 ÿ6 ÿ1 Eq. (2), we obtain 10 K for MgB2 [23], 3:0 10 K for SiC [24], ÿ6 ÿ1 and 6:7 10 K for sapphire [25]. Therefore, the E2g m jDj ! thermal expansion mismatch will lead to a tensile strain 2 ÿ 2 7%: (3) m j j ! in films on SiC, and a compressive strain on sapphire. The E2g D
From Eq. (1) and using representative numbers, 0:9 0:2 1:1,and^ 0:2, T ! c ÿ f ÿ5:5% 10:4% 5%: (4) Tc ! This value is quite close to the experimental enhancement of Tc. This analysis identifies the underlying mechanism of the tensile-strain-induced increase in Tc as the de- crease in the E2g phonon frequency. It leads to the large (7%) increase in the band coupling strength, which, being leveraged in f ; ^ , more than compensates the 5% lowering of the temperature/energy scale governed by !. The result is consistent with that of Yildirim and Gu¨lseren on the pressure effects in MgB2 [28], although this previous work did not address biaxial strain. The E2g phonon frequencies of MgB2 films of different thicknesses were measured by Raman scattering using a SPEX Triplemate spectrometer in backscattering geome- try with the 514.5 nm Ar laser excitation. Figure 4 shows the spectra of two films: a 3400 A thick film on SiC and a 800 A thick film on sapphire. The thicker film on SiC has larger tensile strain than the thinner film on FIG. 3. AFM images of (a) a 75 A and (b) a 900 A thick sapphire and consequently has a higher Tc0 as compared MgB2 films grown on (0001) SiC substrates at 720 C. (c) Low- to the film on sapphire. The broad peak around 600 cmÿ1 magnification and (d) high-resolution cross-section TEM im- ages of a 2100 A thick MgB film on (0001) SiC taken along the has been assigned to the E2g mode [27]. A clear difference 2 between the peak energies for the two films is seen. The 1100 direction. The dotted line in (d) indicates the MgB2=SiC interface, and the small arrows show the lattice deformation higher Tc film on SiC has a softer E2g mode than the ÿ1 due to dislocations. lower Tc film on sapphire by about 20 cm , i.e., 3.3%.
147006-3 147006-3 PHYSICAL REVIEW LETTERS week ending VOLUME 93, NUMBER 14 1 OCTOBER 2004
on SiC, 3400Å DOE under Grants No. DE-FG02-03ER46063 (G. S.) on Al O , 800Å and No. DE-FG02-01ER45907 (X. X. X.). 200 2 3
[1] J. Nagamatsu et al., Nature (London) 410, 63 (2001). 100 [2] R. J. Cava, H.W. Zandbergen, and K. Inumaru, Physica (Amsterdam) 385C, 8 (2003).
Intensity (arb.units) [3] J. S. Slusky et al., Nature (London) 410, 343 (2001). [4] S. Lee et al., Physica (Amsterdam) 397C, 7 (2003). [5] Y. Bugoslavsky et al., Nature (London) 410, 563 (2001). 0 [6] M. Monteverde et al., Science 292, 75 (2001). 200 300 400 500 600 700 800 -1 [7] D. G. Hinks, H. Claus, and J. D. Jorgensen, Nature Raman shift (cm ) (London) 411, 457 (2001). [8] N. Hur et al., Appl. Phys. Lett. 79, 4180 (2001). FIG. 4. Raman spectra of two MgB2 films: 3400 A thick on [9] A.V. Pogrebnyakov et al., Appl. Phys. Lett. 82,4319 SiC and 800 A thick on sapphire. The former shows a lower E2g (2003). phonon frequency. [10] H. J. Choi et al., Nature (London) 418, 758 (2002). [11] W. Pickett, Nature (London) 418, 733 (2002). [12] S. Souma et al., Nature (London) 423, 65 (2003). In conclusion, we have shown unambiguously that [13] A.Y. Liu, I. I. Mazin, and J. Kortus, Phys. Rev. Lett. 87, 087005 (2001). Tc0 in MgB2 films deposited by HPCVD increases [14] J. M. An and W. E. Pickett, Phys. Rev. Lett. 86, 4366 with biaxial tensile strain. The highest Tc0 obtained (2001). in the MgB2 films, 41.8 K, is well above the bulk [15] X. H. Zeng et al., Nat. Mater. 1, 35 (2002). value. First-principles calculations show that the Tc in- [16] X. H. Zeng et al., Appl. Phys. Lett. 82, 2097 (2003). crease is due to the softening of the bond-stretching E2g [17] J. A. Floro et al., MRS Bull. 27, No. 1, 19 (2002). phonon mode, which is confirmed by Raman scattering [18] R. Koch, J. Phys. Condens. Matter 6, 9519 (1994). [19] W. D. Nix and B. M. Clemens, J. Mater. Res. 14, 3467 measurements. There have been reports that a reduced Tc (1999). in MgB2 corresponds to a higher E2g phonon frequency [27,29]. Our result is the first example that lowering the [20] B.W. Sheldon, K. H. A. Lau, and A. Rajamani, J. Appl. Phys. 90, 5097 (2001). E2g phonon frequency can increase Tc to above the bulk [21] Z. K. Liu, D. G. Schlom, Q. Li, and X. X. Xi, Appl. Phys. value. In this work, the Tc enhancement seems to saturate Lett. 78, 3678 (2001b). at large film thicknesses. If one could reduce the E2g [22] P. Chaudhari, J. Vac. Sci. Technol. 9, 520 (1972). phonon frequency further using other techniques, the [23] C. Buzea and T. Yamashita, Supercond. Sci. Technol. 14, Tc of MgB2 should be even higher than what we have R115 (2001). shown here. [24] Z. Li and R. Brandt, J. Am. Chem. Soc. 69, 863 (1986). We thank Bill Nix for useful discussions on the coales- [25] Y.S. Touloukian et al., Thermophysical Properties of cence tensile strain and E. Wertz for his help with the Matter (Plenum, New York, 1977), Vol. 13. temperature calibration. The work is supported in part by [26] V. Milman and M. C. Warren, J. Phys. Condens. Matter 13, 5585 (2001). ONR under Grants No. N00014-00-1-0294 (X. X. X.) and [27] A. F. Goncharov and V.V. Struzhkin, Physica No. N0014-01-1-0006 (J. M. R.), by NSF under Grants No. (Amsterdam) 385C, 117 (2003). DMR-0306746 (X. X. X. and J. M. R.), No. DMR- [28] T. Yildirim and O. Gu¨lseren, J. Phys. Chem. Solids 63, 9876266 and No. DMR-9972973 (Q. L.), No. DMR- 2201 (2002). 0103354 (G. S. and X. X. X.), No. DMR-0114818 [29] T. Masui, S. Lee, and S. Tajima, Phys. Rev. B 70, 024504 (W. E. P.), and No. DMR-0245702 (J. C. H. S.), and by (2004).
147006-4 147006-4 APPLIED PHYSICS LETTERS VOLUME 85, NUMBER 11 13 SEPTEMBER 2004
Properties of MgB2 thin films with carbon doping A. V. Pogrebnyakova) and X. X. Xi Department of Physics, Department of Materials Science and Engineering and Materials Research Institute, The Pennsylvania State University, University Park, Pennsylvania 16802 J. M. Redwing, V. Vaithyanathan, D. G. Schlom, and A. Soukiassian Department of Materials Science and Engineering and Materials Research Institute, The Pennsylvania State University, University Park, Pennsylvania 16802 S. B. Mi and C. L. Jia Institut für Festkörperforschung, Forschungszentrum Jülich GmbH, D-52425 Jülich, Germany J. E. Giencke and C. B. Eom Department of Materials Science and Engineering and Applied Superconductivity Center, University of Wisconsin, Madison, Wisconsin 53706 J. Chen, Y. F. Hu, Y. Cui, and Qi Li Department of Physics and Materials Research Institute, The Pennsylvania State University, University Park, Pennsylvania 16802 (Received 29 March 2004; accepted 15 June 2004)
We have studied structural and superconducting properties of MgB2 thin films doped with carbon during the hybrid physical-chemical vapor deposition process. A carbon-containing precursor metalorganic bis(methylcyclopentadienyl)magnesium was added to the carrier gas to achieve carbon doping. As the amount of carbon in the film increases, the resistivity increases, Tc decreases, and the upper critical field increases dramatically as compared to clean films. The self-field Jc in the carbon doped film is lower than that in the clean film, but Jc remains relatively high to much higher magnetic fields, indicating stronger pinning. Structurally, the doped films are textured with columnar nano-grains and highly resistive amorphous areas at the grain boundaries. The carbon doping approach can be used to produce MgB2 materials for high magnetic-field applications. © 2004 American Institute of Physics. [DOI: 10.1063/1.1782258] 1 ͑ ͒ The 39 K superconductor MgB2 has attracted tremen- and 21.6 °C to transport MeCp 2Mg to the reactor. The dous interest for its potential in high magnetic-field applica- amount of carbon doping depends on the secondary hydro- ͑ ͒ tions. In particular, high critical current density Jc near the gen flow rate through the MeCp 2Mg bubbler, which was 2 depairing limit has been observed in MgB2, and unlike high varied between 25 and 200 sccm to vary the flow rate of ͑ ͒ temperature superconductors grain boundaries in MgB2 do MeCp 2Mg into the reactor from 0.0065 to 0.052 sccm. 3 not behave like weak links. Although clean MgB2 shows The chemical compositions of a series of films were deter- 4 low upper critical field Hc2, in high resistivity MgB2 films mined by wavelength dispersive x-ray spectroscopy to estab- 5 Hc2 is substantially higher; the critical current density Jc in lish a correlation between the carbon concentrations in the 6 ͑ ͒ clean MgB2 is suppressed quickly by magnetic field, but films and the hydrogen flow rates through the MeCp 2Mg defects and impurities have been shown to enhance bubbler. The nominal atomic concentrations determined by pinning.7,8 Recently, we have developed a hybrid physical- this approach are used as the carbon concentrations presented chemical vapor deposition (HPCVD) technique which pro- in this letter. The films were deposited on (001) 4H–SiC 10 duces in situ epitaxial MgB2 films with Tc above 40 K. substrates at 720 °C. The thickness of the films was around Because of the highly reducing H2 ambient and the high 2000 Å. purity B2H6 source used in the HPCVD process, the tech- The resistivity (in log scale) versus temperature curves for MgB films with different carbon doping levels are nique produces very clean MgB2 thin films. In this letter, we 2 shown in Fig. 1 a . The carbon doping causes a dramatic describe the properties of HPCVD MgB2 films doped with ( ) increase in the resistivity, whereas the T of the film is sup- carbon. The Hc2 values in these “dirtier” films are signifi- c cantly higher and the vortex pinning is significantly stronger pressed much more slowly. For example, with a carbon con- than those in the clean films. centration of 24 at. %, the residual resistivity increases from the undoped value of less than 1 ⍀ cm to ϳ200 ⍀ cm, In situ epitaxial growth of MgB2 films by HPCVD has been described in detail previously.9 For carbon but Tc only decreases from over 41 to 35 K. The dependen- cies of residual resistivity and T on the carbon concentration doping, we added bis(methylcyclopentadienyl)magnesium c in the doped MgB films are plotted in Fig. 1(b). T is sup- [͑MeCp͒ Mg], a metalorganic magnesium precursor, to the 2 c 2 pressed to below 4.2 K at a nominal carbon concentration of carrier gas of 300 sccm hydrogen at 100 Torr. The flow of the 42 at. % when the residual resistivity is 440 m⍀ cm. This is boron precursor gas, 1000 ppm diborane ͑B H ͒ in H , was 2 6 2 very different from those in carbon-doped single crystals, kept at 150 sccm. A secondary hydrogen flow was passed where T is suppressed to 2.5 K at a residual resistivity of ͑ ͒ c through the MeCp 2Mg bubbler which was held at 760 Torr ⍀ 11 50 cm when 12.5 at. % of carbon is doped into MgB2. This discrepancy indicates that only a small portion of the a) Electronic mail: [email protected] carbon in the films is doped into the MgB2 structure and the
0003-6951/2004/85(11)/2017/3/$22.002017 © 2004 American Institute of Physics Downloaded 01 Nov 2005 to 128.118.88.140. Redistribution subject to AIP license or copyright, see http://apl.aip.org/apl/copyright.jsp 2018 Appl. Phys. Lett., Vol. 85, No. 11, 13 September 2004 Pogrebnyakov et al.
FIG. 1. (a) Resistivity vs temperature curves for MgB2 films of different carbon doping. (b) Residual resistivity (closed circles) and T (open circles) c FIG. 3. a X-ray diffraction –2 scans for MgB films with carbon dop- as a function of carbon concentration for films plotted in (a).In(a), from ( ) 2 bottom to top, the nominal carbon concentrations of the curves are 0, 7.4, ing. From top to bottom, the nominal carbon concentrations are 0, 7.4, 15, 15, 22, 29, 34, 39, 42, and 45 at. %. 22, 28, 39, 42, and 45 at. %. The spectra are shifted vertically for clarity. The peaks labeled with an asterisk are due to the SiC substrate peaks. (c) The c-axis lattice constant (open triangles) and a-axis lattice constant rest most likely forms high resistance grain boundaries giv- (closed squares) of the carbon doped MgB2 films as a function of nominal ͑ ͒ 12 carbon concentration. ing rise to poor connectivity of the Mg B1−xCx 2 grains. The granular structure of the carbon-doped MgB2 films is confirmed by TEM. Figure 2(a) is a cross-sectional TEM In the planar-view image in Fig. 2(c), the change of contrast image of a film with 22 at. % nominal carbon concentration indicates an equiaxial in-plane morphology of the columnar taken along the ͓¯110͔ direction of the substrate. It shows that grains, and an amorphous phase is also observed between the the film consists of columnar nano-grains (the contrast grains. We were not able to determine the composition of the ͑ ͒ amorphous areas, but it is most likely that the large portion changes laterally, but not vertically) of Mg B1−xCx 2 with a preferential c-axis orientation. The selected area electron dif- of carbon that is not doped into MgB2 is contained in these areas. The insert in Fig. 2(c) is a typical diffraction pattern fraction pattern taken from the MgB2 /SiC interface area in Fig. 2(b) shows two types of features, diffraction spots and taken along the film normal. The strong hexagonal- arcs. The spots belong to the single crystal SiC substrate distributed spots show that the hexagonal-on-hexagonal in- plane relationship between the columnar grains and SiC (SC) and the arcs to the MgB2 film (MB). The arcs consist of many fine spots originating from individual columnar grains dominates, while the diffraction rings reveal grains that are which show a deviation of their c axis from the film normal. randomly in-plane oriented. Figure 3(a) shows –2 scans of an undoped MgB2 film and films doped with different amounts of carbon. Compared to the undoped films, the MgB2 00l peaks are suppressed as carbon concentration increases, and dramatically when the carbon concentration is above ϳ30 at. %. Meanwhile, as shown in Fig. 3(b), both the c and a axes expand until about 30 at. %, above which the c lattice constant decreases and the a lattice constant increases dramatically. This behavior is different from that in carbon-doped single crystals, where the a axis lattice constant decreases but that of c axis remains almost constant for all the carbon concentration.11 The peak marked by “?” is likely 101 MgB2, the most intense diffrac- tion peak of MgB2. It becomes stronger as the carbon con- centration increases, indicating an increased presence of ran- domly oriented MgB2. The peaks marked by “boron carbide,” according to extensive pole figure analysis, are most likely due to B4C, B8C, or B13C2. Their intensities also FIG. 2. TEM results of a film with 22 at. % nominal carbon concentration. increase with carbon concentration. From the TEM and x-ray (a) Cross-sectional image taken along the ͓¯110͔ direction of the substrate. diffraction results, we conclude that below about 30 at. %, a (b) Selected area electron diffraction taken from the MgB2 /SiC interface ͑ ͒ small portion of carbon is doped into the Mg B1−xCx 2 co- area. (c) Planar-view image showing nano-grains of carbon doped MgB2 and an amorphous phase between the grains. The insert is the select area lumnar, c-axis-oriented nano-grains, and the rest goes into electron diffraction pattern taken along the film normal. the grain boundaries consisting of highly resistive amor- Downloaded 01 Nov 2005 to 128.118.88.140. Redistribution subject to AIP license or copyright, see http://apl.aip.org/apl/copyright.jsp Appl. Phys. Lett., Vol. 85, No. 11, 13 September 2004 Pogrebnyakov et al. 2019
While the undoped film has high self-field critical current densities, they are suppressed quickly by magnetic field due to the weak pinning. For the film doped with 11 at. % nomi- nal carbon concentration, Jc values are relatively high in much higher magnetic fields. This indicates a significantly enhanced vortex pinning in carbon doped MgB2 films. In conclusion, MgB2 thin films were doped with carbon by adding bis(cyclopentadienyl)magnesium to the carrier gas during the HPCVD process. The residual resistivity increases rapidly while Tc decreases much more slowly with carbon doping. Structural analyses show that only part of the carbon is doped into the MgB2 lattice and the rest forms highly FIG. 4. Upper critical field as a function of temperature for an undoped film ͑ ͒ and two doped films with 7.4 at. % and 22 at. % nominal carbon concentra- resistive foreign phases in the grain boundaries. Hc2 T and ͑ ʈ ͒ tions, respectively. The closed symbols are for parallel field Hc2 and the its slope near Tc increase dramatically as compared to the ͑ Ќ ͒ open symbols are for perpendicular field Hc2 . clean films, which is consistent with the multiband supercon- ductor model of Gurevich and indicates a dirtier band 14 phous phases or boron carbides. Above about 30 at. %, the upon carbon doping. The critical current density in mag- ͑ ͒ Mg B1−xCx 2 nano-grains are completely separated from netic field increases markedly from that in the clean film due each other by highly resistive phases, become more ran- to stronger vortex pinning. domly oriented, and their lattice constants relax. This is con- sistent with the result in Fig. 1(b). The work is supported in part by ONR under grant Nos. The upper critical field H was measured using a Quan- c2 N00014-00-1-0294 (X.X.X.) and N0014-01-1-0006 (J.M.R.), tum Design PPMS system with a 9 T superconducting mag- by NSF under grant Nos. DMR-0306746 (X.X.X. and net. Figure 4 shows the results for an undoped, 7.4 at. %, and J.M.R.), DMR-9876266 and DMR-9972973 (Q.L.), and 22 at. % carbon doped films. The value of H is defined by c2 through the MRSEC for Nanostructure Materials (C.B.E.), 50% of the normal-state resistance R͑H ͒=0.5R͑T ͒.Itcan c2 c and by DOE under grant No. DE-FG02-03ER46063 be clearly seen that carbon doping changes the downward Ќ ͑ ͒ (D.G.S.). curvature in Hc2 T for the undoped film to an upward cur- 1J. Nagamatsu, N. Nakagawa, T. Muranaka, Y. Zenitani, and J. Akimitsu, vature in the carbon doped films. Both the slope, dHc2 /dT, Nature (London) 410,63(2001). near Tc and the low temperature Hc2 increase with carbon 2S. Y. Xu, Q. Li, E. Wertz, Y. F. Hu, A. V. Pogrebnyakov, X. H. Zeng, X. concentration. In high magnetic field measurements, Brac- X. Xi, and J. M. Redwing, Phys. Rev. B 68, 224501 (2003). 3 cini et al. have shown that carbon-doped MgB2 films as de- D. C. Larbalestier, L. D. Cooley, M. O. Rikel, A. A. Polyanskii, J. Jiang, ͑ ͒ scribed here have extraordinary Hc2 0 values as high as S. Patnaik, X. Y. Cai, D. M. Feldmann, A. Gurevich, A. A. Squitieri, M. T. 13 ͑ ͒ Naus, C. B. Eom, E. E. Hellstrom, R. J. Cava, K. A. Regan, N. Rogado, 70 T. The transport Jc H at different temperatures, deter- M. A. Hayward, T. He, J. S. Slusky, P. Khalifah, K. Inumaru, and M. mined by a 1 V criterion from 20–50 m bridges, for an Haas, Nature (London) 410, 186 (2001). 4 undoped and a carbon doped MgB2 film are shown in Fig. 5. P. C. Canfield and G. Crabtree, Phys. Today 56,34(2003). 5A. Gurevich, S. Patnaik, V. Braccini, K. H. Kim, C. Mielke, X. Song, L. D. Cooley, S. D. Bu, D. M. Kim, J. H. Choi, L. J. Belenky, J. Giencke, M. K. Lee, W. Tian, X. Pan, A. Siri, E. E. Hellstrom, C. B. Eom, and D. Larbalestier, Supercond. Sci. Technol. 17, 278 (2004). 6X. H. Zeng, A. V. Pogrebnyakov, M. H. Zhu, J. E. Jones, X. X. Xi, S. Y. Xu, E. Wertz, Q. Li, J. M. Redwing, J. Lettieri, V. Vaithyanathan, D. G. Schlom, Z. K. Liu, O. Trithaveesak, and J. Schubert, Appl. Phys. Lett. 82, 2097 (2003). 7Y. Bugoslavsky, G. K. Perkins, X. Qi, L. F. Cohen, and A. D. Caplin, Nature (London) 410, 563 (2001). 8C. B. Eom, M. K. Lee, J. H. Choi, L. Belenky, X. Song, L. D. Cooley, M. T. Naus, S. Patnaik, J. Jiang, M. O. Rikel, A. A. Polyanskii, A. Gurevich, X. Y. Cai, S. D. Bu, S. E. Babcock, E. E. Héllstrom, D. C. Larbalestier, N. Rogado, K. A. Regan, M. A. Hayward, T. He, J. S. Slusky, K. Inumaru, M. Haas, and R. J. Cava, Nature (London) 411, 558 (2001). 9X. H. Zeng, A. V. Pogrebnyakov, A. Kotcharov, J. E. Jones, X. X. Xi, E. M. Lysczek, J. M. Redwing, S. Y. Xu, Q. Li, J. Lettieri, D. G. Schlom, W. Tian, X. Q. Pan, and Z. K. Liu, Nat. Mater. 1,35(2002). 10A. V. Pogrebnyakov, J. M. Redwing, J. E. Jones, X. X. Xi, S. Y. Xu, Q. Li, V. Vaithyanathan, and D. G. Schlom, Appl. Phys. Lett. 82, 4319 (2003). 11S. Lee, T. Masui, A. Yamamoto, H. Uchiyama, and S. Tajima, Physica C 397,7(2003). 12J. Rowell, Supercond. Sci. Technol. 16, R17 (2003). 13V. Braccini, A. Gurevich, J. Giencke, M. Jewell, C. Eom, D. Larbalestier, A. Pogrebnyakov, Y. Cui, B. T. Liu, Y. F. Hu, J. M. Redwing, Q. Li, X. X. Xi, R. Singh, R. Gandikota, J. Kim, B. Wilkens, N. Newmann, J. Rowell, FIG. 5. Critical current density as a function of magnetic field ͑Hʈc͒ and B. Moeckly, V. Ferrando, C. Tarantini, D. Marr, M. Putti, C. Ferdeghini, temperature for (a) an undoped film and (b) a film doped with with 11 at. % R. Vaglio, and E. Haanappel (unpublished). nominal carbon concentration. 14A. Gurevich, Phys. Rev. B 67, 184515 (2003).
Downloaded 01 Nov 2005 to 128.118.88.140. Redistribution subject to AIP license or copyright, see http://apl.aip.org/apl/copyright.jsp PHYSICAL REVIEW B 69, 174101 ͑2004͒
Absence of low-temperature phase transitions in epitaxial BaTiO3 thin films
D. A. Tenne,1,* X. X. Xi,1,2,3 Y. L. Li,2 L. Q. Chen,2,3 A. Soukiassian,2 M. H. Zhu,1,† A. R. James,1,‡ J. Lettieri,2 D. G. Schlom,2,3 W. Tian,4 and X. Q. Pan4 1Department of Physics, the Pennsylvania State University, University Park, Pennsylvania 16802, USA 2Department of Materials Science and Engineering, the Pennsylvania State University, University Park, Pennsylvania 16802, USA 3Materials Research Institute, the Pennsylvania State University, University Park, Pennsylvania 16802, USA 4Department of Materials Science and Engineering, the University of Michigan, Ann Arbor, Michigan 48109, USA ͑Received 14 January 2004; published 3 May 2004͒
We have studied phase transitions in epitaxial BaTiO3 thin films by Raman spectroscopy. The films are found to remain in a single ferroelectric phase over the temperature range from 5 to 325 K. The low- ͑ ͒ temperature phase transitions characteristic of bulk BaTiO3 tetragonal-orthorhombic-rhombohedral are absent in the films. X-ray diffraction shows that the BaTiO3 films are under tensile strain due to the thermal expansion mismatch with the buffer layer. A phase-field calculation of the phase diagram and domain structures in
BaTiO3 thin films predicts, without any priori assumption, that an orthorhombic phase with in-plane polariza- tion is the thermodynamically stable phase for such values of tensile strain and temperature, consistent with the experimental Raman results.
DOI: 10.1103/PhysRevB.69.174101 PACS number͑s͒: 77.84.Dy, 77.55.ϩf, 78.30.Hv, 63.70.ϩh
I. INTRODUCTION phase transition from cubic to tetragonal phase above room temperature. In the present paper we focus on the tempera- Phase transitions in thin film ferroelectrics, which differ ture range 5–325 K. substantially from those in corresponding bulk materials, are of significant scientific and technological interest.1,2 Much effort has been devoted in recent years to the thermodynam- II. EXPERIMENTAL RESULTS ics of domain structures and phase diagrams in ferroelectric The BaTiO films studied in this work were grown by thin films.3–6 Different approaches have been employed by 3 5 pulsed laser deposition on (001)SrTiO3 and (001)LaAlO3 various groups. For example, Pertsev et al. predicted a num- substrates. Between the substrate and BaTiO3 film a 0.3 m ber of domain stability maps for BaTiO3 and PbTiO3 using 6 thick conducting layer of SrRuO3 was deposited to screen thermodynamic calculations whereas Li et al. proposed a Raman signal from the substrate.24 The thickness of the phase-field model to predict the dependence of domain struc- BaTiO films was 1 m. The deposition conditions of the tures on temperature and strain. Experimentally, enhance- 3 BaTiO3 and SrRuO3 layers were the same as those described ment of the paraelectric-ferroelectric phase transition tem- previously.25 Both layers were deposited at the substrate tem- perature, Tc , by strain has been reported in epitaxial perature of 750 °C. The structural properties of the ferroelectric films.7–10 Recently, Streiffer et al. observed BaTiO3 /SrRuO3 samples were characterized by x-ray dif- nanoscale 180° stripe domains in ferroelectric PbTiO3 thin ͑ ͒ 11 fraction and transmission electron microscopy TEM . films. Strain-induced ferroelectricity in SrTiO3 thin films ͑ 12 13 BaTiO3 single crystal was obtained from MTI Corp. Rich- has been reported by Fuchs et al. and Tikhomirov et al. mond, CA͒. Raman spectra were recorded using a SPEX In this paper, we present a Raman scattering study of BaTiO3 Triplemate spectrometer equipped with a liquid-nitrogen- films in the temperature range 5–320 K, which shows that cooled multichannel coupled-charge-device detector. The the tetragonal-orthorhombic and orthorhombic-rhombohedral 514.5 nm line of an Arϩ laser was used for excitation, and phase transitions characteristic for the bulk BaTiO3 disap- the laser power density was kept at a low level pear completely in the thin films. Detailed structural charac- (Շ30 W/cm2) to avoid heating of the sample. Raman spec- terizations of the films suggest that a tensile strain locks the tra from the films were measured in backscattering geometry films in a single ferroelectric phase. A comparison with the both in parallel z(x,x)¯z and perpendicular z(x,y)¯z polariza- phase-field model calculation confirms that the orthorhombic ͒ phase with in-plane polarization is the stable phase in this tion configurations (z direction is normal to the film plane . strain and temperature range. The measurements in the crystal were in backscattering ge- ometries both along the c axis of the initial tetragonal phase In bulk BaTiO3, the crystal structure changes from cubic to tetragonal at 403 K, then to orthorhombic at 278 K, and to ͓z(x,x)¯z and z(x,y)¯z geometries͔ and perpendicular to it rhombohedral at 183 K.14 The cubic phase is paraelectric and ͓x(z,z)¯x and x(z,y)¯x geometries͔. ͑ ͒ the other phases are ferroelectric. Raman spectroscopy is a Figure 1 shows Raman spectra of a a BaTiO3 single ͑ ͒ powerful tool to study phase transitions as the structural crystal and b a BaTiO3 film grown on a SrTiO3 substrate, changes alter the vibrational spectra. Lattice dynamics and both measured in the parallel polarization configuration, as a Raman spectra of single crystal and polycrystalline BaTiO3 function of temperature. The spectra of the single crystal have been studied in detail.15–22 Raman studies of thin show clear changes attributed to the phase transitions. At low 23 BaTiO3 films have mostly focused on the ferroelectric temperatures, in the rhombohedral phase, the most intensive
0163-1829/2004/69͑17͒/174101͑5͒/$22.5069 174101-1 ©2004 The American Physical Society D. A. TENNE et al. PHYSICAL REVIEW B 69, 174101 ͑2004͒
FIG. 1. Temperature evolution of Raman spectra of ͑a͒ a BaTiO single crystal and ͑b͒ a BaTiO film grown on SrTiO sub- 3 3 3 FIG. 2. Temperature dependence of the A1 (TO2) and A1 (TO3) strate with a SrRuO buffer layer, measured in parallel polarization ͑ ͒ 3 phonon frequencies for BaTiO3 single crystal solid symbols and ¯ ͑ geometry z(x,x)z. Dashed-dotted lines are guides to eye. Star in the film grown on SrTiO3 substrate with SrRuO3 buffer layer open spectra of the film indicates features due to SrRuO3 layer symbols͒.
16–18 lines observed in the polarized Raman spectrum are at 173, BaTiO3. The heavy overdamping of the soft mode in Ϫ1 187, 242, 485, 522, and 714 cm , and are attributed to BaTiO3 has been attributed to the order-disorder character of 26 TO1 ,LO1 ,TO2 ,LO2 ,TO3 ,LO3 modes of A1 symmetry, the phase transitions. In the crystal overdamping is ob- Ϫ respectively.20 The peak at 310 cm 1 is due to mixed served in the orthorhombic and tetragonal phases and in- 20 LO2 –TO3 phonon of E symmetry. The shoulder at creases with temperature, in agreement with earlier reported 261 cmϪ1 was shown by micro-Raman experiments to be results.17,18 In thin films overdamping has also been observed 22 due to the near-domain-wall regions of BaTiO3 crystal. to increase with temperature. However, no qualitative differ- ϳ When the BaTiO3 crystal was warmed to 185 K, the Ra- ence has been observed in the low-frequency Raman spectra man spectra went through large changes, most noticeable in of the films, as well as the crystals in the orthorhombic and the frequency range 150–300 cmϪ1. Sudden jumps to a tetragonal phases; there are no sharp changes in the over- Ϫ1 higher frequency, as large as 35 cm took place for the A1 damping behavior which would indicate the O-T phase tran- (TO2) phonon line. A less dramatic, but still well pro- sitions. Therefore, the behavior of the overdamped soft mode ϳ Ϫ1 nounced jump ( 10 cm ) occurred for the A1 (TO3) line. cannot be used to determine the phase transition temperature. These changes are indicative of the rhombohedral- So we focus on two hard phonon modes, TO2 and TO3, orthorhombic phase transition. When the crystal was heated which exhibit most noticeable changes upon the ϳ further to 280 K, another sudden jumps of the A1 TO2 and rhombohedral-orthorhombic and orthorhombic-tetragonal ͑ Ϫ1 ͒ TO3 phonon frequencies 25 and 4 cm , respectively and phase transitions in BaTiO3 crystals. changes in the Raman spectra occurred, which indicates the The temperature dependence of the frequencies of the A1 orthorhombic-tetragonal phase transition. The observed fre- (TO2) and A1 (TO3) phonon modes for both the crystal and quency shifts are well detectable in Raman spectra despite film shown in Fig. 2 clearly demonstrates the difference be- large linewidths TO2 and TO3 phonons, and are in agreement tween the crystal and film behavior. Large jumps in the A1 15,21,22 with the reported results for BaTiO3 crystals. (TO2) phonon frequency are clearly seen in the crystal, but it The temperature dependence of the Raman spectra of the is constant in the film. The A1 (TO3) phonon frequency ͑ ͒ BaTiO3 film, shown in Fig. 1 b , is markedly different from shows smaller jumps in the crystal, whereas no such jumps that of the single crystal. It does not show sharp changes in are visible in the film. This behavior indicates that the the temperature range of 5–300 K. When temperature in- BaTiO3 film does not undergo any phase transition in the creases, the phonon lines broaden, but their frequencies are temperature range 5–325 K. The same behavior was ob- either almost constant (TO2) or change only very slightly served for films grown both on SrTiO3 and LaAlO3 sub- and gradually (TO3). The two intensive sharp lines observed strates. The tetragonal-orthorhombic and orthorhombic- Ϫ1 at 173 and 187 cm at low temperatures ͑rhombohedral rhombohedral phase transitions of the bulk BaTiO3 are ͒ phase in the BaTiO3 crystal are broader and their relative completely absent in the thin films. intensity is much smaller in the films at all temperatures A quantitative determination of film strain is critical for measured. In fact, the positions, relative intensities and line the understanding of the complete disappearance of the low- shapes of phonon modes in the spectra of the films are simi- temperature phase transitions. Figure 3 shows the x-ray dif- lar to those of the orthorhombic phase of the single crystal. fraction results for a BaTiO3 /SrRuO3 film on LaAlO3 sub- ͓ ͑ ͔͒ The low-frequency range of Raman spectra ͑below strate. The -2 scan Fig. 3 a shows that both the BaTiO3 Ϫ1 100 cm ) both in the crystals and thin films contains a and SrRuO3 layers grow with c axis normal to the substrate. feature of overdamped soft mode characteristic for The full width at half maximum ͑FWHM͒ values in 2 are
174101-2 ABSENCE OF LOW-TEMPERATURE PHASE... PHYSICAL REVIEW B 69, 174101 ͑2004͒
FIG. 4. Domain stability map for BaTiO3 films under biaxial strain. The stars indicate the strain at room temperature, determined from x-ray data, and estimated strain at 5 K.
ϭ and c0 4.031 Å, slightly different from those in the bulk BaTiO . Then from the measured lattice constants of our FIG. 3. X-ray diffraction patterns of BaTiO3 /SrRuO3 film on 3 ͑ ͒ BaTiO films the strain in the thin film (aϪa )/a is calcu- LaAlO3 substrate: a –2 scan; LaAlO3 and SrRuO3 peaks are 3 0 0 marked by stars and down triangles, respectively. ͑b͒ scan of the lated to be a tensile strain of 0.58%. Although the lattice ϭ ͓ ͔ BaTiO3 202 peak; 0° is aligned to be parallel to the 100 constants of both LaAlO3 and SrTiO3 substrates and of in-plane direction of the LaAlO3 substrate. SrRuO3 buffer layer are smaller than that of BaTiO3, the misfit strain in our BaTiO3 film is fully relaxed at the growth ͒ temperatures because the film thickness far exceeds the criti- small (0.48° for the 003 peak . The scan of the BaTiO3 cal thickness for strain relaxation. The tensile strain arises ͓ ͑ ͔͒ 202 peak Fig. 3 b shows that the layers are also in-plane from the mismatch between the thermal expansions of ͓ ͔ ʈ͓ ͔ aligned with the substrate lattice with 100 BaTiO3 100 SrRuO and BaTiO . The thermal expansion coefficient of ʈ 3 3 SrRuO3 ͓100͔ LaAlO3. The FWHM value for the 202 peak ͑ ͒ 28 SrRuO3 Ref. 27 is smaller than that of BaTiO3. Conse- is 0.8°. From these scans, the lattice parameters of SrRuO3 quently, when a BaTiO3 film is cooled to room temperature buffer layers are determined to be aϭ3.93 Å for the after the deposition, it contracts more than the SrRuO3 buffer pseudocubic lattice, equal to the bulk values. For BaTiO3 layer and a tensile stress is imposed on it. It is difficult to ϭ ϭ films the measured lattice constants are a b 4.01 calculate the strain between 5–300 K because the thermal Ϯ0.01 Å ͑in plane͒, and cϭ4.004Ϯ0.001 Å. Similar x-ray expansion coefficients of SrRuO3 and BaTiO3 are tempera- diffraction data were obtained for films grown on SrTiO3 ture dependent. It is, however, expected to increase when the substrates, with the same epitaxial orientation, the same temperature is lowered, which is schematically indicated by SrRuO3 and BaTiO3 lattice constants, and even smaller the dotted line in Fig. 4. The end point at 5 K ͑0.78%͒ was FWHM values (0.38° for 003 and 0.7° for 202 peaks͒. estimated by assuming the average thermal mismatch be- TEM studies showed that the films consist of columnar tween SrRuO3 and BaTiO3 to be temperature independent. grains with in-plane grain size of ϳ100 nm. Electron dif- fraction results indicated the BaTiO grains are oriented in 3 III. PHASE-FIELD MODELING AND DOMAIN STABILITY same direction and confirmed the epitaxial relationship be- MAP tween the substrate, SrRuO3 and BaTiO3 layers described above. High crystalline quality of the films allows to pre- We compare the film strain with the calculated phase dia- sume the structural defects are unlikely to cause so drastic gram of BaTiO3 films under biaxial substrate constraint and changes in the phase transition behavior. Although we cannot over the temperature range from room temperature to 0 K. completely exclude other factors which can affect the film The widely quoted phase diagram of Pertsev et al.5 assumes properties, such as composition fluctuations or grain size ef- a simple domain structure with a particular domain wall ori- fect, we believe the strain in the thin films is the most prob- entation as a priori for the thermodynamic calculations. As a able cause of suppressing the low-temperature phase transi- matter of fact, for BaTiO3 films, three different domain sta- tions characteristic for bulk BaTiO3. bility maps have been presented, which depend on the as- From the material parameters of BaTiO3, the same as sumptions of possible domain states and domain-wall used by Koukhar et al. in Ref. 5, we calculated the lattice orientations.5 In this work, we used a phase-field approach in ϭ constants of the unstrained BaTiO3 film to be a0 3.987 Å which no priori assumption on the occurrence of certain
174101-3 D. A. TENNE et al. PHYSICAL REVIEW B 69, 174101 ͑2004͒ phases and domain structures is made.6 A brief description of cannot absolutely rule out the possibility of existing mixed the technique is as follows. phases in the films at room temperature, if the lower uncer- A ferroelectric domain structure is represented by the spa- tainty limit is taken into account. In the case of mixed ϭ tial distribution of the polarization field P (P1 ,P2 ,P3), and orthorhombic-tetragonal phases there would be 90° domain its temporal evolution is described by the time dependent walls at the boundaries between the orthorhombic phase with Ginzburg-Landau ͑TDGL͒ equations, in-plane polarization and tetragonal phase polarized along c axis. However, the TEM images of our films showed the P ͑x,t͒ ␦Fץ i ϭϪ ͑ ϭ ͒ ͑ ͒ existence of many antiphase boundaries, but no 90° domains L␦ ͒ i 1,2,3 , 1 ץ t Pi͑x,t were found. Therefore, the films are most likely to be in the orthorhombic phase with in-plane polarization. This is con- where L is the kinetic coefficient related to the domain-wall firmed by the Raman spectroscopy results ͑Figs. 1 and 2͒ as mobility, and F is the total free energy of the system. x the form of the spectra, the phonon line positions and relative ϭ (x1 ,x2 ,x3) is a rectangular coordinate system and x3 is intensities in the BaTiO film are close to those of the ortho- normal to the film-substrate interface. The total free energy 3 rhombic phase of BaTiO3 single crystal. Because the tensile of the system includes the bulk free energy, the domain wall strain in the film increases when the temperature decreases ͑ ͒ energy, and the elastic energy. The TDGL equation 1 was from room temperature to 5 K, as schematically shown by solved numerically using the semi-implicit Fourier-spectral 29 the dotted line in Fig. 4, the film remains in the orthorhombic method for the time-stepping and spatial discretization. The phase without undergoing any phase transitions. macroscopic constraint of the substrate to the film is mea- ¯⑀ ϭ¯⑀ ϭ⑀ ¯⑀ ϭ sured by the average strain 11 22 and 12 0. The con- IV. CONCLUSION tinuities of the deformation and stresses on film-substrate interface provide the local/microscopic constraint. By evolv- In summary, Raman spectroscopy in BaTiO3 films grown ing Eq. ͑1͒ from an initial paraelectric state with small ran- by pulsed laser deposition on SrTiO3 and LaAlO3 substrates dom perturbation, the thermodynamically stable domain covered by SrRuO3 buffer layers in the temperature range 5–325 K shows that the phase transitions between the differ- structure and phase at given temperature and substrate con- ͑ straint strain were obtained at the end of a simulation. ent ferroelectric phases tetragonal-orthorhombic- ͒ rhombohedral , characteristic of bulk BaTiO3, are com- The results of the phase-field calculation for BaTiO3 films are plotted in a phase diagram ͑domain stability map͒ shown pletely absent in the films. This behavior is explained by the in Fig. 4. The phases represented by various symbols are: presence of tensile strain in the BaTiO3 film caused by ther- solid squares—the tetragonal phases with polarization nor- mal mismatch with the underlying SrRuO3 layer. This is con- ϭ Ϯ firmed by a thin film phase diagram obtained from a phase- mal to the film/substrate interface, i.e., P (0,0, Pt); solid diamonds—the distorted rhombohedral phase with P field calculation, which predicts, without any priori ϭ Ϯ Ϯ Ϯ assumption, that the thermodynamically stable phase be- ( P1r , P1r , P3r); solid circles—the orthorhombic phase with polarization parallel to the film/substrate inter- tween 5–325 K for the experimentally determined value of ϭ Ϯ Ϯ tensile strain in the BaTiO3 film is an orthorhombic phase face, i.e., P ( Po , Po,0); and half-solid circles—the dis- ϭ Ϯ Ϯ with in-plane polarization. torted orthorhombic phase with P ( P1o,0, P3o)orP ϭ Ϯ Ϯ (0, P1o , P3o). There are small regions, represented by ACKNOWLEDGMENTS crossed diamonds, within which more than one ferroelectric phases coexist. From this domain map, it is seen that a ten- This work was partially supported by DOE under Grant ͑ ͒ sile strain of 0.58% places a BaTiO3 film in the orthorhom- Nos. DFFG02-84ER45095 Xi and DE-FG02-97ER45638 bic phase at room temperature. It should be noted that the ͑Schlom͒, by DARPA under Grant No. MDA972-01-1-0023 uncertainty in the in-plane lattice parameter determined prom ͑Xi͒ and by NSF under Grant Nos. DMR-01-22638 ͑Chen͒, x-ray data ͑0.01 Å͒ produces rather large uncertainty in the DMR-9875405 and DMR-9871177 ͑Pan͒, and DMR- value of strain ͑shown by error bars on Fig. 4͒. Therefore, we 0103354 ͑Schlom, Chen, Xi͒.
*Electronic address: [email protected]; On leave from the Institute of 3722 ͑2000͒; V.G. Koukhar, N.A. Pertsev, and R. Waser, Phys. Semiconductor Physics, Novosibirsk, Russia. Rev. B 64, 214103 ͑2001͒. †Currently at Tsinghua University, Beijing, China. 6 Y.L. Li, S.Y. Hu, Z.K. Liu, and L.Q. Chen, Appl. Phys. Lett. 78, ‡Current address: Solid State Physics Laboratory, Lucknow Rd., 3878 ͑2001͒; Y.L. Li, S.Y. Hu, Z.K. Liu, and L.Q. Chen, Acta Timarpur, Delhi-110 054, India. Mater. 50, 395 ͑2002͒; L.Q. Chen, Annu. Rev. Mater. Research 1 S.K. Streiffer, C. Basceri, C.B. Parker, S.E. Lash, and A.I. Kin- 32,113͑2002͒. gon, J. Appl. Phys. 86, 4565 ͑1999͒. 7 Y. Yoneda, T. Okabe, K. Sakaue, H. Terauchi, H. Kasatani, and K. 2 D. Damjanovic, Rep. Prog. Phys. 61, 1267 ͑1998͒. Deguchi, J. Appl. Phys. 83, 2458 ͑1998͒. 3 J.S. Speck and W. Pompe, J. Appl. Phys. 76, 466 ͑1994͒ 8 C.L. Li, Z.H. Chen, D.F. Cui, Y.L. Zhou, H.B. Lu, C. Dong, F. 4 A.L. Roytburd, J. Appl. Phys. 83, 228 ͑1998͒; S.P. Alpay and A.L. Wu, and H. Chen, J. Appl. Phys. 86, 4555 ͑1999͒. Roytburd, ibid. 83, 4714 ͑1998͒. 9 K. Iijima, T. Terashima, K. Yamamoto, K. Hirata, and Y. Bando, 5 N.A. Pertsev, A.G. Zembilgotov, and A.K. Tagantsev, Phys. Rev. Appl. Phys. Lett. 56, 527 ͑1990͒. Lett. 80, 1988 ͑1998͒; N.A. Pertsev and V.G. Koukhar, ibid. 84, 10 Yu. I. Yuzyuk, P. Simon, I.N. Zakharchenko, V.A. Alyoshin, and
174101-4 ABSENCE OF LOW-TEMPERATURE PHASE... PHYSICAL REVIEW B 69, 174101 ͑2004͒
E.V. Sviridov, Phys. Rev. B 66, 052103 ͑2002͒; Yu.I. Yuzyuk, ͑1982͒; H. Presting, J.A. Sanjurjo, and H. Vogt, ibid. 28, 6097 J.L. Sauvajol, P. Simon, V.L. Lorman, V.A. Alyoshin, I.N. Za- ͑1983͒. kharchenko, and E.V. Sviridov, J. Appl. Phys. 93, 9930 ͑2003͒. 20 J.D. Freire and R.S. Katiyar, Phys. Rev. B 37, 2074 ͑1988͒. 11 S.K. Streiffer, J.A. Eastman, D.D. Fong, C. Thompson, A. 21 U.D. Venkateswaran, V.M. Naik, and R. Naik, Phys. Rev. B 58, Munkholm, M.V. Ramana Murty, O. Auciello, G.R. Bai, and 14 256 ͑1998͒. G.B. Stephenson, Phys. Rev. Lett. 89, 067601 ͑2002͒. 22 M. Osada, M. Kakihana, S. Wada, T. Noma, and W.-S. Cho, Appl. 12 D. Fuchs, C.W. Schneider, R. Schneider, and H. Rietschel, J. Phys. Lett. 75, 3393 ͑1999͒. Appl. Phys. 85, 7362 ͑1999͒. 23 L.H. Robins, D.L. Kaiser, L.D. Rotter, P.K. Schenck, G.T. Stauf, 13 O. Tikhomirov, H. Jiang, and J. Levy, Phys. Rev. Lett. 89, 147601 and D. Rytz, J. Appl. Phys. 76, 7487 ͑1994͒. ͑2002͒. 24 V.I. Merkulov, J.R. Fox, H.-C. Li, W. Si, A.A. Sirenko, and X.X. 14 Landolt-Bo¨rnstein Numerical Data and Functional Relationships Xi, Appl. Phys. Lett. 72, 3291 ͑1998͒. in Science and Technology, New Series, Group III ͑Springer, 25 D.A. Tenne, A.M. Clark, A.R. James, K. Chen, and X.X. Xi, Berlin, 1981͒,p.67,Vol.16. Appl. Phys. Lett. 79, 3836 ͑2001͒. 15 C.H. Perry and D.B. Hall, Phys. Rev. Lett. 15, 700 ͑1965͒. 26 R.E. Cohen and H. Krakauer, Phys. Rev. B 42, 6416 ͑1990͒. 16 M. DiDomenico, Jr., S.P.S. Porto, and S.H. Wemple, Phys. Rev. 27 J.J. Neumeier, A.L. Cornelius, and K. Andres, Phys. Rev. B 64, Lett. 19, 855 ͑1967͒; M. DiDomenico, Jr., S.H. Wemple, S.P.S. 172406 ͑2001͒. Porto, and R.P. Bauman, Phys. Rev. 174, 522 ͑1968͒. 28 Y.S. Touloukian, R.K. Kirby, R.E. Taylor, and T.Y.R. Lee, Ther- 17 A. Scalabrin, A.S. Chaves, D.S. Shim, and S.P.S. Porto, Phys. mal Expansion, Nonmetallic Solids, Thermophysical Properties Status Solidi B 79, 731 ͑1977͒. of Matter ͑Plenum, New York, 1977͒, Vol. 13. 18 G. Burns and F.H. Dacol, Phys. Rev. B 18, 5750 ͑1978͒. 29 L.Q. Chen and J. Shen, Comput. Phys. Commun. 108, 147 19 H. Vogt, J.A. Sanjurjo, and G. Rossbroich, Phys. Rev. B 26, 5904 ͑1998͒.
174101-5 Überblick
Festkörperphysik
Oxide – Tausendsassas für die Elektronik Oxide bieten eine spektakuläre Vielfalt an funktionalen Eigenschaften Jochen Mannhart und Darrell G. Schlom
Dank gewaltiger Fortschritte der Dünnschicht- technik ist es seit einem Jahrzehnt möglich, aus einer Vielzahl von Materialien hochwertige Schichten und Multilagenstrukturen zu wach- sen. Die für elektronische Bauelemente zur Verfügung stehende Materialpalette hat sich da- durch wesentlich verbreitert und ist inzwischen außerordentlich reichhaltig. Die mit Oxiden erzielbaren Funktionalitäten sind sehr vielfältig Ultraschall-Bild und umfassen mit Piezoelektrizität, Sensor- eines Babies im eigenschaften und elektrooptischem Verhalten Mutterleib, aufge- sogar die Übergangsbereiche zur Mechanik, nommen mithilfe von Ultraschall- Optik, Chemie und Biologie. Die Halbleiter- sendern und -emp- industrie nutzt schon jetzt funktionale Oxide, fängern aus dem beispielsweise Ferroelektrika für nichtflüchtige Oxid Pb(Zr,Ti)O3. Speicher, und im Jahr 2007 sollen in den Tran- Die Aufnahme wurde digital bear- sistoren der Laptops Hafnium-basierte Oxide beitet (Quelle: GE als Gate-Dielektrika eingesetzt werden. Healthcare Techno- logies.)
in aktuelles, großes Problem der Halbleiter- Gate-Barrieren-Stapel aus SiO2 and Oxinitrid (Si-O- industrie ist die begrenzte Leistungsfähigkeit N) aber nicht mehr, falls t etwa 1 nm unterschreitet,
von SiO2, das als Gate-Isolator in den Si- da dann der durch die SiO2-Barriere tunnelnde Strom E 2 MOSFETs (Metall-Oxide-Semiconductor Field-Effect eine Dichte von über 1000 A/cm aufweist. Die Tunnel- Transistors) verwendet wird, den „Arbeitspferden“ der stromdichte steigt für jedes Zehntel Nanometer, um das meisten Computerchips. Falls bis 2007 keine alternati- die Barrierendicke verringert wird, um eine Größen- ven Materialien für die Gate-Isolatoren zur Verfügung ordnung an und führt zu unakzeptablem Energiever- stehen, verlangsamt sich die gegenwärtig phänomenale brauch, zu Kühlproblemen sowie zu einer unerwünsch- Entwicklung der Computer-Leistung massiv, da die In- ten Rückwirkung des DS-Kanals auf das Gate. Dieses tegrationsdichte der Transistoren dann nicht mehr dem Problem lässt sich nur mit Materialien als Gate-Isolato- berühmten Mooreschen „Gesetz“ folgen könnte [1]. ren lösen, welche schon bei großen Barrierendicken die Ein Blick auf einen MOSFET eines aktuellen Pen- vorteilhaften dielektrischen Eigenschaften ultradünner tium-4-Prozessors verdeutlicht die SiO2-Schichten aufweisen. Wie Physik, die diesem Problem zugrun- Kompakt unten beschrieben, sollen deswe- de liegt (Abb. 1). In diesen MOS- gen in den kommenden Jahren in Prof. Dr. Jochen Oxide sind als gute Isolatoren unersetz- FETs ist die nitrierte SiO2-Schicht, den hochintegrierten Schaltungen Mannhart, Bereich welche die Gate-Kontakte von den lich für die Mikroelektronik, weisen aber Hafnium-basierte Oxide als Gate- Elektronische darunterliegenden Drain-Source- auch zahlreiche andere interessante Eigen- Isolatoren mit hohen Dielektrizi- Korrelationen und Magnetismus (EKM), Kanälen isoliert, nur 1,2 nm dick. schaften auf. tätskonstanten eingesetzt werden. Institut für Physik, Um die FETs mit möglichst geringer Für Bauelemente stehen heute oxidische Dieser Einsatz komplexer Oxide Universität Augs- Arbeitsspannung betreiben zu kön- Schichtsysteme aus u. a. magnetischen, ferro- in der Elektronik ist kein Einzelfall. burg, 86135 Augs- nen und dennoch im Drain-Sour- elektrischen, piezoelektrischen und sog. Für Anwendungen in der Elektronik burg; Prof. Darrell multiferroischen Materialien sowie Hoch- G. Schlom, PhD, Ma- ce-(DS)-Kanal die für den Betrieb sind die Oxide sogar ausgesprochen terials Science and notwendige Ladungsträgerdichte zu temperatur-Supraleitern zur Verfügung. beliebt. Woran liegt dies? Zum gro- Engineering, Penn erhalten, muss das Verhältnis von Die vielfältigen Anwendungen reichen ßen Teil lässt sich das auf die große State University, 108 Dielektrizitätskonstante und Dicke von Gate-Barrieren für die Mikroelektronik Elektronegativität des Sauerstoffs, Materials Research der Gate-Barriere t möglichst groß, über supraleitende Filter für den Mobil- den im Vergleich zu anderen An- Institute Building, funk, FRAM-Speicher, Infrarotdetektoren 2– University Park, PA t also möglichst klein sein. Ande- ionen kleinen Radius des O -Ions 16802-5005, USA rerseits funktionieren die heutigen bis hin zu selbstreinigenden Fliesen. und seine doppelte Ladung zurück-
Physik Journal © 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim 1617-9439/05/0606-45 4 (2005) Nr. 6 45 Überblick führen. Aufgrund dieser besonderen Eigenschaften angeregt werden können. Aufgrund der Spin-Ordnung lassen sich viele Oxide, selbst solche mit mehreren sind die Oxide wichtige Kandidaten für die Spintronik Kationensorten, unter praktikablen Synthesebedingun- [2, 3], deren Bauelemente, z. B. Leseköpfe von Festplat- gen in thermodynamisch stabilen Phasen herstellen. tenlaufwerken, den Spin der Elektronen nutzen. Ob sich Da diese, auch wegen der hohen Polarisierbarkeit des sogar, wie im Konzept der „Orbitronic“ vorgeschlagen Sauerstoffs, zudem eine außerordentliche Vielfalt elek- [3], Anregungen der Orbital-Anordnung im Kristallgitter tronischer Eigenschaften aufweisen, sind sie besonders für Bauelemente nutzen lassen, wird die Zukunft zeigen. attraktiv für Bauelemente. Oxide haben aufgrund ihrer großen Energielücken Aufgrund der Größe und der Bindungseigenschaften das Zeug zu hervorragenden Isolatoren. Lange Zeit der Sauerstoffionen sowie des breiten Spektrums der zur begnügte man sich damit, sie als solche in elektrischen Verfügung stehenden Metallionen kristallisieren die Oxi- Komponenten einzusetzen. Die elektronischen Eigen- de in einer Fülle von Strukturen (s. schaften der Oxide lassen sich jedoch durch Dotieren Info kasten „Kristallstrukturen“). Sie dramatisch verändern, wobei sehr hohe Dotierkonzen- bilden unter anderem Rutile (z. B. trationen von mehreren 1021/cm3 gebräuchlich sind [4].
CrO2) und Spinelle (z. B. Fe3O4), Dotierte Oxide können ausgesprochen gute Leiter sein. häufig auch Perowskite (z. B. Einen ganz besonderen Fall bilden hierbei die Kuprat-
BaTiO3) oder perowskitverwandte Supraleiter, die bei Temperaturen bis zu 135 K, unter Gitter (z. B. YBa2Cu3O7). Die Man- Druck sogar bis 153 K, supraleitend sind. Die hohen nigfaltigkeit der Kristallstrukturen, Dotierungen bringen, im Gegensatz zu Halbleitern, die Tendenz des Sauerstoffs, um nicht nur Ladungsträger in Bänder ein, sondern ändern Metallionen mit d-Valenzelektronen zugleich auch die spektrale Dichte der Elektronen, Oktaeder oder Tetraeder auszubil- die Korrelationsparameter und die Suszeptibilitäten. den, sowie die Variationsmöglich- Dadurch können Phasenübergänge, zum Beispiel Me- keiten der Sauerstoffbesetzung und tall-Isolator-Übergänge, induziert werden. Somit lassen der Dotierung bilden sich in vielför- sich die Eigenschaften solcher Oxide in extrem weiten Abb. 1: migen elektronischen Eigenschaften Bereichen durch Ändern ihrer Ladungsträgerdichte Transmissions-Elektronenmikroskopische (TEM) Aufnahme des Querschnitts eines ab. Dieser Vielfalt übergeordnet ist variieren, was nicht nur durch Dotierung, sondern bei- CMOS-Transistors in einem Pentium-4- das charakteristische Auftreten von spielsweise auch durch elektrische Felder oder durch Prozessor (Quelle: Intel Corporation). Energielücken im Elektronenvolt- Bestrahlung mit Licht möglich ist (s. u.). Bereich. Diese werden durch den Da Sauerstoff-Defekte die Oxide in der Regel dotie- Einfluss der elektrostatischen Energie (Madelung-Ener- ren und Sauerstoff zudem leicht diffundiert, wird in der gie) und der Energien der chemischen Bindung auf die Oxidelektronik der Kontrolle der Sauerstoffkonzentrati- elektronischen Energien der freien Metall- und Sauer- on viel Aufmerksamkeit gewidmet. Die einzelnen Oxide stoff-Ionen erzeugt. Bestimmt durch die Energielücken verhalten sich hierbei sehr unterschiedlich. In manchen sind eine größere Zahl der Oxide Bandisolatoren. Die Manganat-Schichten lässt sich der Sauerstoffgehalt bei- elektronischen und magnetischen Eigenschaften vieler spielsweise oft nur mit Aufwand präzise einstellen. Die
Oxide sind auch durch starke elektronische Korrelati- Hoch-Tc-Supraleiter sind hier robuster und zudem lang- onen geprägt (s. Infokasten „Elektronische Korrelatio- zeitstabil; epitaktische Schichten, wie sie zum Beispiel nen“). Diese Korrelationen können in den komplexen in Hochfrequenzfiltern eingesetzt werden, erfüllen sogar Kristallstrukturen Ladungsordnungen und sogar Ord- die hohen Anforderungen der Raumfahrttechnik. nungen des Spinsystems oder der Orbitale im Kristall- Wenngleich von den ionischen Materialien bisher gitter generieren, wobei auch Quasiteilchenzustände bevorzugt Oxide für elektronische Anwendungen ge-
Kristallstrukturen