About this Journal Submit a Manuscript Table of Contents
ISRN Materials Science
Volume 2013 (2013), Article ID 710798, 7 pages
http://dx.doi.org/10.1155/2013/710798
Research Article

Influence of Oxygen on Zinc Oxide Films Fabricated by Ion-Beam Sputter Deposition

1Department of Physics, Fu Jen Catholic University, Hsinchuang, New Taipei City 24205, Taiwan
2Graduate Institute of Applied Science and Engineering, Fu Jen Catholic University, Hsinchuang, New Taipei City 24205, Taiwan

Received 25 June 2013; Accepted 30 July 2013

Academic Editors: Z. Jiang and H. Saxén

Copyright © 2013 Jin-Cherng Hsu and Yueh-Sheng Chiang. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

In this study, zinc oxide films were deposited by an ion-beam sputter deposition in various oxygen partial pressures at room temperature. The films changed the structures from amorphous to polycrystalline with increasing the oxygen partial pressure (). The optimal was found at Torr because the film prepared at the oxygen partial pressure had the lowest resistivity and the highest transparence in the visible light region. The lowest resistivity results from a great number of oxygen vacancy sites formed on the polycrystalline surface as exposed to the atmosphere. Moreover, the film has the highest XRD peak intensity, smallest FWHM diffraction peak, smallest -spacing, and smallest biaxial stress.

1. Introduction

ZnO has a hexagonal wurtzite structure with a wide direct band gap of about 3.4 eV. It is naturally an n-type semiconductor due to its deviation from stoichiometry as the presence of oxygen vacancies or zinc interstitials [1]. It is important to fabricate an optimal ZnO film for development of optoelectronic devices due to its high transparence in the visible light region and its low resistivity [2]. In the previous studies of Ti, Zr, Si, and Al metallic oxide films prepared by ion-beam sputter deposition (IBSD) [3], the oxygen partial pressure () can affect their deposition rates, optical properties, and surface morphologies. Depositing at each optimal , the metallic oxide films have higher transparence, higher refractive indices, lower extinction coefficients, and lower surface roughness. They are all amorphous structures and nonconductive films [4]. In this study, an optimal for zinc metallic target was also found. However, the as-deposited zinc oxide film is polycrystalline structure and has conductive property due to free electron carrier, which induces the optical band gap shift by the Burstein-Moss (BM) effect [5] and band-gap-narrowing (BGN) effect [6]. The BM effect is the blue-shift phenomenon of the band edge where the Fermi level merges into the conduction band with the increase of the carrier concentration. Yet, the BGN effect is a red-shift phenomenon of the band edge because the electron-electron repulsive interaction and the localization of the electron wave function is weakened by the screening of the potential in the high carrier concentration [7]. Apparently, the band gap shift is affected by the combination of the BM and BGN effects as increasing carrier concentrations. Sakai et al. have also investigated that the carrier concentration strongly affects the optical band gap resulting from the oxygen vacancies controlled by the oxygen flow rate during the growth of the polycrystalline ZnO film by a reactive plasma deposition [7]. In this study, the optical band gaps and the ZnO structures were determined by controlling the various oxygen partial pressures during the IBSD at room temperature.

In our previous study of IBSD, the adsorbate, water, and hydroxyl () groups are formed on the surface of the as-deposited ZnO films as exposed to the atmosphere. The influences of the adsorbate on the optical property have been illustrated [8]. Furthermore, the electric properties affected by the adsorbate were also investigated in this study.

2. Experiment

Figure 1 shows schematic drawing of the ion beam sputtering system for this study. The system with a Kaufmann type ion source manufactured by Veeco Inc. was equipped with 3 cm diameter molybdenum grids. A target mounted on a water-cooled copper block was pure zinc metal (99.99%) slab with a size of 12 × 10 × 1 cm. Silicon wafer substrate for XRD measurement and one inch diameter B270 glass for optical measurement were cleaned with ethanol in an ultrasonic bath and blow-dried with dry nitrogen gas and then were mounted on a rotatable holder. The chamber was evacuated to a base pressure of 5 × 10−6 Torr by a cryo-pump before deposition. Ar gas was fed into ion source as a working gas at flow rate of ~10 sccm. The ion-beam voltage and ion-beam current were kept at 1000 V and 20 mA, respectively. An incident angle between the ion beam and the surface normal of the zinc target was 45°. The target was presputtered by the Ar ion beam for half an hour to clear the surface in the metallic mode before each deposition. Then, oxygen gas was fed into chamber by mass flow meter during deposition, where was changed from 2 × 10−5 to 2 × 10−4 Torr. The as-deposited sample was labeled by ZnO- as the sample was prepared at  Torr. The total pressure was regulated at about 3 × 10−4 Torr by feeding Ar gas into ion source during deposition.

