Combinatorial approach for ferroelectric material libraries ... - PNAS

0 downloads 0 Views 2MB Size Report
Jan 23, 2007 - combinatorial method depends on the preparation of a reliable high-density ... trafine mist is generated with nebulizer (frequency: 1.65 MHz).
Combinatorial approach for ferroelectric material libraries prepared by liquid source misted chemical deposition method Ki Woong Kim*, Min Ku Jeon*, Kwang Seok Oh*, Tai Suk Kim*, Yun Seok Kim†, and Seong Ihl Woo*‡ *Department of Chemical and Biomolecular Engineering and Center for Ultramicrochemical Process Systems, and †Department of Materials Science and Engineering, Korea Advanced Institute of Science and Technology, 373-1 Guseong-dong, Yuseong-gu, Daejeon 305-701, Republic of Korea Communicated by Gabor A. Somorjai, University of California, Berkeley, CA, November 24, 2006 (received for review August 24, 2006)

Combinatorial approach for discovering novel functional materials in the huge diversity of chemical composition and processing conditions has become more important for breakthrough in thin film electronic and energy-conversion devices. The efficiency of combinatorial method depends on the preparation of a reliable high-density composition thin-film library. The physico-chemical properties of each sample on the library should be similar to those of the corresponding samples prepared by one-by-one conventional methods. We successfully developed the combinatorial liquid source misted chemical deposition (LSMCD) method and demonstrated its validity in screening the chemical composition of Bi3.75LaxCe0.25-xTi3O12 (BLCT) for high remanent polarization (Pr). LSMCD is a cheap promising combinatorial screening tool. It can control the composition up to ppm level and produce homogeneous multicomponent library. LSMCD method allows us to prepare BLCT thin-film library at the variation of 0.4 mol% of La. Maximum 2Pr is 35 ␮C/cmⴚ2 at x ⴝ 0.21. The intensity of (117) XRD peak is quantitatively related to 2Pr. Newly developed scanning piezoelectric deformation measurement for nano-sized samples using scanning probe microscope (SPM) is also found out to be reliable for determining the relative ranking of Pr value rapidly. BLCT thin film array 兩 microbeam x-ray diffraction 兩 scanning probe microscope

C

ombinatorial methods consist of library design, preparation of library, evaluation of its functional properties, and its physicochemical characteristics. These methods are considered to be very efficient and quickly provide us the quantitative activity structure relationship to establish the criteria for designing novel functional materials, which have been extensively used for the discovery of functional electronic thin film materials, novel polymers, and catalysts since the mid-1990s (1–10). Thin-film combinatorial libraries have been synthesized mainly by RF-sputtering, pulsed laser ablation, and chemical vapor deposition. Thin-film libraries containing various compositions were fabricated by adopting multimask techniques (3–6). Ten binary masks or five quaternary masks should be used to prepare 1,024 (210) different samples on a wafer. These mask operations complicated the operation of the instrument as well as creating the contamination problems. These problems were solved by adopting only two masks, one x–y movable shutter and one shadow mask (7–11). The x–y movable shutter moves at a constant velocity in a linear or stepwise way to control the deposition time in one direction of wafer. The shadow mask containing regularly spaced holes is located just below the moving shutter, which fabricates a high-density discrete library in a very simple manner. Recently, off-axis codeposition method was reported to prepare thin-film metal alloy library (12, 13). By controlling the distance between the target guns and the substrate, multitarget materials can be deposited on the substrate with concentration gradient simultaneously. This method does not need to use a movable shutter or many masks. Only shadow mask is necessary when the fabrication of discrete library is

1134 –1139 兩 PNAS 兩 January 23, 2007 兩 vol. 104 兩 no. 4

