Evolution of structure of SiO_{2}nanoparticles upon cooling from the melt
- Nguyen Thi Xuan Huynh^{1}Email author,
- Vo Van Hoang^{2}Email author and
- Hoang Zung^{3}
Received: 20 April 2008
Accepted: 30 October 2008
Published: 30 October 2008
Abstract
Evolution of structure of spherical SiO_{2} nanoparticles upon cooling from the melt has been investigated via molecular-dynamics (MD) simulations under non-periodic boundary conditions (NPBC). We use the pair interatomic potentials which have weak Coulomb interaction and Morse type short-range interaction. The change in structure of SiO_{2} nanoparticles upon cooling process has been studied through the partial radial distribution functions (PRDFs), coordination number and bond-angle distributions at different temperatures. The core and surface structures of nanoparticles have been studied in details. Our results show significant temperature dependence of structure of nanoparticles. Moreover, temperature dependence of concentration of structural defects in nanoparticles upon cooling from the melt toward glassy state has been found and discussed.
PACS Codes: 61.43.Bn; 78.55.Qr; 78.67.Bf
Introduction
Silica nanoparticles have potential applications in many fields including ceramics, chromatography, catalysis and chemical mechanical polishing [1]. In recent years, SiO_{2} nanoparticles have been investigated by means of experimental techniques such as NMR (nuclear magnetic resonance), SAXS (small angle X-ray scattering) [2], light absorption [3], FTIR (Fourier transform infrared) spectra and photoluminescence [4–7], etc. According to the diffraction data [8], liquid and amorphous silica have Z_{si-O}= 4 and Z_{O-Si}= 2, i.e. the main structural element of the network is a slightly distorted SiO_{4} tetrahedron and the adjacent tetrahedra are linked to each other through the shared vertices. This means that Si atoms having 1, 2, 3 or 5-fold coordinations and O atoms have 1 or 3-fold coordinations can be considered as structural defects in SiO_{2} nanoparticles, which are similar to those observed in the bulk [9, 10]. Among them, the most serious defects are structural units with Z_{si-O}= 2, Z_{si-O}= 3 and Z_{si-O}= 1. It was found that the structural defects can play an important role in their structure and properties including photoluminescence [6], catalysis and micro-electronics [11]. In particular, the structural defects lead to the narrowing of band gap and the formation of localized states within band tail, which were experimentally measured from light-absorption study of silica nanoparticles [3]. However, very few theoretical researches have been done for SiO_{2} nanoparticles despite the observation of interesting properties [11–16]. Therefore, it motivates us to carry out the research on liquid and amorphous SiO_{2} nanoparticles by MD simulation via the microstructural analysis. Namely, we investigate the structural evolution of SiO_{2} nanoparticles upon cooling from the melt toward glassy state (or amorphous one).
Calculation
Where q_{ i }and q_{ j }represent the charges of ions i and j, for Si atom q_{ Si }= 1.30e and for O atom q_{0} = -0.65e (e is the elementary charge unit); r denotes the interatomic distance between atoms i and j; the other parameters of the Morse potentials can be found in Refs. [17–19]. MD simulations were done in a spherical SiO_{2} particle with three different sizes of 2 nm, 4 nm and 6 nm corresponding the real density of 2.20 g/cm^{3}for amorphous SiO_{2} which have corresponding numbers of atoms of 276 (92 silicon atoms and 184 oxygen ones), 2214 (738 silicon atoms and 1476 oxygen ones) and 7479 (2493 silicon atoms and 4986 oxygen ones), respectively. We used the Verlet algorithm with the MD time step of 1.60 fs. Each model contained the number of Si and O atoms in accordance to the SiO_{2} stoichiometry. Firstly, N atoms are randomly placed in a sphere of fixed radius and the NPBC model has been relaxed for 5 × 10^{4} MD steps at 7000 K in order to get a good equilibrated liquid model. We use the simple non-slip with non-elastic reflection behavior boundary (i.e. if during the relaxation atoms move out of the spherical boundary they have to be placed back to the surface and it is called NPBC). Then, the temperature of the system was decreased linearly in time as T = T_{0} - αt where α = 4.2956 × 10^{13}K/s is the cooling rate, T_{0} is the initial temperature of 7000 K and t is the cooling time. This cooling process was continued until the temperature of the system was equal to 350 K. It is essential to notice that melting and evaporation points for bulk SiO_{2} are T_{ m }= 1923K and T_{ boil }= 2503K, respectively [1]. This means that the initial configuration at T = 7000K corresponds to the superheated liquid and that only due to using of NPBC it is possible to reach such a condition for SiO_{2} nanoparticles otherwise the evaporation occurs. Note that using potentials have been successfully used in MD simulations of both structure and thermodynamic properties of silica [17–19], and in the investigation of the structure changes in cristobalite and silica glass at high temperatures [18]. These potentials reproduced well the melting temperature of cristobalite and the glass phase transition temperature of silica glass and calculated data were more accurate than those observed in other simulation works in which the traditional interatomic potentials with more strong electrostatic interaction have been used such as BKS or TTAM potentials [20–22]. Although the potentials described above were proposed for the bulk SiO_{2}, they also described well structural features and surface energy of amorphous SiO_{2} nanoparticles compared with those observed for amorphous SiO_{2} nanoclusters by using BKS ones [15]. Melting point of nanoparticles is size dependent in that it strongly reduces with decreasing nanoparticle size [23], and it is out of scope of the present work. In order to investigate the evolution of structure of SiO_{2} nanoparticles upon cooling from the melt, we saved a number of configurations at finite temperatures (i.e. at temperatures ranged from 7000 K to 350 K with the temperature interval of 700 K). and then we relaxed them for 5 × 10^{4} MD steps before calculating static properties. In order to calculate the coordination number and bond-angle distributions in SiO_{2} nanoparticles, we adopted the fixed values R_{Si-Si}= 3.30 Å, R_{Si-O}= 2.10 Å and R_{O-O}= 3.00 Å. Here R corresponds to the position of the minimum after the first peak in PRDFs for the amorphous state at a real density of 2.20 g/cm^{3}. In order to improve statistics, the results have been averaged over four and three independent runs for nanoparticles with the size of 2 nm and 4 nm, respectively. The single run was done for the size of 6 nm due to large number of atoms in the model. In addition, since surface structure plays an important role in the structure and properties of nanoparticles we also focus attention to the surface of SiO_{2} nanoparticles. Therefore, we need a criterion to decide which atoms belong to the surface and which ones belong to the core of nanoparticles. There is no common principle for such choice of surface or core of the amorphous substances. The definition of thickness of the surface in [13] is somewhat arbitrary: all atoms that were within 5 Å of the hull just touched the exterior of the droplet, and were considered to belong to the surface, atoms that had the distance between 5 Å and 8 Å from the hull belong to the transition zone and the remaining atoms belong to the interior. In contrast, no definition of surface was clearly presented for the amorphous Al_{2}O_{3} thin film; they used the top 1 Å or 3 Å layer of the amorphous thin film for surface structural studies [24]. From structural point of view, it can be considered that atoms belong to the surface if they could not have full coordination for all atomic pairs in principle and in contrast, atoms belong to the core if they could have full coordination for all atomic pairs in principle like those located in the bulk. Therefore, in the present work atoms located in the outer shell of SiO_{2} spherical nanoparticle with thickness of 3.30 Å (i.e. the largest radius of the coordination spheres found in the system) belong to the surface of the nanoparticles and the remaining atoms belong to the core.
Results and discussion
Structural characteristics of 4 nm SiO_{2} nanoparticle upon cooling from the melt.
T (K) | r _{ ij } | g _{ ij } | Z _{ ij } | |||||||
---|---|---|---|---|---|---|---|---|---|---|
Si-Si | Si-O | O-O | Si-Si | Si-O | O-O | Si-Si | Si-O | O-Si | O-O | |
7000 | 2.98 | 1.51 | 2.61 | 2.90 | 9.29 | 2.79 | 2.27 | 3.00 | 1.56 | 3.97 |
3500 | 3.02 | 1.52 | 2.58 | 4.52 | 14.93 | 4.08 | 2.73 | 3.59 | 1.80 | 4.83 |
350 | 3.02 | 1.52 | 2.54 | 8.86 | 32.22 | 9.30 | 3.76 | 3.93 | 1.96 | 5.96 |
Exp. [8] | 3.12 | 1.62 | 2.65 | 4.00 | 4.00 | 2.00 | 6.00 |
Percentage of Si and O atoms with corresponding coordination numbers for SiO_{2} nanoparticle with the size of 4 nm upon cooling from the melt to the amorphous state.
T (K) | Z _{Si-O} | Z _{O-Si} | |||
---|---|---|---|---|---|
2 | 3 | 4 | 1 | 2 | |
7000 | 21.68 | 53.25 | 22.13 | 46.25 | 51.33 |
6300 | 16.58 | 55.51 | 25.93 | 43.52 | 53.86 |
5600 | 13.47 | 54.27 | 30.41 | 41.13 | 56.30 |
4900 | 7.90 | 53.93 | 36.50 | 35.46 | 62.12 |
4200 | 4.74 | 50.77 | 42.82 | 30.90 | 67.10 |
3500 | 1.99 | 37.40 | 59.44 | 22.09 | 75.71 |
2800 | 0.72 | 23.98 | 74.75 | 13.87 | 84.70 |
2100 | 0.72 | 11.20 | 87.85 | 6.98 | 92.16 |
1400 | 0.72 | 6.82 | 92.32 | 4.97 | 94.08 |
350 | 0.72 | 6.01 | 93.13 | 4.95 | 93.72 |
Parameters fitted with the Arrhenius law of the probability for the occurrence of structural defects upon cooling from the melt.
Size | Z_{Si-O}= 2 | Z_{Si-O}= 3 | Z_{O-Si}= 1 | |||
---|---|---|---|---|---|---|
A _{ ij } | E_{ ij }(K) | A _{ ij } | E_{ ij }(K) | A _{ ij } | E_{ ij }(K) | |
2 nm | 1.884 | 13052.693 | 1.781 | 4870.674 | 0.954 | 4033.876 |
4 nm | 2.223 | 16215.199 | 2.435 | 6515.887 | 1.162 | 5919.718 |
6 nm | 3.273 | 18964.670 | 2.782 | 7391.125 | 1.303 | 6814.464 |
Data in [9] | 58.600 | 31100.000 | 8.900 | 24760.000 |
It is interesting to discuss about the influence of the choice of interatomic potentials and boundary conditions on the obtained results. First, basing on the results obtained for the bulk SiO_{2} by using different interatomic potentials in literature one can infer that static properties of models are less sensitive dependent on the interatomic potentials used in simulation compared those of the dynamic or thermodynamic ones [9, 10, 29]. Similar results can be suggested for SiO_{2} nanoparticles, however, it is also of interest to expand the research in this direction. Quite other situation for the influence of the boundary on the properties of nanoparticles may be occurred. Indeed, according to the model developed in [30], the size dependence of glass transition temperature of nanosized materials might be related to the boundary conditions. The fact, we obtained SiO_{2} nanoparticles by using two different boundary conditions, i.e. elastic reflection boundary conditions and the non-elastic reflection ones [15]. We found that structural properties of amorphous SiO_{2} nanoparticles obtained by using two different boundary conditions are similar to each other and a significant discrepancy was found for the surface structure. Applying elastic reflection boundary conditions leads to an increase of structural defects at the surface of 2 nm nanoparticle compared with those observed in non-elastic ones. In contrast, no systematic changes have been found for 4 nm nanoparticle. In addition, we also found the boundary effects on the density profile of 4 nm nanoparticle. Such relatively small effects on static properties might lead to dramatic changes in the dynamics of atomic species in nanoparticles [15].
In addition, pressure-induced structural changes in nanoparticles are also of great interest. In particular, we found pressure-induced transition from low density form to high density one of amorphous TiO_{2} nanoparticles [31]. Transition pressure of amorphous TiO_{2} nanoparticles is size dependent in that the smaller nanoparticle size the higher the transition pressure like those commonly found in practice for nanocrystalline substances. Moreover, transition pressure of nanoparticles is higher than that of the bulk counterpart due to the surface effects [31]. On the other hand, so-called layer structure of amorphous TiO_{2} nanoparticles is enhanced with increasing density of a model. The phenomenon oxygen atoms have a tendency to concentrate at the surface of amorphous nanosized substances has been found and it is strongly enhanced by pressurization in addition to the same but weaker tendency of Ti atoms. It causes a strong violation of local charge neutrality at the surface of nanoparticles, and consequently, in order to achieve the local charge neutrality Ti and O atoms also have a tendency to respectively concentrate in the inner shells. It leads to the forming so-called shell structure of compressed nanoparticles (i.e. respectively ordering spherical shells with Ti-rich and O-rich layers). Similar results for amorphous SiO_{2} nanoparticles can be suggested since we found the same tendency caused by compression of initial 4 nm SiO_{2} nanoparticle [32].
Conclusion
i) Calculations show that upon cooling from the melt toward amorphous state, tetrahedral network structure is forming in SiO_{2} nanoparticles. Amorphous SiO_{2} nanoparticles have a slightly distorted tetrahedral network structure if their size is large enough in that Si atoms are mainly surrounded by four O atoms like those observed in the bulk.
ii) The probability for the occurrence of structural defects in SiO_{2} nanoparticles corresponding with Z_{si-O}= 2, Z_{si-O}= 3 and Z_{O-Si}= 1 have been described well by an Arrhenius law at high temperatures. However, prefactors A_{ ij }and activation energies E_{ ij }of the Arrhenius law for our models differ from those observed in the bulk.
iii) During the temperature decreasing, the increase of numerical density in the core of nanoparticle and its decrease in the surface layer are observed.
iv) Moreover, SiO_{ x }stoichiometry in the core and in the surface shell is also different from each other, i.e. while surface has a positive deviation from the SiO_{2} stoichiometry, the core has a negative one over wide range of high temperature region (i.e. in the liquid state). The situation changes slightly at low temperatures. This means that in the high temperature region the core behaves oxygen-deficiency defects due to a negative deviation from the SiO_{2} stoichiometry, while the surface shell has oxygen-excess defects due to a positive deviation from the SiO_{2} one. This phenomenon may cause the appearance of additional defects in the SiO_{2} nanoparticles.
v) In the present work, simulations were performed just at the ambient pressure density of 2.20 g/cm^{3} for amorphous SiO_{2}, it is of interest to carry out the study at different densities ranged from low to high ones. Indeed, the pressure effects on static properties and dynamics of amorphous SiO_{2} have been under intensive investigation [19, 33]. In addition, it is better if one can carry out the quantum MD simulation for SiO_{2} nanoparticles since the quantum effects become more important in the systems at nanoscale.
Declarations
Acknowledgements
This work was accomplished in the Computational Physics Lab of the College of Natural Sciences – National University of HochiMinh City – Vietnam.
Authors’ Affiliations
References
- Iler RK: The chemistry of silica (Wiley, New York, 1979).Google Scholar
- Green DL, Lin JS, Lam Y-F, Hu MZ-C, Schaefer DW, Harris MT: J Colloid Interface Sci. 2003, 266: 346-10.1016/S0021-9797(03)00610-6.View ArticleGoogle Scholar
- Altman IS, Lee D, Chung JD, Song J, Choi M: Phys Rev B. 2001, 63: R161402-10.1103/PhysRevB.63.161402.View ArticleADSGoogle Scholar
- Glinka YD, Lin S-H, Chen Y-T: Phys Rev B. 2002, 66: 035404-10.1103/PhysRevB.66.035404.View ArticleADSGoogle Scholar
- Glinka YD, Lin S-H, Chen Y-T: Appl Phys Lett. 1999, 75: 778-10.1063/1.124510.View ArticleADSGoogle Scholar
- Colder A, Huisken F, Trave E, Ledoux G, Guillois O, Reynaud C, Hofmeister H, Pippel E: Nanotechnology. 2004, 15: L1-10.1088/0957-4484/15/3/L01.View ArticleADSGoogle Scholar
- Glinka YD, Lin S-H, Chen Y-T: Phys Rev B. 2000, 62: 4733-10.1103/PhysRevB.62.4733.View ArticleADSGoogle Scholar
- Mozzi RL, Warren BE: J Appl Cryst. 1969, 2: 164-10.1107/S0021889869006868.View ArticleGoogle Scholar
- Horbach J, Kob W: Phys Rev B. 1999, 60: 3169-10.1103/PhysRevB.60.3169.View ArticleADSGoogle Scholar
- Hoang VV: Defect and Diffussion Forum. 2005, 242–244: 77-10.4028/www.scientific.net/DDF.242-244.77.View ArticleGoogle Scholar
- Bromley ST, Flikkema E: Comp Mater Sci. 2006, 35: 382-10.1016/j.commatsci.2004.08.018.View ArticleGoogle Scholar
- Schweigert IV, Lehtinen KEJ, Carrier MJ, Zachariah MR: Phys Rev B. 2002, 65: 235410-10.1103/PhysRevB.65.235410.View ArticleADSGoogle Scholar
- Roder A, Kob W, Binder K: J Chem Phys. 2001, 114: 7602-10.1063/1.1360257.View ArticleADSGoogle Scholar
- Song J, Choi M: Phys Rev B. 2002, 65: 241302-10.1103/PhysRevB.65.241302.View ArticleADSGoogle Scholar
- Hoang VV, Phys J: Chem B. 2007, 111: 12649-View ArticleGoogle Scholar
- Hoang VV: NANO. 2007, 2: 301-10.1142/S1793292007000659.View ArticleGoogle Scholar
- Huff NT, Demiralp E, Cagin T, Goddard WA, Non-Cryst J: Solids. 1999, 253: 133-Google Scholar
- Takada A, Richet P, Catlow CRA, Price GD, Non-Cryst J: Solids. 2004, 345–346: 224-Google Scholar
- Hoang VV, Hai NT, Zung H: Phys Lett A. 2006, 356: 246-10.1016/j.physleta.2006.03.052.View ArticleADSGoogle Scholar
- BWH van Beest, Kramer GJ, RA van Santen: Phys Rev Lett. 1990, 64: 1955-10.1103/PhysRevLett.64.1955.View ArticleADSGoogle Scholar
- Tsuneyuki S, Tsukada M, Aoki H, Matsui Y: Phys Rev Lett. 1988, 61: 869-10.1103/PhysRevLett.61.869.View ArticleADSGoogle Scholar
- Tsuneyuki S, Matsui Y: Phys Rev Lett. 1995, 74: 3197-10.1103/PhysRevLett.74.3197.View ArticleADSGoogle Scholar
- Alcoutlabi M, McKenna GB, Phys J: Condens Matter. 2005, 17: R461-10.1088/0953-8984/17/15/R01.View ArticleADSGoogle Scholar
- Adiga SP, Zapol P, Curtiss LA: Phys Rev B. 2006, 74: 064204-10.1103/PhysRevB.74.064204.View ArticleADSGoogle Scholar
- Binder K, Kob W: Glassy Materials and Disordered Solids (World Scientific, Singapore, 2005).Google Scholar
- Hoang VV, Oh SK: Physica B. 2004, 352: 73-10.1016/j.physb.2004.08.011.View ArticleADSGoogle Scholar
- Diebold U: Surf Sci Rep. 2003, 48: 53-10.1016/S0167-5729(02)00100-0.View ArticleADSGoogle Scholar
- Stesmans A, Clemer K, Afanas'es VV, Phys J: Condens Matter. 2005, 17: L393-10.1088/0953-8984/17/38/L01.View ArticleADSGoogle Scholar
- Hoang VV: Eur Phys J B. 2006, 54: 291-10.1140/epjb/e2007-00005-y.View ArticleADSGoogle Scholar
- McCoy JD, Curro JG: J Chem Phys. 2002, 116: 9154-10.1063/1.1481379.View ArticleADSGoogle Scholar
- Hoang VV: J Phys D: Appl Phys. 2007, 40: 7454-10.1088/0022-3727/40/23/029.View ArticleADSGoogle Scholar
- Hoang VV, Zung H: Amorphous-amorphous phase transition in amorphous simulated SiO2 nanoparticles, Tech. Proc. of the. 2007, Santa Clara, California, U.S.A., NSTI Nanotechnology Conference and Trade Show, May 20–24, 2007Google Scholar
- Hoang VV, Zung H, Hai NT: J Phys: Condens Matter. 2007, 19: 116104-10.1088/0953-8984/19/11/116104.ADSGoogle Scholar
Copyright
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.