710798.fig.001
Figure 1: Schematic drawing of the ion-beam sputtering system.

The total transmittance spectrum () of the substrate was measured by Varian Cary-5E spectrometer before mounting. The B270 glass substrate has two surfaces. Each surface has the same transmittance (). An incident light from the left side of the substrate is repeatedly reflected between the two surfaces schematically shown in Figure 2. Then, the transmittance , the sum of the incoherent irradiances at the right side of substrate, can be deduced below [9]: where is the reflection of the single surface and equals () due to its free absorption. Then, can be represented by

710798.fig.002
Figure 2: The two relative transmittances, and , are schematically shown by the sums of the incoherent irradiances in the right side of the B270 substrate, where and are the transmissions of the uncoated and coated surfaces of the substrate and   and are the reflectance of the two surfaces, respectively.

If the as-deposited ZnO is a low optical absorption material in the visible region, the transmittance of the film deposited on the one of substrate’s surface is ; then In (3), is the transmittance measured by Cary 5E spectrometer after deposition and is also the sum of the incoherent irradiances in the right side of the substrate including the repeatedly reflected light schematically shown in Figure 2. It can be similarly deducted from (1). Minus values of in (3) reduce the irradiances interference from the other un-coated substrate’s surface. Measuring the transmittances of and the and calculating the transmittances by the above three formulas, the optical absorption () of the ZnO film is more accurately evaluated by the equation of Beer-Lambert’s law [10], , where small capital is the film’s thickness.

The surface morphologies of the films were investigated by an atomic force microscope, Dimension 3100 made by Digital Instruments. The charge carrier concentrations, sheet resistivities, and the mobilities of the ZnO films were measured by the van der Pauw method with an HL55WIN Hall System. The crystalline structures were characterized by using a Cu line () in X-ray diffractometer made by Rigaku Multiflex. An M-2000 spectroscopic ellipsometer made by J. A. Woollam Co. was used to measure the thicknesses of the thin ZnO films listed in Table 1.

tab1
Table 1: The resistivities and thickness of the films deposited at various oxygen partial pressures.

3. Results and Discussion

Figure 3 shows that oxygen partial pressures affect the deposition rates of ZnO films. The rate begins to decrease at  Torr as zinc suboxide forming on the Zn metallic target surface [4]. As soon as  Torr, oxidization layer of zinc metal formed on the target surface decreases the sputtering yield. A sharp drop of deposition rate is about three times less than the rare at the low . Castellano [11] also has explained this phenomenon by ideal gas theory to predict the sharp drop of deposition rate in narrow range of at which a film forms from metal to oxide. The basic mechanism of this transition is the trapping of oxygen molecules on the target surface [12]. Figure 3 shows that the deposition rate similarly decreases and the surface roughness of ZnO film decreases with increasing because each deposited particle averages higher momentum transport rate from the energetic argon ion beam at the lower deposition rate. Then, the deposited particles with the higher surface mobility aggregate smoothly on the substrate.

710798.fig.003
Figure 3: Deposition rates and surface roughness of the ZnO films versus oxygen partial pressure.
3.1. Optical Properties

Figure 4 shows the optical transmittance spectra of the ZnO films between 250 and 750 nm. The transmittances of ZnO-2 to -8 films deposited at  Torr are lower than 80% in visible region due to the heavier optical absorption of the suboxide film. The transmittances of ZnO-10 to -20 deposited at lower deposition rate at  Torr are larger than 80% as the transparent films approach to stoichiometric [13, 14]. However, the transmittances decrease sharply at absorption edges of about 350 nm wavelength. When increasing , the red-shift phenomena of the absorption edges obviously happen at the ZnO-10 to -20 samples. For describing precisely, the optical band gap () of ZnO film is determined by the equation of for the direct transition, where , , and are a constant, the incident photon energy, and the absorption coefficient, respectively. The absorption coefficient is calculated by Beer-Lambert’s law with the equation of [14]. Then, are obtained by Tauc’s plot [15] with extrapolating the linear portion to the photon energy axis as plotted the versus in Figure 5. The absorption edges red-shift from 3.41 to 3.31 eV with increasing .