preferred. However, this method cannot control the thickness of each sample on the library and is not suitable when the physicochemical properties vary with thickness. The library fabrication methods discussed above need to prepare single- or multicomponent target materials, which is sometimes very expensive, and it can be difficult to obtain a uniform composition. It is also very difficult to control the deposition rate of each component in multicomponent target arising from the difference in its inherent sputtering property. Sputtering rate depends on many variables, so it is too difficult to control precisely. Therefore, accurate control of composition within ⫾3% in the library cannot be obtained. To overcome these problems, some research groups applied the wet process including spray pyrolysis and ink-jet technique for combinatorial library synthesis. However, these processes did not accomplish a finely compositional tuned thin-film library because of their inherent weakness (14–17). Therefore, in our study, liquid source misted chemical deposition (LSMCD) method was developed to prepare a thin-film library. Adequate amount of metal precursor compounds was dissolved in suitable solvent. This mixed solution is equivalent to the target in the sputtering method. Composition control is more accurate with LSMCD than with conventional sputtering, because exact weighing can be obtained even in ppm level. Furthermore, the cost of metal precursor solution (a few dollars) is just the cost of chemicals and much cheaper than target (average cost for 4-inch target is ⬇$5,000). Hundreds of targets are required for the fabrication of a very diverse library. The material cost to fabricate a thin-film library with LSMCD is negligible compared with conventional sputtering method. Ultrafine mist is generated with nebulizer (frequency: 1.65 MHz) and transported by Ar flow to the substrate. Suitable thermal treatments are required to obtain multicomponent metal alloy or metal oxide thin-film (18–23). The variation of deposition time with x–y movable shutter will generate the concentration gradient in a similar manner to combinatorial sputtering. It is possible to fabricate homogeneous liquid-phase combinatorial library by liquid-phase mixing in LSMCD, whereas in sputtering method, each component of library have to be mixed by solid-state diffusion, resulting in the formation of more homogeneous multicomponent compound with LSMCD. It can be concluded Author contributions: K.W.K. and S.I.W. designed research; K.W.K. and T.S.K. performed research; K.W.K., K.S.O., and Y.S.K. contributed new reagents/analytic tools; K.W.K., M.K.J., and S.I.W. analyzed data; and K.W.K. and S.I.W. wrote the paper. The authors declare no conflict of interest. Freely available online through the PNAS open access option. Abbreviations: LSMCD, liquid source misted chemical deposition; BLT, (Bi,La)4Ti3O12; BCT, Bi4-xCexTi3O12; BLCT, Bi3.75LaxCe0.25-xTi3O12; P–E, polarization-electric field; SPM, scanning probe microscope, AFM, atomic force microscope. ‡To

whom correspondence should be addressed. E-mail: [email protected].

This article contains supporting information online at www.pnas.org/cgi/content/full/ 0610146104/DC1. © 2007 by The National Academy of Sciences of the USA

www.pnas.org兾cgi兾doi兾10.1073兾pnas.0610146104

0.25/0

0.228/0.022

0.207/0.043

0.185/0.065

0.163/0.087

0.142/0.108

0.12/0.13

0.098/0.152

(0.25/0)

(0.228/0.022)

(0.207/0.043)

(0.184/0.066)

(0.16/0.09)

(0.143/0.107)

(0.12/0.13)

(0.098/0.152)

188nm

189nm

192nm

195nm

196nm

194nm

193nm

188nm

0.236/0.014

0.214/0.036

0.193/0.057

0.171/0.079

0.149/0.101

0.128/0.122

0.106/0.144

0.084/0.166

(0.231/0.019)

(0.215/0.035)

(0.194/0.056)

(0.171/0.079)

(0.147/0.103)

(0.128/0.122)

(0.109/0.141)

(0.084/0.166)

192nm

191nm

195nm

200nm

200nm

205nm

200nm

196nm

0.222/0.028

0.2/0.05

0.178/0.072

0.157/0.093

0.135/0.115

0.114/0.136

0.092/0.158

0.07/0.18

(0.223/0.027)

(0.2/0.05)

(0.176/0.074)

(0.159/0.091)

(0.133/0.117)

(0.114/0.136)

(0.089/0.161)

(0.07/0.18)

191nm

