Abstract: Electron emission properties of cathode surfaces affect considerably the electrical characteristics of glow discharges. Using a heavy-particle hybrid model ill 2 dimensions, we investigate the influence of the secondary electron emission coefficient gamma oil the calculated discharge characteristics for both 'clean' and 'dirty' cathode surface conditions, and assuming a constant gamma parameter as well. The effect of the gas heating and the role of the heavy particles reflected from the cathode on this process is also studied.

Abstract: Our aim is to develop and optimize a model for a dielectric barrier discharge used as a chemical reactor for gas treatment. In order to determine the optimum operating conditions, we have studied the influence of the gas flow rate, reactor geometry and applied voltage parameters on the discharge characteristics. For this purpose, a two-dimensional time-dependent fluid model has been applied to an atmospheric pressure DBD in helium with nitrogen impurities, in a cylindrical geometry. The numerical model is based on the continuity and flux equations for each type of particles treated, the electron energy equation and the Poisson equation. The gas flow is incorporated in the flux equations as a source term. The set of coupled partial differential equations is solved by the so-called modified strongly implicit method. The background gas flow is numerically treated separately, assuming in the model that there is no influence of the plasma on the flow. Indeed, the gas convection velocity is calculated using the commercial code Fluent and it is used as input into the 2D fluid model. The plasma characteristics have been studied in terms of gas flow rate, applied voltage amplitude and frequency, and geometrical effects. The electric currents as a function of time for a given applied potential have been obtained, as well as the number densities and fluxes of plasma species.

Abstract: The reaction behavior of species that may affect the growth of ultrananocrystal line and nanocrystalline diamond ((U)NCD) films is investigated by means of molecular dynamics simulations. Impacts of CHx (x = 0 – 4), C2Hx (x=0-6), C3Hx (x=0-2), C4Hx (x = 0 – 2), H, and H-2 on clean and hydrogenated diamond (100)2 x 1 and (111) 1 x 1 surfaces at two different substrate temperatures are simulated. We find that the different bonding structures of the two surfaces cause different temperature effects on the sticking efficiency. These results predict a temperature-dependent ratio of diamond (100) and (111) growth. Furthermore, predictions of which are the most important hydrocarbon species for (U)NCD growth are made.

Abstract: In this paper, the initial mechanisms of nanoparticle formation and growth in radiofrequency acetylene (C2H2) plasmas are investigated by means of a comprehensive self-consistent one-dimensional (1D) fluid model. This model is an extension of the 1D fluid model, developed earlier by De Bleecker et al. Based on the comparison of our previous results with available experimental data for acetylene plasmas in the literature, some new mechanisms for negative ion formation and growth are proposed. Possible routes are considered for the formation of larger (linear and branched) hydrocarbons C2nH2 (n = 3, 4, 5), which contribute to the generation of C2nH− anions (n = 3, 4, 5) due to dissociative electron attachment. Moreover, the vinylidene anion (H2CC−) and higher anions (n = 24) are found to be important plasma species.

Abstract: A hybrid model, called the hybrid plasma equipment model, was used to study Ar/Cl(2) inductively coupled plasmas used for the etching of Si. The effects of substrate bias, source power and gas pressure on the plasma characteristics and on the fluxes and energies of plasma species bombarding the substrate were observed. A comparison with experimentally measured etch rates was made to investigate how the etch process is influenced and which plasma species mainly account for the etch process. First, the general plasma characteristics are investigated at the following operating conditions: 10% Ar 90% Cl(2) gas mixture, 5mTorr total gas pressure, 100 sccm gas flow rate, 250W source power, -200V dc bias at the substrate electrode and an operating frequency of 13.56MHz applied to the coil and to the substrate electrode. Subsequently, the pressure is varied from 5 to 80mTorr, the substrate bias from -100 to -300V and the source power from 250 to 1000W. Increasing the total gas pressure results in a decrease of the etch rate and a less anisotropic flux to the substrate due to more collisions of the ions in the sheath. Increasing the substrate bias has an effect on the energy of the ions bombarding the substrate and to a lesser extent on the magnitude of the ion flux. When source power is increased, it was found that, not the energy, but the magnitude of the ion flux is increased. The etch rate was more influenced by a variation of the substrate bias than by a variation of the source power, at these operating conditions. These results suggest that the etch process is mainly affected by the energy of the ions bombarding the substrate and the magnitude of the ion flux, and to a lesser extent by the magnitude of the radical flux.