710798.fig.004
Figure 4: Transmittance of the ZnO films deposited at different oxygen partial pressures.
710798.fig.005
Figure 5: The optical band gaps of the ZnO-10 to -20 samples determined by extrapolating the linear part of the curve to the photon energy axis.

The red-shift phenomenon is also proved by Tan et al. They have fabricated the ZnO films by metal organic chemical vapor deposition at the growth temperatures from 200 to 500°C and found their films transfer from amorphous to polycrystalline structure at the treatment processes [16]. Sun et al. have also studied that the has a close relationship with the amount of oxygen vacancy sites that induce the electron carriers in the polycrystalline ZnO film deposited by laser ablating Zn target in various [17]. In our study, the ZnO-2 to -8 films belong to amorphous structures, and the polycrystalline structures form in the ZnO-10 to -20 films illustrated by XRD as shown in Figure 6. The structures gradually transfer from amorphous to polycrystalline structure with increasing during deposition. The ZnO-10 structure belongs to the mixture of the amorphous and polycrystalline structures because a little XRD peak begins to grow at of Torr. That is indirectly proved by the fact that its of 3.41 eV is larger than that of zinc oxide material having a polycrystalline structure with a direct band gap of about 3.37 eV [18]. Because the polycrystalline structure has oxygen vacancy sites to induce the electron carriers, the amount of the carrier also increases with increasing . Until  Torr, the amount decreases again due to the fact that oxygen atoms occupy their oxygen vacancies during deposition. Illustrating ZnO-10 to -15 samples from Figures 5 and 7, the absorption edges red-shift from 3.41 to 3.32 eV with respect to the increase sheet concentrations from 3.01 × 1015 to 1.82 × 1016 cm−2 with increasing . The red-shift phenomenon is the electron-electron repulsive interaction from the many electron carriers due to the band-gap-narrowing (BGN) effect [7], that is, the optical band gaps red-shift resulting from the increases of electron carrier concentrations in the films. However, the sheet concentration reduces to 6.15 × 1015 cm−2; the optical band gap still slightly decreases to 3.31 eV due to BM effect at  Torr. Lu et al. [19] and Sakai et al. [7] have mentioned that the carrier concentration in the film affects the shift of the optical absorption edge for ZnO films in both the BM effect and the BGN effect as a function of , where described the optical band gap as the minimum energy needed to excite an electron from the valence band to the conduction band. is the energy gap between the conduction and valence band edges in the pure undoped crystal. is the energy blue-shift due to BM effect. is the energy red-shift due to BGN effect. The band gap shift determined by a combination of the two effects is a function of . Thus, there is a critical , 1.5 × 10−4 Torr, to determine the band gap () blue- or red-shift.

710798.fig.006
Figure 6: The (002) peaks of XRD clearly measured from ZnO-10 to -20 films.
710798.fig.007
Figure 7: Mobility, sheet concentration, and sheet resistivities of the ZnO films plotted as a function of oxygen partial pressure.
3.2. Electrical Properties

Figure 7 shows the measurements of Hall effect for the mobility and sheet concentration of the ZnO films. The sheet concentrations increase with increasing until 1.2 × 10−4 Torr, because the electron carriers are induced by the increasing amount of oxygen vacancies and intrinsic defects of ZnO films. The mobilities decrease from 10.4, 4.08 to 2.16 cm2/V·s because of the electron scattering and grain-boundary scattering in polycrystalline film [20]. Yet, the sheet conductivity is proportional to the product of the mobility and carrier concentration of electrons if the conductivity is primarily due to electron. In the absence of the largest amount of sheet concentration and mobility, ZnO-12 film still has the highest sheet conductivity. Because the sheet resistivity is defined as the inverse of the sheet conductivity, ZnO-12 has the least sheet resistivity of 256 Ω/square as shown in Table 1. Figure 8 also shows that the resistivities decrease as increasing until 1.2 × 10−4 Torr and increase again as  Torr, and ZnO-12 has the least resistivity of 2.4 mΩ·cm.