194nm

196nm

202nm

204nm

208nm

201nm

195nm

0.208/0.042

0.186/0.064

0.164/0.086

0.143/0.107

0.121/0.129

0.099/0.151

0.078/0.172

0.056/0.194

(0.201/0.049)

(0.189/0.061)

(0.162/0.088)

(0.143/0.107)

(0.121/0.129)

(0.095/0.155)

(0.078/0.172)

(0.056/0.194)

194nm

196nm

197nm

203nm

207nm

205nm

200nm

195nm

0.194/0.056

0.172/0.078

0.15/0.1

0.129/0.121

0.107/0.143

0.085/0.165

0.064/0.186

0.042/0.208

(0.195/0.055)

(0.172/0.078)

(0.151/0.099)

(0.129/0.121)

(0.111/0.139)

(0.088/0.162)

(0.068/0.182)

(0.042/0.208)

197nm

195nm

197nm

202nm

201nm

204nm

200nm

194nm

0.18/0.07

0.158/0.092

0.136/0.114

0.115/0.135

0.093/0.157

0.071/0.179

0.05/0.2

0.028/0.222

(0.183/0.067)

(0.158/0.092)

(0.136/0.114)

(0.115/0.135)

(0.093/0.157)

(0.076/0.174)

(0.052/0.198)

(0.028/0.222)

190nm

192nm

197nm

200nm

206nm

200nm

196nm

192nm

0.166/0.084

0.144/0.106

0.122/0.128

0.101/0.149

0.079/0.171

0.057/0.193

0.036/0.214

0.014/0.236

(0.165/0.085)

(0.146/0.104)

(0.118/0.132)

(0.103/0.147)

(0.079/0.171)

(0.057/0.193)

(0.036/0.214)

(0.014/0.236)

192nm

194nm

196nm

196nm

201nm

200nm

195nm

186nm

0.152/0.098

0.13/0.12

0.108/0.142

0.087/0.163

0.065/0.185

0.043/0.207

0.022/0.228

0/0.25

(0.153/0.097)

(0.13/0.12)

(0.112/0.148)

(0.09/0.16)

(0.065/0.185)

(0.047/0.203)

(0.02/0.23)

(0/0.25)

187nm

190nm

193nm

194nm

198nm

196nm

193nm

184nm

The top number is the chemical composition Ce兾La designed to prepare. The blanket value is chemical composition (Ce兾La) measured by wavelength dispersive spectroscopy, and the bottom number is the thickness measured by SEM cross-section image. The location of 8 ⫻ 8 samples on the wafer is equivalent to the location of the square in the table.

that LSMCD enables us to fabricate a thin-film library with more accurate control of composition at a low cost. The price of the LSMCD instrument is also much cheaper than a conventional PVD instrument. These advantages of LSMCD over the sputtering method will be demonstrated in fabricating a multicom ponent ferroelectric oxide thin-film library for the optimization of its electrical properties as follows. Bismuth layered perovskite ferroelectric films, SrBi2Ta2O9 (24) and (Bi,La)4Ti3O12 (BLT) (25), have drawn an interest for ferroelectric random-access memory because of the fatigue-free properties over 1010 read/write cycles, unlike Pb(Zr,Ti)O3 (PZT) (26, 27). However, these materials have smaller remanent poKim et al.

larization (Pr) than PZT. Numerous attempts have been made to substitute lanthanide group atoms for Bi site in Bi4Ti3O12 to increase Pr. It is well known that Pr depends on the chemical composition (28–32). Woo et al. first presented that the Cesubstituted Bi4Ti3O12 (Bi4-xCexTi3O12, BCT) fabricated by LSMCD exhibited large Pr (2Pr ⫽ 20 ␮Ccm⫺2 at 15 V, x ⫽ 0.75) and superior fatigue endurance over 1010 read/write cycles. Hence, it is quite reasonable to substitute Bi with both La and Ce for ferroelectric materials with high Pr. We prepared a Bi3.75LaxCe0.25-xTi3O12 (BLCT, 0 ⱕ x ⱕ 0.25) array with LSMCD to optimize the chemical composition for high Pr. Bi/(La⫹Ce) was fixed to 3.75/0.25 because our structural refinement result PNAS 兩 January 23, 2007 兩 vol. 104 兩 no. 4 兩 1135