Abstract: This paper presents the computational study of the metal-organic (MO) CVD of titanium dioxide (TiO2) films grown using titanium tetraisopropoxide (TTIP) as a precursor and nitrogen as a carrier gas. The TiO2 films are deposited under atmospheric pressure. The effects of the precursor concentration, the substrate temperature, and the hydrolysis reaction on the deposition process are investigated. It is found that hydrolysis of the TTIP decreases the onset temperature of the gas-phase thermal decomposition, and that the deposition rate increases with the precursor concentration and with the decrease of substrate temperature. Concerning the mechanism responsible for the film growth, the model shows that at the lowest precursor concentration, the direct adsorption of the precursor is dominant, while at higher precursor concentrations, the monomer deposition becomes more important.

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: 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.

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: Optical (atomic absorption spectroscopy, AAS; atomic emission spectroscopy, AES; atomic fluorescence spectroscopy, AFS; and optogalvanic spectroscopy) and mass spectrometric (magnetic sector, quadrupolemass analyzer, QMA; quadrupole ion trap, QIT; Fourier transform ion cyclotron resonance, FTICR; and time-of-flight, TOF) instrumentation are well suited for coupling to the glow discharge (GD). The GD is a relatively simple device. A potential gradient (5001500 V) is applied between an anode and a cathode. In most cases, the sample is also the cathode. A noble gas (e.g. Ar, Ne, and Xe) is introduced into the discharge region before power initiation. When a potential is applied, electrons are accelerated toward the anode. As these electrons accelerate, they collide with gas atoms. A fraction of these collisions are of sufficient energy to remove an electron from a support gas atom, forming an ion. These ions are, in turn, accelerated toward the cathode. These ions impinge on the surface of the cathode, sputtering sample atoms from the surface. Sputtered atoms that do not redeposit on the surface diffuse into the excitation/ionization regions of the plasma where they can undergo excitation and/or ionization via a number of collisional processes. GD sources offer a number of distinct advantages that make them well suited for specific types of analyses. These sources afford direct analysis of solid samples, thus minimizing the sample preparation required for analysis. The nature of the plasma also provides mutually exclusive atomization and excitation processes that help to minimize the matrix effects that plague so many other elemental techniques. Unfortunately, the GD source functions optimally in a dry environment, making analysis of solutions more difficult. These sources also suffer from difficulties associated with analyzing nonconductingsamples. In this article, first, the principles of operation of the GD plasma are reviewed, with an emphasis on how those principles relate to optical spectroscopy and mass spectrometry. Basic applications of the GD techniques are considered next. These include bulk analysis, surface analysis, and the analysis of solution samples. The requirements necessary to obtain optical information are addressed following the analytical applications. This section focuses on the instrumentation needed to make optical measurements using the GD as an atomization/excitation source. Finally, mass spectrometric instrumentation and interfaces are addressed as they pertain to the use of a GD plasma as an ion source. GDsources provide analytically useful gas-phase species from solid samples. These sources can be interfaced with avariety of spectroscopic and spectrometric instruments for both quantitative and qualitative analysis.

Abstract: In this paper, a new implementation of the Metropolis Monte Carlo (MMC) algorithm is presented. When combining the MMC model with a molecular dynamics (MD) code, crystal growth by plasma-enhanced chemical vapor deposition can be simulated. As the MD part simulates impacts of growth species onto the surface on a time scale of picoseconds, the MMC algorithm simulates the slower adatom surface processes. The implementation includes a criterion for the selection of atoms that are allowed to be displaced during the simulation, and a criterion of after how many MMC cycles the simulation is stopped. We performed combined MD-MMC simulations for hydrocarbon species that are important for the growth of ultrananocrystalline diamond (UNCD) films at partially hydrogenated diamond surfaces, since this implementation is part of a study of the growth mechanisms of (ultra)nanocrystalline diamond films. Exemplary for adatom arrangements during the growth of UNCD, the adatom surface behavior of C and C2H2 at diamond (111)1 × 1, C and C4H2 at diamond (111)1 × 1 and C3 at diamond (100)2 × 1 has been investigated. For all cases, the diamond crystal structure is pursued under the influence of MMC simulation. Additional longer time-scale MD simulations put forward very similar structures, verifying the MMC algorithm. Nevertheless, the MMC simulation time is typically one order of magnitude shorter than the MD simulation time.

Abstract: We present here some of our modeling efforts for PECVD growth of nanostructured carbon materials with focus on amorphous hydrogenated carbon. Experimental data from an expanding thermal plasma setup were used as input for the simulations. Attention was focused both on the film growth mechanism, as well as on the hydrocarbon reaction mechanisms during growth of the films. It is found that the reaction mechanisms and sticking coefficients are dependent on the specific surface sites, and the structural properties of the growth radicals. The film growth results are in correspondence with the experiment. Furthermore, it is found that thin a-C:H films can be densified using an additional H-flux towards the substrate.

Abstract: Molecular dynamics simulations were used to investigate the size-dependent melting mechanism of nickel nanoclusters of various sizes. The melting process was monitored by the caloric curve, the overall cluster Lindemann index, and the atomic Lindemann index. Size-dependent melting temperatures were determined, and the correct linear dependence on inverse diameter was recovered. We found that the melting mechanism gradually changes from dynamic coexistence melting to surface melting with increasing cluster size. These findings are of importance in better understanding carbon nanotube growth by catalytic chemical vapor deposition as the phase state of the catalyst nanoparticle codetermines the growth mechanism.