710798.fig.008
Figure 8: Resistivities of the ZnO films plotted as a function of oxygen partial pressure.

Nevertheless, the decrease of the sheet concentration of ZnO-20 deposited at higher  Torr is attributed to improvement of the stoichiometry that may decrease the concentration of the oxygen vacancies and intrinsic defects. Moreover, the mobility increases a little again at the as shown in Figure 7. It is for that the electron scattering dominates the mobilities of films resulting from the nearly invariable gain-boundary scattering on the almost same grain sizes of the films described by XRD in the following section.

3.3. XRD Measurement of the ZnO Films

It seems that the electrical property is not only attributed to the film’s oxidation. ZnO film has a strong tendency for self-organized growth as the film has a strong preferred -axis growth orientation [15]. When the films are about 100 nm thick prepared at  Torr, the films obviously have (002) XRD peaks locating at double diffraction angle () around 34.17° as shown in Figure 6. Table 2 shows the fitted results of the XRD peaks with the intensity, double diffraction angle (), and FWHM (). The diffraction angle () can be represented by -spacing (), between the atomic (002) planes, calculated by Bragg’s diffraction equation, , and the FWHM can be represented by grain size () of the ZnO crystal, by Scherrer’s formula [21], as a function of shown in Figure 9. According to the results of Figures 8 and 9, the higher the crystallinity is the lower the resistivity at  Torr is. The correlation between the most intensity of XRD and the lowest resistivity of 2.4 mΩ·cm for ZnO-12 sample definitely indicates that the crystalline structure plays an important role in the electrical behavior of the ZnO film. Besides, the crystalline ZnO-12 film has the least of 2.623 Å and the largest of 17.5 nm shown in Figure 9 due to the best crystallization at the highest diffraction angle, , shown in Table 2.

tab2
Table 2: The XRD spectra of the ZnO samples, from ZnO-10 to ZnO-20, fitted by Gauss functions and represented by intensities, double diffraction angles, and FHWMs.
710798.fig.009
Figure 9: The XRD intensities, -spacings, and grain sizes of the ZnO-10 to -20 samples versus oxygen partial pressure.

Besides, the biaxial strains of the ZnO films in -axis direction are determined by the previous XRD as shown in Figure 10. The residual stress () based on the biaxial strain model was calculated by the formula,  GPa, where , the lattice constant of ZnO film, was calculated by the following equation [22]: and is the unstrained lattice parameter of 5.205 Å from a ZnO powder sample [23], where , , and are Miller indices of the ZnO (002) crystallinity. ZnO-12 film also has the least biaxial stress of −1.848 GPa, which indicates that the film has a much denser crystalline structure than the others [24]. The narrowest -spacing, largest grain, size and denser structure of ZnO-12 film may have the advantage of reducing the electron conduction path in the polycrystalline structure and decreasing the resistivity.

710798.fig.0010
Figure 10: The -lattice constants and biaxial stresses of the ZnO-10 to -20 samples ploted as a function of oxygen partial pressure.

Furthermore, water molecular is easily absorbed by oxygen vacancy sites on the ZnO film [5, 25]. The electric property may relate to not only the polycrystalline structure but also to the adsorbate on the surface. Schoenes et al. have reported that ZnO film had a larger conductivity resulting from oxygen vacancy sites or Zn interstitial atoms [24, 26]. ZnO-12 has the highest XRD peak, largest grains, and denser structure. Exposing the ZnO-12 film to the atmosphere, the oxygen vacancy and interstitial sites on the crystalline surface absorbed the water and hydroxyl groups. The greater amounts of the absorbents cover the surface of ZnO-12 film as the two subpeaks locating at the binding energies of 532.7 and 531.7 eV have the largest areas in XPS spectra. That has been proved by our previous study [8]. It is deduced that the ZnO crystalline formation and oxygen vacancy sites on the ZnO surface is associated with oxygen partial pressure during deposition, and the electron carriers produced by this crystalline structure benefit lower resistance of ZnO film as exposed to the atmosphere.

4. Conclusion