CHEMISTRY

Table 1. The chemical composition and thickness of BLCT thin-film 8 ⴛ 8 samples on a library annealed at 700°C for 1 h in furnace under oxygen atmosphere

with BLT powder showed that 2Pr was maximum when the Bi/La was 3.75/0.25 (33–35). Fast and reliable measurement of remanent polarization is required for a continuous-gradient library and ␮m-sized samples, because the conventional capacitance–voltage measurement method is slow and requires electronic wirings. A scanning probe microscope can apply the electric field into the sample on the library using a sharpened tip as movable electrode and induced the local deformation of samples. The correlation between the deformation length and Pr will be discussed. Results and Discussion A BLCT (0 ⱕ x ⱕ 0.25) thin-film library was prepared by LSMCD equipped automated shutter (see Materials and Methods). The chemical compositions of BLCT samples on the library measured by wavelength dispersive spectroscopy and thickness of each sample on the library measured by SEM cross-section image are shown in Table 1. The real location of each sample is the same as the location of Table 1; 0.4 mole % of composition control was demonstrated, as was the uniform thickness, 200 nm ⫾ 10 nm. Microbeam XRD (Bruker D8Discover GADDS CS, CuK␣ radiation) was used to analyze the crystallographic phase and the preferred orientation of as-prepared thin-film samples on the library. As shown in Fig. 1A, regardless of La content, all samples on the library did not contain any peaks of Bi2O3 (2␪ ⫽ 28°) as well as the other mono- or binary oxides. The intensities of (00n) peaks arising from the extent of crystal growth in the c-axis did not change much with La content. However, the intensity of (117) diffraction peak normalized by that of Pt (111) peak increased with the increase of La content and reached the maximum at x ⫽ 0.21. In the region of x ⬎ 0.21, there was an abrupt decrease in this value (Fig. 1B). The crystal growth in the a and b axes is related to the intensity of the (117) peak. It was reported that Ps is only 4␮Ccm⫺2 along the c axis and 50␮Ccm⫺2 along the a or b axis measured with single crystal Bi4Ti3O12 (36, 37). It was also reported that BLT thin film grown epitaxially in the a axis shows a Ps value of 50␮Ccm⫺2 (39). Fig. 2 shows the grain size as a function of La content. As shown in Fig. 2, grain size increased with the increase of La content. This can improve the ferroelectricity and piezoelectricity of thin film. However, when the grain size is larger than 0.5 ␮m, the leakage current density increased from 10⫺6 Acm⫺2 to 10⫺4 Acm⫺2 in the region of La content greater than x ⫽ 0.2 [supporting information (SI) Fig. 7]. 2Pr values of 64 samples on the library were measured by the polarization-electric field (P–E) hysteresis loop at the applied voltage of 10 V as shown in Fig. 3. Twenty Pt top electrodes (diameter: 200 ␮m) were deposited on each sample with a thickness of 100 nm to form a capacitor structure. As shown in Fig. 3, the remanent polarization increased with La content (x). Pr increased slightly between x ⫽ 0 and x ⫽ 0.10, but increased rapidly between x ⫽ 0.1 and 0.21 and decreased abruptly when x ⬎ 0.23. The maximum 2Pr (35␮Ccm⫺2) was obtained at x ⫽ 0.21. Comparing Fig. 3 with Fig. 1B, it could be suggested that the remanent polarization of BLCT thin film (x ⱕ 0.25) be quantitatively related to the intensity of normalized (117) peak of BLCT. As shown in Fig. 4, BLCT shows excellent fatigue endurance at 5 V over 1010 cycles similar to other bismuthlayered ferroelectrics. This is due to the increase of structural stability induced by decrease in oxygen vacancies (27). Substitution of La and Ce by Bi decreased oxygen vacancies. Four different BLCT thin films (2 cm ⫻ 2 cm) were fabricated at the same deposition condition of combinatorial experiments discussed above by varying the La content (x ⫽ 0.1, 0.12, 0.18, and 0.194) to validate the results obtained by 8 ⫻ 8 thin film samples on the combinatorial library. As shown in Fig. 5A, the 1136 兩 www.pnas.org兾cgi兾doi兾10.1073兾pnas.0610146104

