Atomistic modelling of Si nanoparticles synthesis

Silicon remains the most important material for electronic technology. Presently, some efforts are focused on the use of Si nanoparticles--not only for saving material, but also for improving the efficiency of optical and electronic devices, for instance, in the case of solar cells coated with a film of Si nanoparticles. The synthesis by a bottom-up approach based on condensation from low temperature plasma is a promising technique for the massive production of such nanoparticles, but the knowledge of the basic processes occurring at the atomistic level is still very limited. In this perspective, numerical simulations can provide fundamental information of the nucleation and growth mechanisms ruling the bottsom-up formation of Si nanoclusters. We propose to model the low temperature plasma by classical molecular dynamics by using the reactive force field (ReaxFF) proposed by van Duin, which can properly describe bond forming and breaking. In our approach, first-principles quantum calculations are used on a set of small Si clusters in order to collect all the necessary energetic and structural information to optimize the parameters of the reactive force-field for the present application. We describe in detail the procedure used for the determination of the force field and the following molecular dynamics simulations of model systems of Si gas at temperatures in the range 2000-3000 K. The results of the dynamics provide valuable information on nucleation rate, nanoparticle size distribution, and growth rate that are the basic quantities for developing a following mesoscale model.