Study on coalescent properties of ZnO ... - Semantic Scholar

Report 6 Downloads 59 Views
Microelectronics Journal 37 (2006) 722–727 www.elsevier.com/locate/mejo

Study on coalescent properties of ZnO nanoclusters using molecular dynamics simulation and experiment Te-Hua Fang a, Win-Jin Chang b,*, Jun-Wei Chiu c a

Institute of Mechanical and Electromechanical Engineering, National Formosa University, Yunlin 632, Taiwan, ROC b Department of Mechanical Engineering, Kun-Shan University, Tainan 710, Taiwan, ROC c Department of Mechanical Engineering, Southern Taiwan University of Technology, Tainan 710, Taiwan, ROC Received 1 September 2005; accepted 18 December 2005 Available online 21 February 2006

Abstract The coalescent properties of ZnO clusters were studied through experiment and molecular dynamics simulation in combination with the tightbinding potential and ZnO potential. The results from the simulation show that the linearly relationship between the melting temperature and the function of atom numbers of NK1/3 was obtainable. Extrapolating the result yield at a melting point of 2130 K for N/N (i.e. NK1/3/0) was slightly lower than the bulk value of 2248 K. In addition, the neck diameter of two ZnO clusters was a function of temperature during coalescence. The contact length was influenced by the coalescence temperature and time, when a cluster was simulated being deposited onto a substrate. The experimental results showed that the grain size increased when the coalescence temperature and sintering time were increased. q 2006 Elsevier Ltd. All rights reserved. Keywords: Molecular dynamics; Tight-binding potential; Coalescent properties; ZnO clusters

1. Introduction Nanoclusters and its coalescence properties are currently attracting lots of attention due to their unique characteristics [1–5]. In the past years, many researchers studied the coalescence properties of metal clusters using molecular dynamic simulation. Lewis et al. [6] developed a moleculardynamics simulation using the embedded-atom method to study the melting, freezing, and coalescence of gold nanoclusters. Qi et al. [7] investigated the melting and freezing of Ni nanoclusters using molecular dynamics with the manybody force field. Zhao et al. [8] studied the coalescence of three identical size nanoclusters through molecular dynamics simulation with an analytic embedded-atom potential. Shim et al. [9] studied the thermal stability of gold nanoparticles, using molecular dynamics simulation in combination with the modified embedded-atom-method potential. Recently, Bilalbegovic [10] studied aggregates of gold clusters by molecular dynamics simulation. Ding et al. [11] utilized molecular dynamics simulation to show that the

* Corresponding author. E-mail address: [email protected] (W.-J. Chang).

0026-2692/$ - see front matter q 2006 Elsevier Ltd. All rights reserved. doi:10.1016/j.mejo.2005.12.007

coalescence of iron nanoclusters occurs at temperatures lower than the cluster melting point, and that the difference between coalescence and melting temperatures increased as the cluster size decreased. Arcidiacono et al. [12] investigated the coalescence process of two gold nanoparticles using a host of initial temperatures and measuring the starting radii in a vacuum with the help of molecular dynamics simulation. Zinc oxide (ZnO) has a wide band gap and good piezoelectric properties. It will have an important application in the field of nanodevices [13–16]. However, very little literature studies coalescent properties of ZnO clusters. In this paper, the effect of the temperature and the structure on the coalescent properties of ZnO nanoclusters was studied using molecular dynamics simulation with the tight-binding potential, ZnO potential and an experimental method. 2. Simulation procedure The interaction forces between ZnO nanocluster zinc atoms, oxygen atoms and zinc atom to oxygen atom were complicated. The many-body tight-binding potential [17] is adopted to simulate the interatomic energy between zinc atoms and expressed as: X Ec ZK ðEbi C Eri Þ (1) i

T.-H. Fang et al. / Microelectronics Journal 37 (2006) 722–727 Table 1 Interaction parameters between zinc atoms for tight-binding potential Parameter

Zn–Zn

Unit

r0 x q A p

2.653 8.9!10K1 4.602 1.477!10K1 9.689

10K1 eV Dimensionless eV Dimensionless

atomic force microscopy (AFM) operated in a noncontact mode and scanning electron microscope (SEM). Coalescent ZnO clusters are insulators and must become conductive to yield scanning electron microscope images through coating a layer of Au. Since the electron beam also generates X-rays emission characteristics of the elements present, energy dispersive spectroscopy (EDS) of the X-rays provides a means of elemental identification. EDS can provide rapid qualitative or with adequate standards, quantitative analysis of the elemental composition with a sampling depth of 1–2 mm.

Table 2 Potential parameters between oxygen atoms Parameter

O–O

Unit

Aij Bij Cij Dij

1.52976!10K22 4.56288 5.127!10K25 0

10K1 nmg/fs2 10 nmK1 10K8 nm8 g/fs2 10K10 nm10 g/fs2