Fig. 1. Structural analysis by using microbeam XRD. (A) The microbeam XRD spectra of BLCT thin-film samples on the library as a function of La content, x. These spectra were scanned in ␪-2␪ type in the range of 2␪ between 17 and 53° under Cu K␣ radiation (␭ ⫽ 1.5405 Å) operated at 40 mA and 40 kV. (B) The ratio of instensity of (117) diffraction peak to that of Pt (111) diffraction peak was plotted as a function of La content.

change in 2Pr of BLCT individually prepared shows a similar trend to that of BLCT prepared combinatorially, even though 2Pr values of individual BLCT samples are slightly lower than those of combinatorial array. Such an offset was due to the thickness difference between samples on the library (200 nm) and individual (240 nm) samples. Even though we need the identical deposition condition for samples on the library and individual samples, it was found out that the deposition rate for samples on the library is 20% lower than that for individual samples, because the shadow mask and moving mask used for samples on the library hindered mist from approaching to the substrate. As shown in Fig. 5 B and C, the thicker individual ferroelectric samples have larger coercive field strength. As the coercive field strength increases, higher electric field strength is necessary to switch the dipole moment. Therefore, individual sample shows a little bit lower remanent polarization than sample on the library at the same chemical composition under the same electric field (500 kV/cm). This result indicated that the preparation of combinatorial array by LSMCD is reliable for the rapid screening of various thin films having various compositions. A high-density library (2,500 samples per wafer) can be prepared with existing technology, but the evaluation of their Kim et al.

Fig. 2. The surface morphology of BLCT thin films as a function of La content. x ⫽ 0 (A), x ⫽ 0.08 (B), x ⫽ 0.17 (C), and x ⫽ 0.25 (D) in BLCT combinatorial array, where x is La content.

functional properties and physico-chemical characterization cannot be performed rapidly as long as conventional procedures are used. It is quite desirable to develop new methods to measure the functional and physico-chemical properties of thin film samples on the library in a very short time. P–E hysteresis characteristics of BLCT thin film samples on the library were measured with RT66A ferroelectric tester one by one for each sample in this study, which is not suitable for high-density library because it takes long time. Here, we report the rapid indirect measurement of remanent polarization by measuring the piezo-electricity with a scanning probe microscope (SPM) (38–41). A metallic atomic force microscope (AFM) tip was used as a movable top electrode. The electric field was applied between the AFM tip and bottom electrode to measure the topographic deformation. A laser beam was focused on the upper side of the cantilever and photo-diode to detect an angle of reflection from the upper side

of the cantilever as shown in Fig. 6A. This enables us to detect the topographic deformation up to 0.01 nm. Local piezoresponse on the grain or domain of each samples on the library can be detected nondestructively just by scanning the AFM tip using XY manipulator, which will decrease the measurement time extensively compared with conventional measurement of the P–E hysteresis loop. Tip bias was continuously changed from ⫺10 V to ⫹10 V, and its topographic deformation of typical BLCT film library is shown as a curve of butterfly shape as shown in Fig. 6B. Fig. 6C shows the relationship between topographic deformation (H ⫽ H1 ⫹ H2) and remanent polarization (2Pr) in BLCT array as a function of La content (x). H1 and H2 are lift heights with a tip bias, ⫺10 V and ⫹10 V, respectively. Piezoelectricity and ferroelectricity show the similar tendency as a function of La content. Maximum H was 0.9 nm at x ⫽ 0.2, indicating that relative ranking of Pr can be estimated rapidly on the basis of topographic deformation. This method can be used reliably for measuring Pr of discrete or continuous high-density ferroelectric combinatorial array indirectly. In conclusion, we fabricated a BLCT thin-film library by using LSMCD equipped with an automated shutter, and mapped the structure and remanent polarization as a function of La content. This library-synthetic method allows us to control the chemical composition accurately (ppm level in the library). Structural analysis by microbeam XRD clearly shows that (117) crystalline growth has an important role in improving the remanent polarization. Newly developed scanning piezoelectric deformation measurement is found out to be reliable for determining the relative ranking of Pr value rapidly. Materials and Methods Preparation of Precursor Solutions. For the preparation of BLT

