Abstract: In this paper, an overview is given of modeling activities going on in our research group, for describing the plasma chemistry and plasmasurface interactions in reactive plasmas. The plasma chemistry is calculated by a fluid approach or by hybrid Monte Carlo (MC)fluid modeling. An example of both is illustrated in the first part of the paper. The example of fluid modeling is given for a dielectric barrier discharge (DBD) in CH4/O2, to describe the partial oxidation of CH4 into value-added chemicals. The example of hybrid MCfluid modeling concerns an inductively coupled plasma (ICP) etch reactor in Ar/Cl2/O2, including also the description of the etch process. The second part of the paper deals with the treatment of plasmasurface interactions on the atomic level, with molecular dynamics (MD) simulations or a combination of MD and MC simulations.

Abstract: In this letter, we investigate the effect of various levels of nitrogen impurity on the electrical performance of an atmospheric pressure dielectric barrier discharge in helium. We illustrate the different current profiles that are obtained, which exhibit one or more discharge pulses per half cycle and evaluate their performance in ionizing the discharge and dissipating the power. It is shown that flat and broad current profiles perform the best in ionizing the discharge and use the least amount of power per generated charged particle.

Abstract: In this theoretical study, the electron emission from a copper surface under ultrashort pulsed laser irradiation is investigated using a one-dimensional particle in cell model. Thermionic emission as well as multiphoton photoelectron emission were taken into account. The emitted electrons create a negative space charge above the target; consequently the generated electric field reduces the electron emission by several orders of magnitude. The simulations indicate that the space charge effect should be considered when investigating electron emission related phenomena in materials under ultrashort pulsed laser irradiation of metals.

Abstract: The influence of the composition on the crystallographic properties of deposited Mg(M)O (with M=Al, Cr, Ti, Y, and Zr) films is studied. For a flexible control of the composition, dual reactive magnetron sputtering was used as deposition technique. Two different approaches to predict the composition are discussed. The first is an experimental way based on the simple relationship between the deposition rate and the target-substrate distance. The second is a route using a Monte Carlo based particle trajectory code. Both methods require a minimal experimental input and enable the user to quickly predict the composition of complex thin films. Good control and flexibility allow us to study the compositional effects on the growth of Mg(M)O films. Pure MgO thin films were grown with a (111) preferential out-of-plane orientation. When adding M to MgO, two trends were noticed. The first trend is a change in the MgO lattice parameters compared to pure MgO. The second tendency is a decrease in the crystallinity of the MgO phase. The experimentally determined crystallographic properties are shown to be in correspondence with the predicted properties from molecular dynamics simulations.

Abstract: The phase transition processes induced by ultrashort, 100 fs pulsed laser irradiation of Au, Cu, and Ni are studied by means of a combined atomistic-continuum approach. A moderately low absorbed laser fluence range, from 200 to 600 J/m2 is considered to study phase transitions by means of a local and a nonlocal order parameter. At low laser fluences, the occurrence of layer-by-layer evaporation has been observed, which suggests a direct solid to vapor transition. The calculated amount of molten material remains very limited under the conditions studied, especially for Ni. Therefore, our results show that a kinetic equation that describes a direct solid to vapor transition might be the best approach to model laser-induced phase transitions by continuum models. Furthermore, the results provide more insight into the applicability of analytical superheating theories that were implemented in continuum models and help the understanding of nonequilibrium phase transitions.

Abstract: Mesoporous carbon, which can be replicated from mesoporous silica and whose surface is hydrophobic, can be an ideal catalyst for the esterification of fatty acids. Here we report an easy and low cost way to prepare sulfonic acid group-functionalized mesoporous carbon. A sample of calcined mesoporous silica SBA-15 was added to an aqueous sucrose solution followed by drying and calcination at different temperatures. In contrast to existing procedures, the obtained hybrid Si/C material was then first sulfonated in H2SO4, before the final removal of the silica template in order to stabilize the porous structure towards the liquid phase sulfonation treatment. Thus the silicacarbon composites, instead of the mesoporous carbon, were successfully sulfonated to introduce SO3H groups, while keeping the ordered mesoporous structure intact. The influence of carbonization temperature was investigated, suggesting an optimum temperature of 873 K. The SO3H group-functionalized mesoporous carbon, denoted as CMK-3-873-SO3H, was characterized by means of XRD, N2 physisorption, SEM, FT-IR, elemental analysis and TEM. It followed that a uniform mesoporous carbon was obtained with an average pore size of 3.89 nm, a specific surface of 807 m2/g and a SO3H group loading of 0.39 meq/g of dry material. Compared with other solid acid catalysts, the resulting material shows enhanced activity in the acid-catalyzed esterification of oleic acid with methanol, and can be used repeatedly. The increased catalytic performance is attributed to the hydrophobic surface and larger pore size of the new catalyst. It can effectively accommodate long chain fatty acids and reject formed water, making the active sites easily accessible.