In this study, optical and electronic properties of ZnO films were controlled by changing oxygen partial pressure in IBSD of zinc metallic target with ion-beam voltage of 1000 V and ion-beam current of 20 mA at room temperature. The optical band gaps of the as-deposited ZnO films are influenced by the two BM and BGN effects. Moreover, comparing the optical and electronic properties, an optimal was found at 1.2 × 10−4 Torr. At the optimal , the ZnO film has the highest optical transmittance and the lowest resistivity of 2.4 mΩ·cm. The low resistivity attributes to water and hydroxyl groups absorbed by a great number of oxygen vacancy sites on the surface of the polycrystalline structure which has the highest XRD peak intensity, the smallest FWHM of diffraction peak, the least -spacing, and the least biaxial stress. This study for finding the optimal oxygen partial pressure of IBSD also benefits the fabrication of aluminum-doped ZnO film [27].

Acknowledgment

The authors acknowledge the sponsorship of the National Science Council of Taiwan under Grant no. NSC102-2221-E-030-011.

References

  1. T. Yamamoto and H. Katayama-Yoshida, “Solution using a codoping method to unipolarity for the fabrication of p-type ZnO,” Japanese Journal of Applied Physics, vol. 38, no. 2, pp. L166–L169, 1999. View at Scopus
  2. E. Budianu, M. Purica, F. Iacomi, C. Baban, P. Prepelita, and E. Manea, “Silicon metal-semiconductor-metal photodetector with zinc oxide transparent conducting electrodes,” Thin Solid Films, vol. 516, no. 7, pp. 1629–1633, 2008. View at Publisher · View at Google Scholar · View at Scopus
  3. C.-C. Lee, D. T. Wei, J.-C. Hsu, and C.-H. Shen, “Influence of oxygen on some oxide films prepared by ion beam sputter deposition,” Thin Solid Films, vol. 290-291, pp. 88–93, 1996. View at Publisher · View at Google Scholar · View at Scopus
  4. C.-C. Lee, J.-C. Hsu, and D.-H. Wong, “Characteristics of some metallic oxides prepared in high vacuum by ion beam sputtering,” Applied Surface Science, vol. 171, no. 1-2, pp. 151–156, 2001. View at Publisher · View at Google Scholar · View at Scopus
  5. E. Burstein, “Anomalous optical absorption limit in InSb,” Physical Review, vol. 93, no. 3, pp. 632–633, 1954. View at Publisher · View at Google Scholar · View at Scopus
  6. A. P. Roth, J. B. Webb, and D. F. Williams, “Band-gap narrowing in heavily defect-doped ZnO,” Physical Review B, vol. 25, no. 12, pp. 7836–7839, 1982. View at Publisher · View at Google Scholar · View at Scopus
  7. K. Sakai, T. Kakeno, T. Ikari et al., “Defect centers and optical absorption edge of degenerated semiconductor ZnO thin films grown by a reactive plasma deposition by means of piezoelectric photothermal spectroscopy,” Journal of Applied Physics, vol. 99, no. 4, Article ID 043508, 7 pages, 2006. View at Publisher · View at Google Scholar · View at Scopus
  8. J.-C. Hsu, Y.-H. Lin, P. W. Wang, and Y.-Y. Chen, “Spectroscopic ellipsometry studies on various zinc oxide films deposited by ion beam sputtering at room temperature,” Applied Optics, vol. 51, no. 9, pp. 1209–1215, 2012. View at Publisher · View at Google Scholar · View at Scopus
  9. H. A. Macleod, Thin-Film Optical Filters, CRC Press, Taylor & Francis, 4th edition, 2010.
  10. W. T. Lim and C. H. Lee, “Highly oriented ZnO thin films deposited on Ru/Si substrates,” Thin Solid Films, vol. 353, no. 1, pp. 12–15, 1999. View at Publisher · View at Google Scholar · View at Scopus
  11. R. N. Castellano, “Reactive ion beam sputtering of thin films of lead, zirconium and titanium,” Thin Solid Films, vol. 46, no. 2, pp. 213–221, 1977. View at Scopus
  12. C.-C. Lee, J.-C. Hsu, D. T. Wei, and J.-H. Lin, “Morphology of dual beam ion sputtered films investigated by atomic force microscopy,” Thin Solid Films, vol. 308-309, no. 1–4, pp. 74–78, 1997. View at Scopus
  13. N. A. Suvorova, I. O. Usov, L. Stan et al., “Structural and optical properties of ZnO thin films by rf magnetron sputtering with rapid thermal annealing,” Applied Physics Letters, vol. 92, no. 14, Article ID 141911, 2008. View at Publisher · View at Google Scholar · View at Scopus
  14. N. Mehan, V. Gupta, K. Sreenivas, and A. Mansingh, “Effect of annealing on refractive indices of radio-frequency magnetron sputtered waveguiding zinc oxide films on glass,” Journal of Applied Physics, vol. 96, no. 6, pp. 3134–3139, 2004. View at Publisher · View at Google Scholar · View at Scopus
  15. J. Tauc, R. Grigorovici, and A. Vancu, “Optical properties and electronic structure of amorphous germanium,” Physica Status Solidi B, vol. 15, no. 2, pp. 627–637, 1966.
  16. S. T. Tan, B. J. Chen, X. W. Sun et al., “Blueshift of optical band gap in ZnO thin films grown by metal-organic chemical-vapor deposition,” Journal of Applied Physics, vol. 98, no. 1, Article ID 013505, 2005. View at Publisher · View at Google Scholar · View at Scopus
  17. L. Sun, W. Cheng, F. Lin, X. Ma, and W. Shi, “Changes of structure and optical energy gap induced by oxygen pressure during the deposition of ZnO films,” Physica B, vol. 381, no. 1-2, pp. 109–112, 2006. View at Publisher · View at Google Scholar · View at Scopus
  18. M.-H. Zhao, Z.-L. Wang, and S. X. Mao, “Piezoelectric characterization individual zinc oxide nanobelt probed by piezoresponse force microscope,” Nano Letters, vol. 4, no. 4, pp. 587–590, 2004. View at Publisher · View at Google Scholar · View at Scopus
  19. J. G. Lu, S. Fujita, T. Kawaharamura et al., “Carrier concentration dependence of band gap shift in n-type ZnO:Al films,” Journal of Applied Physics, vol. 101, no. 8, Article ID 083705, 2007. View at Publisher · View at Google Scholar · View at Scopus
  20. A. F. Mayadas and M. Shatzkes, “Electrical-resistivity model for polycrystalline films: the case of arbitrary reflection at external surfaces,” Physical Review B, vol. 1, no. 4, pp. 1382–1389, 1970. View at Publisher · View at Google Scholar · View at Scopus
  21. B. D. Cullity, Elements of X-Ray Diffractions, Addison-Wesley, Reading, Mass, USA, 1978.
  22. C. S. Barret and T. B. Massalski, Structure of Metals, Pergamon Press, Oxford, UK, 1980.
  23. R. Cebulla, R. Wendt, and K. Ellmer, “Al-doped zinc oxide films deposited by simultaneous rf and dc excitation of a magnetron plasma: relationships between plasma parameters and structural and electrical film properties,” Journal of Applied Physics, vol. 83, no. 2, pp. 1087–1095, 1998. View at Scopus
  24. J. Schoenes, K. Kanazawa, and E. Kay, “Band and hopping conduction in high-resistivity ZnO,” Journal of Applied Physics, vol. 48, no. 6, pp. 2537–2542, 1977. View at Publisher · View at Google Scholar · View at Scopus
  25. A. Önsten, D. Stoltz, P. Palmgren, S. Yu, M. Göthelid, and U. O. Karlsson, “Water adsorption on znO(0001): transition from triangular surface structures to a disordered hydroxyl terminated phase,” Journal of Physical Chemistry C, vol. 114, no. 25, pp. 11157–11161, 2010. View at Publisher · View at Google Scholar · View at Scopus
  26. M. Suchea, S. Christoulakis, N. Katsarakis, T. Kitsopoulos, and G. Kiriakidis, “Comparative study of zinc oxide and aluminum doped zinc oxide transparent thin films grown by direct current magnetron sputtering,” Thin Solid Films, vol. 515, no. 16, pp. 6562–6566, 2007. View at Publisher · View at Google Scholar · View at Scopus
  27. Y.-Y. Chen, J.-C. Hsu, P. W. Wang, Y.-W. Pai, C.-Y. Wu, and Y.-H. Lin, “Dependence of resistivity on structure and composition of AZO films fabricated by ion beam co-sputtering deposition,” Applied Surface Science, vol. 257, no. 8, pp. 3446–3450, 2011. View at Publisher · View at Google Scholar · View at Scopus