Fig. 3. The characterization of P–E characteristic using a ferroelectric tester (RT66A). The plot of 2Pr of BLCT thin-film samples on the library as a function of La content, x. Pt top electrodes of 200 ␮m diameter were deposited on the thin-film library by using RF-sputtering system to fabricate the capacitor structure for measuring the remanent polarization. The thickness of top electrode was ⬇100 nm.

Kim et al.

(Bi/La ⫽ 3.75/0.25, Bi 20% excess) solution, 2.465 g of Bi(NO3)3䡠6H2O and 0.122 g of La(NO3)3䡠6H2O were dissolved in 10 ml of 2-methoxyethanol (Solution 1). One milliliter of Ti isopropoxide was dissolved in 9 ml of 2-methoxyethanol (Solution 2). Solution 1 and 2 were mixed. After the solution mixing, 10 ml of 2-methoxyethanol and 1 ml of 2-ethylhexanoic acid were added and stirred for 30 min. In preparing BCT (Bi/Ce ⫽ 3.75/0.25, Bi 20% excess), Ce nitrate was added instead of La nitrate, and the other conditions were the same as those of preparing BLT solution. Preparation of BLCT Thin Film Library by LSMCD Technique. Schematic diagrams of LSMCD are shown in SI Fig. 8A. We generated the mists of BLT and BCT metal precursor PNAS 兩 January 23, 2007 兩 vol. 104 兩 no. 4 兩 1137

CHEMISTRY

Fig. 4. Fatigue endurance behavior of BLCT thin film samples (x ⫽ 0.014, 0.172, 0.194). The frequency was 1 MHz, and applied voltage was ⫾ 5 V. Only these compositions were tested because of long test time.

Fig. 5. Confirmation of the reliability of combinatorial experimentation by using LSMCD process. (A) The comparison of 2Pr of BLCT samples deposited by combinatorial experimentation with that of BLCT samples deposited individually (one sample at one time). (B) P–E hysteresis loop of BLCT samples chosen among the BLCT library. (C) P–E hysteresis loop of BLCT fabricated individually at the same condition.

solution by using an ultrasonic nebulizer. First, BCT mist was transported to the top of the computer-operated shutter by an Ar flow of 500 sccm. Below the shutter, shadow mask containing 1138 兩 www.pnas.org兾cgi兾doi兾10.1073兾pnas.0610146104

Fig. 6. Mapping of topographic deformation using the SPM. (A) The schematic drawing of measuring topographic deformation by using SPM. (B) The topographic deformation of samples detected by SPM. Detection time was 1 s per point and DC voltage was changed continuously from ⫺10 V and 10 V. A 2 ␮m ⫻ 2 ␮m region of each library was scanned with a gold-coated AFM tip. (C) The plot of topographical deformation (H ⫽ H1 ⫹ H2) and 2Pr of BLCT thin-film samples on the library as a function of La content.

8 ⫻ 8 holes of 3 mm in diameter was intimately contacted on top of the wafer. BLT mist was deposited on the wafer in a stepwise fashion (deposition time: 20 sec/step) by moving shutter in x-direction. Linear variation of deposition time generated a linear concentration gradient of BLT in x-direction. Then, BLT Kim et al.