Abstract: Nanocrystalline diamond films have generated much interested due to their diamond-like properties and low surface roughness. Several techniques have been used to obtain a high re-nucleation rate, such as hydrogen poor or high methane concentration plasmas. In this work, the properties of nano-diamond films grown on silicon substrates using a continuous DC bias voltage during the complete duration of growth are studied. Subsequently, the layers were characterised by several morphological, structural and optical techniques. Besides a thorough investigation of the surface structure, using SEM and AFM, special attention was paid to the bulk structure of the films. The application of FTIR, XRD, multi wavelength Raman spectroscopy, TEM and EELS yielded a detailed insight in important properties such as the amount of crystallinity, the hydrogen content and grain size. Although these films are smooth, they are under a considerable compressive stress. FTIR spectroscopy points to a high hydrogen content in the films, while Raman and EELS indicate a high concentration of sp2 carbon. TEM and EELS show that these films consist of diamond nano-grains mixed with an amorphous sp2 bonded carbon, these results are consistent with the XRD and UV Raman spectroscopy data.

Abstract: This paper outlines two different numerical simulation approaches, carried out by our group, used for describing hydrocarbon plasmas in their applications for either nanoparticle formation in the plasma or the growth of nanostructured thin films, such as nanocrystalline diamond (NCD). A plasma model based on the fluid approach is utilized to study the initial mechanisms giving rise to nanoparticle formation in an acetylene plasma. The growth of NCD is investigated by molecular dynamics simulations, describing the interaction of the hydrocarbon species with a substrate.

Abstract: We have investigated the plasma processes in an atmospheric pressure glow discharge (APGD) in He used for analytical spectrometry by means of fluid and Monte Carlo (MC) simulations. Typical results include the potential and electric field distributions in the plasma, the density profiles of the various plasma species throughout the discharge, the mean electron energy, as well as the rates of the various collision processes in the plasma, and the relative importance of the different production and loss rates for the various species. The similarities and differences with low-pressure glow discharges are discussed. The main differences are a very small cathode dark space region and a large positive column as well as the dominant role of molecular ions. Some characteristic features of the APGD, such as the occurrence of the different spatial zones in the discharge, are illustrated, with links to experimental observations.

Abstract: A numerical parameter study has been performed for a cylindrical atmospheric pressure dielectric barrier discharge (DBD) in helium with nitrogen impurities using a two-dimensional time-dependent fluid model. The calculated electric currents and gap voltages as a function of time for a given applied potential are presented, as well as the number densities of the various plasma species. This study shows that for the geometry under consideration the applied voltage parameters have a large impact on the electric current profiles and that the discharge current is always determined by the electron and ion conduction currents while the displacement current is nearly negligible. A relative broadening of the current profiles (compared with the duration of the half cycle of the applied voltage) with an increase in the applied frequency is obtained. Nearly sinusoidal current wave forms, usually typical for radio frequency DBDs, are observed while still operating at the frequencies of tens of kilohertz. For the setup under investigation, the Townsend mode of the DBD is observed in the entire range of applied voltage amplitudes and frequencies. It is shown that the average power density dissipated in the discharge increases with rising applied voltage and frequency. An increase in applied voltage frequency leads to an increase in the electron density and a decrease in electron energy, while increasing the voltage amplitude has the opposite effect.

Abstract: In this paper, some modelling approaches to describe direct current (dc) magnetron discharges developed in our research groups will be presented, including an analytical model, Monte Carlo simulations for the electrons and for the sputtered atoms, a hybrid Monte Carlo-fluid model and particle-in-cell-Monte Carlo collision simulations. The strengths and limitations of the various modelling approaches will be explained, and some characteristic simulation results will be illustrated. Furthermore, some other simulation methods related to the magnetron device will be briefly explained, more specifically for calculating the magnetic field distribution inside the discharge, and for describing the (reactive) sputtering.

Abstract: The nucleation and growth of titanium dioxide (TiO2) and zinc oxide (ZnO) thin films on Fe2O3 (hematite), Al2O3 (á-alumina) and SiO2 (á-quartz) are studied by molecular dynamics simulations. The results show the formation of a strong interface region between the substrate and the film in the six systems studied here. A combination of polycrystalline and amorphous phases are observed in the TiO2 films grown on the three substrates. ZnO deposition on the Fe2O3 and Al2O3 crystals yields a monocrystalline film growth. The ZnO film deposited on the SiO2 crystal exhibits less crystallinity. The simulation results are compared with experimental results available in the literature.

Abstract: The atmospheric pressure dielectric barrier discharge in helium is a pulsed discharge in nature. If during the electrical current pulse a glow discharge is reached, then this pulse will last only a few microseconds in operating periods of sinusoidal voltage with lengths of about 10 to 100 µs. In this paper we demonstrate that right before a glow discharge is reached, the discharge very closely resembles the commonly assumed Townsend discharge structure, but actually contains some significant differing features and hence should not be considered as a Townsend discharge. In order to clarify this, we present calculation results of high time and space resolution of the pulse formation. The results indicate that indeed a maximum of ionization is formed at the anode, but that the level of ionization remains high and that the electric field at that time is significantly disturbed. Our results also show where this intermediate structure comes from.

Abstract: In this paper, an overview is given of different modeling approaches used for describing gas discharge plasmas, as well as plasma-surface interactions. A fluid model is illustrated for describing the detailed plasma chemistry in capacitively coupled rf discharges. The strengths and limitations of Monte Carlo simulations and of a particle-in-cell-Monte Carlo collisions model are explained for a magnetron discharge, whereas the capabilities of a hybrid Monte Carlo-fluid approach are illustrated for a direct current glow discharge used for spectrochemical analysis of materials. Finally, some examples of molecular dynamics simulations, for the purpose of plasma-deposition, are given.