4. Results and discussion

Table 3 Potential parameters between zinc and oxygen atoms Parameter

Zn–O

Unit

Aij Bij Cij Dij

8.4868!10K24 2.79252 0 0

10K1 nmg/fs2 10 nmK1 10K8 nm8 g/fs2 10K10 nm10 g/fs2

where Ebi represents the band-structure energy of the i atom and Eri denotes the repulsive energy of the i atom and they are expressed as:   )1=2 ( X K2q rrijK1 0 Ebi ZK x2 e (2) js1

Eri Z

X

 Kp

Ae

rij r0K1

723



Fig. 1 depicts the dependence of the melting temperature on the ZnO atoms cluster size. The dashed line reflects the linearly relationship between the melting temperature and a function of NK13. The melting point decreased linearly due to curvature effects [11] as the number of atoms decreased. The result was in agreement with previous studies [11,19]. The melting point of 400, 800, 1200, 2800 and 6400 atoms were 1345, 1460, 1590, 1720 and 1815 K, respectively. Extrapolating the result for N/N (i.e. NK1/3/0) yields a predicted value for a cluster consisting of very large N and a melting point of 2130 K, which was slightly lower than the bulk value of 2248 K. The discrepancy might be explained because the cluster calculations have a free surface. Fig. 2 illustrates the relationship between the radial distribution function and the radial separation for a ZnO cluster of 1200 atoms at various temperatures. According to the statements of Fig. 1, it is a solid state for a cluster with 1200 atoms at 700 and 1100 K while a liquid cluster exists at 1700 K. For the cluster in the solid state, it can be seen in the figure that the radial distribution function has higher peaks and deeper valleys than those in the liquid state due to the density effect.

(3)

j

where rij is the distance between zinc atoms i and j; r0 is the first-neighbors distance in the lattice and p, q, x and A are constants which they are listed in Table 1. An ionic potential is used to perform the pair interaction sums among oxygen and zinc atoms and it is expressed as [18] F Z Aij expðKBij rÞK

Cij Dij e2 K 8 C sgn 6 r r r

(4)

where Aij, Bij, Cij and Dij are the potential parameters and listed in Tables 2 and 3 for various atoms and sgnZG1 depending on whether the ions are of similar or opposite charge. 3. Experimental details The ZnO nanoparticles were about 50 nm in diameter and were coalesced in a vacuum chamber. The chamber temperature could be controlled at better than 0.5 K and was adjusted from 700 to 1300 K. The surface morphology of the coalescent clusters after the thermal sintering process was observed with

Fig. 1. The melting point of a ZnO cluster consists of different numbers of atoms.

724

T.-H. Fang et al. / Microelectronics Journal 37 (2006) 722–727

Fig. 2. The relationship between the radial distribution function and the radial separation for a ZnO cluster of 1200 atoms at various temperatures.

Fig. 4. The neck diameter of two ZnO clusters consist of same 1200 atoms is function of temperature during coalescence.

Fig. 3 shows the snapshots of ZnO clusters consisting of 1200 and 400 atoms at a temperature of 1500 K during coalescence. As expected, the two clusters coalesced and gradually became similar to a sphere as the time was increased. Fig. 4 depicts the neck diameter of two ZnO clusters consisting of the same 1200 atoms and a function of temperature during coalescence. The neck diameter increased as the time was increased and finally it reached the diameter of the coalescence sphere. The increasing rate of the neck diameter increased while increasing the temperature of coalescence due to the higher energy. Fig. 5 illustrates the initial configurations of different ZnO clusters consisting of 400 atoms. In Fig. 6, the variation of the

average center distance of ZnO clusters consisting of 400 atoms at different time steps and a temperature of 1000 K. The average center distance decreased as the time was increased during coalescence because they gradually formed a sphere. To study the behavior of the coalescence of a spherical cluster deposited on a substrate, the snapshots of a ZnO cluster consisting of 1200 atoms deposited on a ZnO substrate consisting of 3200 atoms at a temperature of 2000 K during coalescence is shown in Fig. 7.

Fig. 3. Snapshots of ZnO clusters consist of 1200 and 400 atoms at temperature of 1500 K during coalescence.

Fig. 5. Initial configurations of different ZnO clusters consist of 400 atoms.

T.-H. Fang et al. / Microelectronics Journal 37 (2006) 722–727

725

Fig. 6. Variation of the average center distance for different ZnO clusters consist of 400 atoms with time step at temperature of 1000 K.

Fig. 8. Variation of the contact length for a ZnO cluster of 1200 atoms deposited on a ZnO substrate of 3200 atoms with time step at various temperatures.

The spherical cluster was gradually deposited onto the substrate as the time was increased. The contact length between the cluster and the substrate at various coalescence temperatures is illustrated in Fig. 8. In the case of 2000 and 3000 K, the atoms are in a liquid state and the contact length quickly reached a constant value, although, a very long time was require at 1000 K since they were in a solid state. The behaviors of coalescence for ZnO nanoparticles were also studied by using an experimental method. Fig. 9 illustrates EDS of ZnO nanoparticle sintered at a temperature of 1100 8C

at different coalescence times. Au is shown in the figure since ZnO was an insulator and must become a conductor for yielding SEM images through coating a layer of Au. The grain size obtained from the experiment increased while increasing the coalescence time and is shown in Fig. 10. Fig. 11(a) and (b) are the AFM image and SEM image, respectively, of the coalescence of ZnO clusters at 1100 8C after 8 and 10 h. It can be observed that a larger grain size was yielded after 10 h than at 8 h. In addition, the grain size also rapidly increased while increasing the coalescence temperature and is shown in Fig. 12. The color of ZnO samples from the experimental observation, not shown here, changed from white

Fig. 7. Snapshots of ZnO cluster of 1200 atoms deposited on ZnO substrate of 3200 atoms at temperature of 2000 K during coalescence.

Fig. 9. EDS of ZnO atoms at temperature of 1100 8C for different coalescence times.

726

T.-H. Fang et al. / Microelectronics Journal 37 (2006) 722–727

Fig. 12. Dependence of grain size on coalescence temperature after 4 h. Fig. 10. Dependence of grain size on coalescence time at temperature of 1100 8C.

to green to brown when the coalescence temperature was increased from 700 to 1300 8C.

5. Conclusions In this paper, the effect of the temperature and the structure on the coalescent properties of ZnO clusters was studied using molecular dynamics simulation and experimental method. The following results were obtained:

(1) The melting point of the ZnO clusters decreased linearly as the simulated atom increased due to a function of NK1/3 of the curvature effect. (2) The neck diameter rate of increase of two simulated ZnO clusters increased as the temperature during coalescence was increased. (3) The contact length between a simulated ZnO cluster and ZnO substrate increased as the coalescence temperature and time were increased. (4) The grain size obtained from the experiment increased when increasing the coalescence time. (5) The grain size rapidly increased while increasing the coalescence temperature.

Acknowledgments This work was partially supported by the National Science Council of Taiwan, under Grant No. NSC 94-2218-E-150-045, NSC 94-2212-E-168-004. References

Fig. 11. (a) AFM image and (b) SEM image of Coalescence of ZnO clusters at 1100 8C after 8 and 10 h, respectively.

[1] K.F. Peters, J.B. Cohen, Y.W. Chung, Phys. Rev. B 57 (1998) 13430. [2] Z.R. Dai, S. Sun, Z.L. Wang, Nano Lett. 1 (2001) 443. [3] K.E.J. Lehtinen, M.R. Zachariah, Phys. Rev. B 63 (2001) 205402. [4] V. Yadha, J.J. Helble, Aerosol Sci. 35 (2004) 665. [5] A. Kara, T.S. Rahman, Surf. Sci. Rep. 56 (2005) 159. [6] L.J. Lewis, P. Jensen, J.L. Barrat, Phys. Rev. B 56 (1997) 2248. [7] Y. Qi, T. C ¸ agin, W.L. Johnson, W.A. Goddard III, J. Chem. Phys. 115 (2001) 385. [8] S.J. Zhao, S.Q. Wang, Z.Q. Yang, H.Q. Ye, J. Phys.: Condens. Matter 13 (2001) 8061. [9] J.H. Shim, B.J. Lee, Y.W. Cho, Surf. Sci. 512 (2002) 262. [10] G. Bilalbegovic, Comput. Mater. Sci. 31 (2004) 181. [11] F. Ding, A. Rose´n, K. Bolton, Phys. Rev. B 70 (2004) 075416.

T.-H. Fang et al. / Microelectronics Journal 37 (2006) 722–727 [12] S. Arcidiacono, N.R. Bieri, D. Poulikakos, C.P. Grigoropoulos, Int. J. Multiphase Flow 30 (2004) 979. [13] N. Raghu, T.R.N. Kutty, Appl. Phys. Lett. 60 (1992) 100. [14] P. Mitra, A. Chatterjee, H. Maiti, Mater. Lett. 35 (1998) 33. [15] N. Kollias, Arch. Dermatol. 135 (1999) 209.

[16] [17] [18] [19]

727

D.R. Clarke, J. Am. Ceram. Soc. 82 (1999) 485. F. Cleri, V. Rosato, Phys. Rev. B 48 (1993) 22. W. Wunderlich, Phys. Status Solidi 170 (1998) 99. M. Dippel, A. Maier, V. Gimple, H. Wider, W.E. Evenson, R.L. Rasera, G. Schatz, Phys. Rev. Lett. 87 (2001) 095505.