Characterization of BLCT Thin Film Library Prepared by LSMCD Technique. To analyze the new crystallographic phase and preferred

orientation of as-prepared thin-film library, we can use the microbeam XRD (D8 DISCOVER with GADDS for combinatorial screening by Bruker-AXS) with small beam size (⬍500 ␮m) and xyz manipulator controlled by computer software. The angle between detector and substrate was 20° and the angle between x-ray gun and substrate was 15°. The XRD mapping was carried out in ␪-2␪ scans under Cu K␣ radiation (␭ ⫽ 1.5405 Å) operated at 40 mA and 40 kV. 1. Jandeleit B, Schaefer DJ, Powers TS, Turner HW, Weinberg WH (1999) Angew Chem Int Ed 38:2494–2532. 2. Yoo YK, Xiang XD (2002) J Phys Condens Matter 14:R49–R78. 3. Xiang X-D, Sun A, Briceno G., Lou Y, Wang KA, Chang H, Wallace-Freeman WG, Chen SW, Schultz PG (1995) Science 268:1738–1740. 4. Briceno G, Chang H, Sun XD, Schultz PG, Xiang XD (1995) Science 270:273–275. 5. Wang JS, Yoo YK, Gao G, Takeuchi I, Sun X-D, Chang H, Xiang X-D, Schultz PG (1998) Science 279:1712–1714. 6. Danielson E, Golden JH, McFarland EW, Reaves CM, Weinberg WH, Wu XD (1997) Nature 389:944–948. 7. Koida T, Komiyama D, Koinuma H, Ohtani M, Lippmaa M, Kawasaki M (2002) Appl Phys Lett 80:565–567. 8. Amis EJ (2004) Nat Mater 3:83–85. 9. Koinuma H, Takeuchi I (2004) Nat Mater 3:429–438. 10. Woo SI, Kim KW, Cho HY, Oh KS, Jeon MK, Kim TS, Tarte NH, Mahmood A (2005) QSAR Comb Sci 24:138–154. 11. Takeuchi I, Chang K, Sharma RP, Bendersky LA, Chang H, Xiang X-D, Stach EA, Song C-Y (2001) J Appl Phys 90:2474–2478. 12. Takeuchi I, Famodu OO, Read JC, Aronova MA, Chang K-S, Craciunescu C, Lofland SE, Wuttig M, Wellstood FC, Knauss L, et al. (2003) Nat Mater 2:180–184. 13. Smith RC, Hoilien N, Roberts J, Campbell SA, Gladfelter WL (2002) Chem Mater 14:474–476. 14. Luo ZL, Geng B, Bao J, Gao C (2006) J Comb Chem 7:942–946. 15. Wang J, Evans JRG (2005) J Mater Res 20:2733–2740. 16. Okamura S, Takeuchi R, Shiosaki T (2002) Jpn J Appl Phys 41: 6714–6717. 17. Li GZ, Yu M, Wang ZL, Lin J, Wang RS, Fang J (2006) J Nanosci Nanotechnol 6:1416–1422. 18. Hoffman MM (1995) Integr Ferroelectr 10:39–53. 19. Chung HJ, Kim JH, Moon WS, Park SB, Hwang CS, Lee MY, Woo SI (1997) Integr Ferroelectr 12:185–197.

Kim et al.

