(72b) Sintering Rate and Mechanism of TiO2 Nanoparticles by Molecular Dynamics Using GPUs
AIChE Annual Meeting
2011
2011 Annual Meeting
Computational Molecular Science and Engineering Forum
Beyond Standard Hardware: GPUs, Cloud Computing and Crowdsourcing
Monday, October 17, 2011 - 12:46pm to 1:02pm
Sintering Rate and Mechanism of TiO2
Nanoparticles
by Molecular Dynamics
B. Buesser, A.J. Gröhn and S.E. Pratsinis
Particle Technology Laboratory, Institute of Process Engineering,
Department of Mechanical and Process Engineering, ETH
Zürich, 8092 Zürich, Switzerland
Titania (TiO2)
nanoparticles have many attractive applications in photovoltaic1 and photocatalytic2 processes, to name a few.
The performance of TiO2 nanoparticles depends considerably on their
size and composition which are determined by the sintering rate during their synthesis.
The sintering rate is crucial for designing gas-phase processes with controlled
product particle size, structure, composition and eventual performance in a
number of applications3.
Seto et al.4 have proposed a sintering
rate similar to Kobata et al.5 for TiO2 and
validated it by accounting for the detailed fluid mechanics of their hot wall
reactor forming rather large TiO2 nanoparticles (dp
= 10 ? 100 nm). Little is known, however, for the sintering rate of small TiO2
nanoparticles (dp < 10 nm) as it is difficult to reliably
measure it.
On the other
hand molecular dynamics (MD) simulations have been used to study the sintering
to full coalescence of metallic and metalloid nanoparticles6,7 though much less has been
done for ceramic ones as their force fields and potentials are difficult to
determine. Koparde and Cummings8 have investigated the
sintering of two TiO2 nanoparticles, with CPU's reaching up to 1 ns
in residence time, by tracking the shrinkage of the center-to-center distance
and the growth of the sintering neck using the force field of Matsui and Akaogi9, while comparing
MD with phenomenological sintering models and investigating the melting of TiO2
nanoparticles10.
Here11, sintering of rutile TiO2
nanoparticles is investigated by MD on graphical processing units (GPU)12 of a consumer graphic
card that was originally developed for the visualization of highly realistic
computer games. The acceleration achieved by the GPU, compared to a single CPU
(Figure 1), allowed simulating TiO2 nanoparticle sintering from
adhesion and neck growth to finally full coalescence over several hundred
nanoseconds for the first time and determining the sintering rate of very small
TiO2 nanoparticles (dp < 5 nm) by monitoring
the evolution of their surface area (Figure 2). For the smallest particle
diameters, the MD-obtained sintering rates were smaller than that predicted by
theory developed for larger particles4,5. Ions on the particle
surface exhibited higher net displacement than bulk ones revealing that surface
diffusion is the dominant sintering mechanism of TiO2 nanoparticles.
An expression for the sintering rate of rutile TiO2 nanoparticles
has been extracted from MD, bridging the gap of knowledge from a few molecules
to several nanometers, the key size range for nanoparticle properties and
performance. This MD-derived sintering rate could facilitate the use of
phenomenological models13 to design processes for
large scale manufacture and processing of small nanoparticles14. The low costs for consumer
GPU hardware and low energy consumption combined with extremely high
computational power make GPU's suitable for scientific computing.
Financial support from the Swiss National Science
Foundation (SNF) grant # 200021-119946/1 is gratefully acknowledged.
References
1. Gratzel M. Photoelectrochemical
cells. Nature. 2001; 414, (6861), 338-344.
2. Formenti M, Juillet F, Meriaudeau
P, Teichner SJ, Vergnon P. Preparation in a hydrogen-oxygen flame of ultrafine
metal oxide particles. Oxidative properties toward hydrocarbons in the presence
of ultraviolet radiation. J. Colloid Interface Sci. 1972; 39, (1),
79-89.
3. Strobel R, Pratsinis SE. Flame
aerosol synthesis of smart nanostructured materials. J. Mater. Chem. 2007;
17, 4743-4756.
4. Seto T, Shimada M, Okuyama K.
Evaluation of sintering of nanometer-sized titania using aerosol method. Aerosol
Sci. Technol. 1995; 23, (2), 183-200.
5. Kobata A, Kusakabe K, Morooka S.
Growth and transformation of TiO2 crystallites in aerosol reactor. AIChE
J. 1991; 37, (3), 347-359.
6. Arcidiacono S, Bieri NR, Poulikakos
D, Grigoropoulos CP. On the coalescence of gold nanoparticles. Int. J.
Multiphas. Flow. 2004; 30, (7-8), 979-994.
7. Zachariah MR, Carrier MJ. Molecular
dynamics computation of gas-phase nanoparticle sintering: a comparison with
phenomenological models. J. Aerosol Sci. 1999; 30, (9), 1139-1151.
8. Koparde VN, Cummings PT. Molecular
dynamics simulation of titanium dioxide nanoparticle sintering. J. Phys.
Chem. B. 2005; 109, (51), 24280-24287.
9. Matsui M, Akaogi M. Molecular
dynamics simulation of the strucutral and physical properties of the four
polymorphs of TiO2. Mol. Simulat. 1991; 6, 239-244.
10. Koparde VN, Cummings PT. Sintering of
titanium dioxide nanoparticles: a comparison between molecular dynamics and
phenomenological modeling. J. Nanopart. Res. 2008; 10, (7), 1169-1182.
11. Buesser B, Gröhn A, Pratsinis SE.
Sintering Rate and Mechanism of TiO2 Nanoparticles by Molecular
Dynamics. J. Phys. Chem. C. 2011; accepted,
12. Garland M, Le Grand S, Nickolls J,
Anderson J, Hardwick J, Morton S, Phillips E, Yao Z, Volkov V. Parallel
Computing Experiences with CUDA. Micro IEEE. 2008; 28, (4), 13-27.
13. Koch W, Friedlander SK. The effect of
particle coalescence on the surface area of a coagulating aerosol. J.
Colloid Interface Sci. 1990; 140, (2), 419-427.
14. Pratsinis SE. Aerosol-based
technologies in nanoscale manufacturing: from functional materials to devices
through core chemical engineering. AIChE J. 2010; 56, (12), 3028-3035.
Figure 1 The Number of timesteps for molecular dynamics of TiO2
nanoparticles undergoing sintering on GPU (circles) and CPU (triangles) as
function of number of atoms.
Figure
2 Evolution of
surface area of two TiO2 nanoparticles (d0 = 3 nm)
at 1800 K.