The chemical composition of prepared thin-film library was obtained by wavelength dispersive spectroscopy (Microspec 3-PC). Surface morphology of thin-film library was analyzed by SEM (Philips 533M). To measure the electrical properties, Pt electrodes with 200 ␮m diameter were deposited on the thin-film library by using the RF-magnetron sputtering system. P–E hysteresis was measured by using a Radiant Technology RT66A ferroelectric tester at 10 V. The leakage current characteristics were measured by using programmable Kethley 617 electrometer with the condition of 0.05 V of step voltage and 0.1-sec delay time. A scanning probe microscope (SPA 400 SPI 3800N, Seiko Instrument Korea, Seoul, South Korea) was used to measure or scan the topographical change and surface roughness of thin film array. When measuring the topographical change, the detection time was 1 sec per point, and voltage was changed from ⫺10 V and 10 V. The scan region was 2 ␮m ⫻ 2 ␮m, and a gold-coated AFM tip (SI-DF3-A, C ⫽ 2.3 N/m) was used as a movable top electrode. The measurement was carried out ⬎20 times per point, and the induced topographical change was averaged. This research was funded by Center for Ultramicrochemical Process Systems (CUPS) sponsored by Korea Science and Engineering Foundation. 20. Woo SI (2001) Materials 12:4–13. 21. Chung HJ, Kim JH, Woo SI (2001) Chem Mater 13:1441–1443. 22. Kim KW, Woo SI, Choi KH, Han KS, Park YJ (2003) Solid State Ionics 159:25–34. 23. Kim KW, Lee SW, Han KS, Chung HJ, Woo SI (2003) Electrochim Acta 48:4223–4231. 24. Araujo CA, Cuchiaro JD, McMillan LD, Scott MC, Scott JF (1995) Nature 374:627–629. 25. Park BH, Kang BS, Bu SD, Noh NW, Jo W (1999) Nature 401:682–684. 26. Scott JF, Araujo CA (1989) Science 246:1400–1405. 27. Vijay DP, Desu SB (1993) J Electrochem Soc 140:2640–2645. 28. Kojima T, Sakai T, Watanabe T, Funakubo H, Saito K, Osada M (2002) Appl Phys Lett 80:2746–2748. 29. Chon U, Kim K-B, Jang HM, Yi G-C (2001) Appl Phys Lett 79:3137–3139. 30. Watanabe T, Funakubo H, Osada M, Noguchi Y, Miyayama M (2002) Appl Phys Lett 80:100–102. 31. Uchida H, Yoshikawa H, Okada I, Matsuda H, Iijima T, Watanabe T, Kojima T, Funakubo H (2002) Appl Phys Lett 81:2229–2231. 32. Kim KT, Kim CI, Kang DH, Shim IW (2002) Thin Solid Films 422:230–234. 33. Jeon MK, Kim Y-I, Sohn JM, Woo SI (2004) J Phys D 37:2588–2592. 34. Jeon MK, Kim Y-I, Nahm S-H, Woo SI (2004) J Kor Phys Soc 45:1240 – 1243. 35. Jeon MK, Kim Y-I, Nahm S-H, Woo SI (2005) J Phys Chem B 109:968–972. 36. Xu Y (1991) in Ferroelectric Materials and Their Applications (Elsevier, Los Angeles), pp 308–320. 37. Lee HN, Hesse D, Zakharov N, Gosele U (2002) Science 296:2006–2009. 38. Ahn CH, Tybell T, Antognazza L, Char K, Hammond RH, Beasley MR, Fischer O, Triscone J-M (1997) Science 276:1100–1103. 39. Paruch P, Tybell T, Triscone J-M (2001) Appl Phys Lett 79:530–532. 40. Hong JI, Song HW, Hong SB, Shin HJ, No KS (2002) J Appl Phys 92:7434–7441. 41. Ahn CH, Rabe KM, Triscone JM (2004) Science 303:488–491.

PNAS 兩 January 23, 2007 兩 vol. 104 兩 no. 4 兩 1139

CHEMISTRY

mist was deposited in the same manner except for the movement of shutter in the opposite direction. After this procedure, the substrate holder was rotated by 90°; BLT and BCT mists were deposited in the same manner except for the deposition time (13 sec/step) as explained in SI Fig. 8B. The above methods complete a single cycle of deposition. We were then able to generate the liquid phase library of same volume but with different Ce/La ratio. After 3 cycles, we prepared the BLCT thin-film library with 200 nm thickness after thermal treatments at 230°C for 2 min, at 400°C for 10 min in air, and at 700°C for 1 h under O2 atmosphere.